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Abstract 



Many stochastic differential equations (SDEs) in the literature have a superlinearly growing nonlinearity in 
their drift or diffusion coefficient. Unfortunately, moments of the computationally efficient Euler-Maruyama 
approximation method diverge for these SDEs in finite time. This article develops a general theory for 
^^ , studying integrability properties such as moment bounds for discrete-time stochastic processes. Using this 

p^ • approach, we establish moment bounds for fully and partially drift-implicit Euler methods and for a class 

• I of new explicit approximation methods which require only a few more arithmetical operations than the 
1 -^ ■ Euler-Maruyama method. These moment bounds are then used to prove strong convergence of the proposed 

j^ ' schemes. Finally, we illustrate our results for several SDEs from finance, physics, biology and chemistry. 
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1 Introduction 

This article investigates integrability and convergence properties of numerical approximation processes for 
stochastic differential equations (SDEs). In order to illustrate one of our main results, the following general 
setting is considered in this introductory section. Let T e (0, cx)), d, ?7i € N := {1, 2, . . .}, let (fi, J-, P) be a prob- 
ability space with a normal filtration {Ft)t£[o,T\i let W : [0,T] x 57 -^- R™ be a standard (Jt)tg[o,T]-Brownian mo- 
tion, let D C M'* be an open set, let /i = (^i, . . . , ^d) : D ^)-W^ and a = {(Ji,j)ie{i,2....,d}.j£{i,2,...,m} ■ D -> R''><'" 
be locally Lipschitz continuous functions and let X : [0, T] x fi ^ £) be an (J^()tg[o_T]-adapted stochastic process 
with continuous sample paths satisfying the SDE 

Xt^Xo+ [ KXs) ds + [ c7{Xs) dW, (1) 

P- almost surely for all t G [0,r]. Here /i is the infinitesimal mean and a ■ a* is the infinitesimal covariance 
matrix of the solution process X of the SDE ^. To guarantee finiteness of some moments of the SDE ((TJ, we 
assume existence of a Lyapunov-type function. More precisely, let q € (0,oo), k e R be real numbers and let 
V: D ^ [l,oo) be a twice continuously differentiable function with E[y(Xo)] < oo and with V{x) > ||x||'' and 






E ( i^ ) (^) ■ ^»(^) + ^ > ^ > ^ ( ^zi^ ) (^) • ^^fc(^) • '^j^fe(^) ^ '^ ■ ^(^) (2) 



for all X £ D. These assumptions ensure E[l/(Xt)] < e*^* • E[V^(Xo)] for all t E [0,T] and, therefore, finiteness of 
the g-th absolute moments of the solution process Xt, t e [0,r], of the SDE (P), i.e., supfg[o,T] '^[ll-'^tll''] < '^■ 
Note that in this setting both the drift coefficient n and the diffusion coefficient a of the SDE ^ may grow 
superlinearly and are, in particular, not assumed to be globally Lipschitz continuous. Our main goal in this 
introduction is to construct and to analyze numerical approximation processes that converge strongly to the 
exact solution of the SDE ([IJ . Strong numerical approximations of the SDE ^ are of particular interest for the 
computation of statistical quantities of the solution process of the SDE dH) through computationally efficient 
multilevel Monte Carlo methods (see Giles [H], Heinrich [53] and, e.g., in Creutzig et al. [2], Hickernell et 
al. [27], Barth, Lang & Schwab [5; and in the references therein for further recent results on multilevel Monte 
Carlo methods). 

Several SDEs from the literature satisfy the above setting (see Subsections l4.2fT4.11l below). For instance, 
the function V{x) = (l + ||a;|p) , x E D, for an arbitrary r E (0,oo) serves as a Lyapunov-type function for 
the stochastic van der Pol oscillator (I269p . for the stochastic Lorenz equation (|275p . for the Cox- IngersoU- Ross 
process (j307p and for the simplified Ait-Sahalia interest rate model (|308D but not for the stochastic Duffing-van 
dcr Pol oscillator (|272l) , not for the stochastic Brusselator (j277p , not for the stochastic SIR model (I283P , not for 
the Lotka-Volterra predator prey model (|303p and, in general, also not for the Langevin equation (I310p . The 
function V{x) = (l + (xi)"^ + 2{x2)'^) , x = (a;i,a;2) E D = R^, for an arbitrary r E (0,oo) is a Lyapunov- 
type function for the stochastic Duffing-van der Pol oscillator (|272p . For the stochastic SIR model (|283l) . 
the function V{x) = (l -I- (xi + X2)'^ + (xs)"^) , x = {xi,X2,X3) E D = (0,oo)^, for an arbitrary r E (0, oo) 
serves as a Lyapunov-type function and for the stochastic Lotka-Volterra system p94p . the function V{x) = 
(l + {vix\ + . . . VdXdY^ , X E D = {0, oo)'*, for an arbitrary r E (0, oo) and an appropriate v = {vi , . . . , Vd) G K'* 
is a Lyapunov-type function. More details on the examples can be found in Section [D 

The standard method for approximating SDEs with globally Lipschitz continuous coefficients is the Euler- 
Maruyama method. Unfortunately, the Euler-Maruyama method often fails to converge strongly to the ex- 
act solution of nonlinear SDEs of the form ((IJ; see [35]. Indeed, if at least one of the coefficients of the 
SDE grows superlinearly, then the Euler-Maruyama scheme diverges in the strong sense. More precisely, let 
Z^ : {0, 1, ... , N} X fl — > R'', TV e N, be Euler-Maruyama approximations for the SDE ([T|) defined recursively 
through Zq := Xq and 



zN r^N 



1 := Z^ + M(Z,f )^ + a{Z:!){W,,,^ - W^) (3) 

for all n e {0, 1, . . . , iV - 1} and all N E N. Here p.: W'- ^ R'* and ct: R"* ^ M'*^" are extensions of fi. and 
a given by fl{x) = 0, a-{x) = for all x E D"^ and by fl{x) = fi{x), a{x) = a{x) for all x E D respectively 
(see also Subsections 13.21 and 13.61 for more general extensions). Theorem 2.1 of [35] (which generalizes Theorem 
2.1 of [35]) then implies in the case d — m — 1 that if there exists a real number e E (0, oo) such that 
\jl{x)\ + \a{x)\ > £|x|(i+^) for aU |j;| > 1/e and if f[a{Xa) ^ 0] > 0, then limAr^oo E[|y^|'^] = oo for all 
r E (0, oo) and therefore limAr-j-oo IE[|Xt — i^^vT] — °° ^^^ ^^^ ''' ^ C^'^l (^^^ ^1^° Sections 4 and 5 in [3S] for 
divergence results for the corresponding multilevel Monte Carlo Euler method). Due to these deficiencies of the 
Euler-Maruyama method, we look for numerical approximation methods whose computational cost is close to 
the Euler-Maruyama method and which converge strongly even in the case of SDEs with superlinearly growing 
coefficients. 



There are a number of strong convergence results for temporal numerical approximations of SDEs of the 
form ([T]) with possibly superlinearly growing coefficients in the literature. Many of these results assume beside 
other assumptions that the drift coefficient fi of the SDE ([1]) is globally one-sided Lipschitz continuous and 
also prove rates of convergence in that case. In particular, if the drift coefficient is globally one-sided Lipschitz 
continuous and if the diffusion coefficient is globally Lipschitz continuous beside other assumptions, then strong 
convergence of the fully drift-implicit Euler method follows from Theorem 2.4 in Hu |:32j and from Theorem 
5.3 in Higham, Mao & Stuart 30 , strong convergence of the split-step backward Euler method follows from 
Theorem 3.3 in Higham, Mao & Stuart [30J, strong convergence of a drift-tamed Euler-Maruyama method 
follows from Theorem 1.1 in |34] and strong convergence of a drift-tamed Milstein scheme follows from Theorem 
3.2 in Gan & Wang [TS]. Theorem 2 and Theorem 3 in Higham & Kloeden [35] generalize Theorem 3.3 and 
Theorem 5.3 in Higham, Mao & Stuart 30 to SDEs with Poisson-driven jumps. In addition. Theorem 6.2 in 
Szpruch et al. [7Tj establishes strong convergence of the fully drift-implicit Euler method of a one-dimensional 
Ait-Sahalia-type interest rate model having a superlinearly growing diffusion coefficient a and a globally one- 
sided Lipschitz continuous drift coefficient fi which is unbounded near 0. Moreover, Theorem 4.4 in Szpruch 
& Mao liOj generalizes this result to a class of SDEs which have globally one-sided Lipschitz continuous drift 
coefficients and in which the function V{x) = 1-1- ||x|p, a; € Z?, is a Lyapunov-type function. A similar method 
is used in Proposition 3.3 in Dereich, Neuenkirch & Szpruch [16 to obtain strong convergence of a drift-implicit 
Euler method for a class of Bessel type processes. Moreover, Gyongy & Millet establish in Theorem 2.10 in 
|23| strong convergence of implicit numerical approximation processes for a class of possibly infinite dimensional 
SDEs whose drift fi and diffusion a satisfy a suitable one-sided Lipschtz condition (see Assumption (CI) in 
[23 for details). Strong convergence of numerical approximations for two-dimensional stochastic Navier-Stokes 
equations are obtained in Theorem 7.1 in Brzezniak, Carelli & Prohl [9] and in Corollary 3.2 in Dorsek [15]. In 
all of the above mentioned results from the literature, the function V{x) = 1-1- Hxp, x G D, is a Lyapunov-type 
function of the considered SDE. A result on more general Lyapunov-type functions is the framework in Schurz |69j 
which assumes general abstract conditions on the numerical approximations. The applicability of this framework 
is demonstrated in the case of SDEs which have globally one-sided Lipschitz continuous drift coefficients and in 
which the function V{x) — 1 + ||x|p, x G D, is a Lyapunov-type function; see [67l |68| \6^. To the best of our 
knowledge, no strong numerical approximation results are known for the stochastic van der Pol oscillator (J269I) . 
for the stochastic Duffing- van der Pol oscillator (|272l) , for the stochastic Lorenz equation (12751) , for the stochastic 
Brusselator (|277p , for the stochastic SIR model (I283P , for the experimental psychology model (I290p , for the Lotka- 
Volterra predator-prey model (j294p and for the instantaneous variance process (|309p in the Lewis stochastic 
volatility model. 

In this article, the following increment-tamed Euler-Maruyama scheme is proposed to approximate the so- 
lution process of the SDE ^ in the strong sense. Let Y'^ : {0,1 . . . ,N} x V, ^ R"^, A^ e N, be numerical 
approximation processes defined through Yf^ :— Xq and 

"+'■ " nrax(l,^||/2(>;f)i+a(r,f)(H^(,^-W^^)||) ^^ 

for all n £ {0, 1, . . . , A — 1} and all A e N. Note that the computation of ^ requires only a few additional 
arithmetical operations when compared to the computation of the Euler-Maruyama approximations ([3]) . More- 
over, we add that the scheme (|3]) is a special case of a more general class of suitable increment-tamed schemes 
proposed in Subsection 13.6.31 below. Next let F^ : [0,T] x il — > M'', A e N, be linearly interpolated continuous- 
time versions of © defined through f;^ := (n + 1 - ^)Y^ + (t " ")^ri+i ^r aU t e [nT/A, (n + 1)T/A], 
n € {0, 1, . . . , A — 1} and all A G N. For proving strong convergence of the numerical approximation processes 
Y^ , A G N, to the exact solution process X of the SDE ((Ij, we additionally assume that there exist real 
numbers 70, 71, c G [0, 00), p G [3, 00) and a three times continuously differentiable extension V : M'' — >■ [1, 00) of 
V: D ^[1,00) such that V{x) > Wxp and 

||/2(x)|| <c\V{x)\i'^], M^)\\nm-M'') <c\V{x)\V-^], \\V^'\^)\\Li^HM^M) < c\V{x)\i''i] 

for all a; G K'^ and all i G {1,2,3}. These assumptions are satisfied in all of the example SDEs from Subsec- 
tions I4.2H4.9I In the case of the squared volatility process (I304p in Subsection 14.101 and in case of the Langevin 
equation (j310p in Subsection 14.111 these assumptions are also satisfied if the model parameters satisfy suitable 
regularity conditions (see Subsections 14.101 and 14 . 1 11 for details). Under these assumptions. Theorem 13.151 below 
shows that 

lim sup E[||Xt-yt^|rl =0 (5) 

for all r G (0,g) satisfying r < -^ — ^_^3 yH) ~ \- Theorem 13.151 therebv proves strong convergence of the 

increment-tamed Euler-Maruyama method (JU) for all example SDEs from Subsections I4.2H4.9I and in parts also 



for the example SDEs from Subsections 14. 10H4. Ill Moreover, using a whole family of Lyapunov-type functions, 
we will deduce from Theorem 13.151 for most of the examples of Section 2] that strong ^''-convergence (O holds 
for all r g (0, oo) (see Corollarv 13 . 1 71 below for details). To the best of our knowledge, Theorem l3.15l is the first 
result in the literature that proves strong convergence of a numerical approximation method for the stochastic 



van der Pol oscillator (I269p . for the stochastic Duffing- van der Pol oscillator (I272p . for the stochastic Lorenz 
equation (j275p . for the stochastic Brusselator (j277p . for the stochastic SIR model (I283p . for the experimental 
psychology model (|290p . for the stochastic Lotka-Volterra predator-prey model (|303p and in the case q > 3 also 
for the instantaneous variance process (I309P in the Lewis stochastic volatility model. 

Theorem 13.151 proves the strong convergence (O in a quite general setting. One may ask whether it is also 
possible to establish a positive strong convergence rate in this setting. There is a strong hint that this is not 
possible in this general setting. More precisely, a future publication will show that in this setting there exist 
SDEs with smooth coefficients whose solution processes are not locally Holder continuous in the strong mean 
square sense with respect to the initial values. This instability suggests that for such SDEs there exist no one- 
step numerical approximation processes which converge in the strong sense with a positive strong convergence 
rate. It remains an open question which conditions on the coefficients /i (more general than globally one-sided 
Lipschitz continuous) and a of the SDE ([1]) are sufficient to ensure strong convergence of appropriate one-step 
numerical approximation processes to the exact solution of the SDE ([1} with the standard strong convergence 
order 1/2 at least. 

Finally, we summarize a few more results of this article. In Section [51 we establish uniform moment bounds 
of approximation processes for SDEs which are typically the first step in proving strong and numerically weak 
convergence results. In particular, Corollary 12.211 in Subsection 12.2.31 proves uniform moment bounds for the 
increment-tamed Euler-Maruyama scheme (j4|). Moreover, Corollarv 12 . 2 71 in Subsection 12 .3 . l1 vields uniform mo- 
ment bounds for the fully drift-implicit Euler scheme and Lemma [2.281 in Subsection 12.3.21 establishes uniform 
moment bounds for partially drift-implicit approximation schemes. These results on uniform moment bounds 
are applications of a general theory which we develop in Subsection l2.1l In this theory (see Propositions 12 . 1 1 and 
12. 71 and Corollaries l2.2ir2.3l and r2.6p we assume a Lyapunov-type inequality to be satisfied by the approximation 
processes on large subevents of the probability space; see inequality (J19p in Corollarv l2.2l One of our main results 
(Theorem 12.131 in Subsection 12.2. ip establishes this Lyapunov-type condition for the Euler-Maruyama approxi- 
mations ([3]) . More precisely, whereas the Euler-Maruyama approximations often do not satisfy a Lyapunov-type 
inequality on events of probability one in the case of superlinearly growing coefficients according to Corollarv l2.17l 
in Subsection 12.2.11 the Euler-Maruyama approximations do satisfy the Lyapunov-type inequality (J19p on large 
subevents of the probability space according to Theorem l2.13l in Subsection l2.2.1l This integrability result on the 
Euler-Maruyama approximation processes can then be transfered to a large class of other one-step approximation 
processes. More precisely. Lemma 12.181 in Subsection 12.2.21 proves that if two general one-step approximation 
schemes are close to each other in the sense of ([55]) (see Lemma 12.181 for the details) and if one approximation 
scheme satisfies the Lyapunov-type inequality (|19p on large subevents, then the other approximation scheme 
satisfies the Lyapunov-type inequality (J19p on large subevents as well. In Section [31 we then proceed to study 
convergence in probability, strong convergence and weak convergence of approximation processes for SDEs. Def- 
inition [XT] in Subsection 13 . 21 specifies a local consistency condition on approximation schemes which is, according 
to Theorem 13.31 in Subsection 13. 3[ sufficient for convergence in probability of the approximation processes to 
the exact solution of the SDE ([1} . This convergence in probability and the uniform moment bounds in Corol- 
lary [5^^ then result in the strong convergence ^ of the increment-tamed Euler-Maruyama approximations ^; 
see Theorem l3.15l in Subsection l3.4.3l for the details. Moreover, we obtain results for approximating moments and 
more general statistical quantities of the exact solution of SDEs of the form ^ in Subsection 13.51 In particular, 
Corollarv l3.23l in Subsection 13. 5. 21 establishes convergence of the Monte Carlo Euler approximations for SDEs of 
the form ([1]). 

1.1 Notation 

Throughout this article, the following notation is used. For a set fl, a measurable space {E,£) and a mapping 
Y : ^ ^ E we denote by (Ja{Y) := {Y~^{A) C il: A G £} the smallest sigma algebra with respect to which 
Y : J7 ^- _E is measurable. Moreover, let Xp G [0,00)7 P G [IjOo), be a family of real numbers such that 






LP(n;K) 

for all martingales Z; N x il — > R, all probability spaces (J7,J-", P) and all p € [l,oo). There indeed exist 
such finite real numbers Xp G [OjCso), p S [l,oo), due to the Burkholder-Davis-Gundy inequality (see, e.g.. 
Theorem 48 in Protter ^1]). Next for two sets A and B we denote by Ai{A, B) the set of all mappings from A 
to B. Furthermore, for natural numbers d, to G N and a d x m-matrix A G R'^x™ ^e denote by A* E R'^xf* the 



transpose of the matrix A. In addition, for d,m E N and arbitrary functions /i ~ (/ii, . . . , fid) : M'^ — > M'^ and 
'^^ Kjk{i,2,...,d},,G{i,2,...,rn} = K)fcg{i,2,...,™} : K'' ^ M'^^" we denotc by ^^^^ : C^R'',R) ^ M{W\R) and 
g^,^a- C2(R'^,R) ^ A^(M'' X M'',R) linear operators defined through 






2 

{fi{x), {V(p){x)) + - trace ((T(a;)cr(a;)* (Hess <y5)(a;)) (7) 



d 



and 






(^M,'^'y5)(2^,2/) := f'{x) n{y) + -^ip"{x){(Tk{y),<7k{y)) 

k=l 



(x) ■ (Ji^kix) ■ CTj^kix) 









(8) 
(a;) -o^^kiv) -OjAv) 



for alls, 2/ G E'^andallfyse C2(M'^,M) where afc : W^ ^ M'^, fc e {1, 2, . . . , m}, fulfinCTfc(2:) = (cri,fe(a:), . . . , crd,fc(a;)) 
for all X S K'^ and all k E {1, 2, . . . , m}. Furthermore, for d G N and a Borel measurable set A E B{M.'^) we 
denote by A^ : B{A) -> [0, oo] the Lebesgue-Borel measure on yl C M'^. In addition, for n,d eN, p E (0, oo] and 
a set A C R we denote by C^(R'', A) the set 

C;{R'^,A) 

j-(n~i) jg locally Lipschitz continuous and there exists a real "j (q\ 

f E C""^(M'', A) : number c e [0, cx)) such that for Agd-almost aU a: G M'' and > 
alHG{l,2,...,n}wehave||/«(x)||i(.)(K.,K)<c|/(a:)|[i-Vpl J 

throughout the rest of this article. Note that this definition is well-defined since Rademacher's theorem proves 
that a locally Lipschitz continuous function is almost everywhere differentiable. Below we also often calculate and 
formulate expressions in the extended positive real numbers [0,oo] = [0, oo) U {oo}. For instance, we frequently 
use the conventions ^ = for all a E [0, oo), ^ = oo for all a E (0, oo] and • oo = 0. 

2 Integrability properties of approximation processes for SDEs 

A central step in establishing strong and numerically weak convergence of approximation processes is to prove 
uniform moment bounds. For this analysis, we propose a Lyapunov-type condition on the one-step function of 
a one-step approximation scheme (see Definition 12.81 in Subsection I2.1.4|) . In Subsections 12.21 and 12.31 we will 
show that many numerical approximation schemes including the Euler-Maruyama scheme, the increment-tamed 
Euler-Maruyama scheme ^ and some implicit approximation schemes satisfy this condition in the case of several 
nonlinear SDEs. Subject of Subsection 12 . 1 1 is to infer from this Lyapunov-type condition on the one-step function 
that the associated numerical approximations have certain uniform integrability properties. 

2.1 General discrete-time stochastic processes 

This subsection introduces a general approach for studying integrability and stability properties of discrete-time 
stochastic processes. We assume a Lyapunov-type estimate on a subevent of the probability space for each time 
step. From this, we derive Lyapunov-type estimates for the process uniformly in the time variable in Subsec- 
tion 12.1.11 This approach is then applied to derive uniform moment bounds in finite time (see Subsection 12.1.31 
see also Proposition 12.11 for infinite time) for a large class of possibly infinite dimensional approximation pro- 
cesses. Note that the state space {E,£) appearing in Propositions 12.11 and 12 . 71 and in Corollaries l2.2lE751 and E^ 
is an arbitrary measurable space. In our examples in Subsections 12.21 and 12.31 below we restrict ourself, however, 
to explicit (see Subsection 12. 2p and implicit (see Subsection 12.31) approximation schemes for finite dimensional 
SDEs driven by standard Brownian motions. Our approach is influenced by ideas in [331 [51 134) : see the end of 
Subsection 12 . 1 .41 for more details on the these articles. 



2.1.1 Lyapunov-type estimates 

The following proposition proves Lyapunov-type estimates for discrete-time stochastic processes which do, in 
general, not hold on the whole probability space but only on a family of typically large subevents of the probability 
space. These subevents are defined in terms of an appropriate Lyapunov-type function V: E ^ [0, oo) on the 
measurable state space {E,£) and in terms of a suitable truncation function C,: [0, oo) — >■ (0,oo]. 

Proposition 2.1. Let p G M, let (fJ, J^, P) he a probability space, let {E,£) be a measurable space, let tn € K, 
n € No, be a non- decreasing sequence, let ( : [0, oo) — !• (0, oo] be a function, letV : i? — > [0, oo) be an £/B{[0, oo))- 
measurable function and let Y : Nq x fl ^>- E , Z : N x ^ ^ 'M. be stochastic processes with E[la„ \Zn\] < oo and 



ln„V{Yn) < e''(*"-*"-iV(r„_i) + lo„Z„ 
for allneN where fi„ := n)JZo{F(yfc) < C(ife+i - tk)} e J^ for all n G Nq. Then 



k=l 

' ePfo-«o) E[V{Yo)] + J2Li eP'*^-*'^ lEjlp, Zi] 



'[("«)1 < E 



fc=0 



C(^fc+i - tk) 



(10) 

(11) 

(12) 



for all n G No . 

Proof of Provosition \2A\ First, observe that assumption (fTO)) . the relation r2„ C rj„_i for all n G N and non- 
negativity of V show that 

ln„V^(>;0 < ln„_ie''(*"-*"-i)F(r„_i) + ln„Z^ (13) 

for all n G N. Estimate ([T^ is equivalent to the inequality 

ln^e-P'-V{Yr,) - lo„_,e-''*"-iy(r„_i) < la.e^"*"^™ 
for all n G N. Next note that (fT4| and the fact Oo = ri imply 



(14) 



la„e-^*"T/(y„) = l^.e-P'-ViYa) + ^ (lo,e-''*^T^(rfe) - In^^.e-P'^-^ViYu-i)) 



fe=i 



<e-^*°y(yo) + Ei"'=e'''*'^'' 



(15) 



fe=i 



for all n G No- This implies pT|) . For proving p^ . note that the relation J7„ C ^n-i for all n G N implies 

{^nT = (f^„-i \ f^„) w ( {^n-iT \ aO - (a.-i \ a.) w ( (a.-i)' ) (i6) 

for all n G N. Iterating equation (fTO)) and using again ilo = ^ shows 



{n^f = { [+J (i7fe\ fifc+i) J l+J ( (r!o)^ ) = l+J (r!fe\ nk+i) 

n— 1 n — 1 

- l+J (rifc n {y(n) > C(ifc+i - tk)}) - l+J {io,y(rfc) > C{tk+i - tfc)} 



(17) 



fc=0 fc=0 

for all n G Nq. Additivity of the probability measure P, Markov's inequality and inequality dill) therefore imply 

k^O L 



n— I n—i 

^[{n^)^']=Y,F[ln,V{Yk)>C{tk+i-tk)] <E 



fc=0 

^E 

fc=0 



C{tk+i ~ tk) 



ePJt.-to) E[V{Yo)] + -£^1 eP'^*"'''^ E[ln,Zi] 

({tk+l — tk) 



for all n G No- This is inequality P^ and the proof of Proposition 12 . 1 1 is thus completed. 



(18) 



D 



Let us illustrate Proposition 12.11 with the following simple implication. If the assumptions of Proposition 12. II 
are fulfilled, if p S (— oo,0), C = oo, E[F(yo)] < c» and if there exist real numbers h,c d (0,oo) such that 
supj.gpj E[Zfe] < ch and i„ = nh for all n E No, then we infer from inequality dill) that limsup„_^Q^ E[V^(y,i)] < 
c/i/(l - exp(p/i)) < i^e!''!'' < oo. 

In many situations, the random variables {Zn)neN appearing in Proposition 12.11 are centered or even appro- 
priate martingale differences. This is subject of the next two corollaries (Corollary 12.21 and Corollary 12. 3p of 
Proposition 12.11 

Corollary 2.2. Let p G M, let (fi, J-", P) be a probability space, let {E,£) be a measurable space, let tn G M, 
n e Nq, be a non- decreasing sequence, let C: [0, oo) — >■ (0, oo] be a function, letV : E ^> [0, oo) be an £/8([0, oo))- 
measurable function and let Y : Nq x fl -^ E be a stochastic process with E,[V{Yq)] < oo and 

^n-,.„{vin)<iit.^^-t,)} ■ ^[ViYn+i) I iYk)kem,...,n}] < 6"^*"+^-*"^ • ViY^) (19) 

P-a.s. for all n G No- Then the stochastic process In^e^''*" V^(y„), n G Nq; is a non-negative supermartingale 
and 

E[lo„n5^„)] < e''(*"-*«)E[T/(ro)], P[(^„)1 < ( E ^fa^^ _ t,) J'^[^(^°)] ' (20) 



E[y(r„)] < e''(*"-*")E[F(ro)] + \\nYn)\\L.inm ( E c{t,+^-t,) )^t^(^o)] 



n-l 



Pitk-ta) 



fe=0 



(1-1/p) 



(21) 



for all n G Nq, p G [l,oo] and all E/B{[Q,oo))- measurable functions V : E ^f [0, oo) with V{x) < V{x) for all 
xG E where f2„ := n^^o{F(Yfc) < C(ifc+i - tfc)} e J" /or aZ/ n G Nq. 



Proof of Corollary \2.2\ The relation r2„+i C fin for all n G No together with the non-negativity of V shows that 
assumption ([T^ is equivalent to the estimate 



E[la„+ie-''*"+^y(r„+i) I iYk)ke{OA,...,n}] < In^e^'^-ViY^) (22) 

P-a.s. for all n G Nq. Combining ([2^ . the assumption E[V^(yo)] < co and the relation 

<^n{iln,e~P'>'V{Yk))ke{o.j,...,n}) (^ MiYk)ke{o,i.....,n}) (23) 

for all n G No proves that the process lsi„e^''*"F(F„), n G No, is a non-negative supermartingale. This implies 
the first inequality in ([20)1 . In addition, this ensures that the stochastic process Z: N x 17 — >■ R given by Z„ = 
in„V(Yn) - E[ln„V{Y„) \ (n)fc6{o,i,....n-i}] P-a.s. for ah n G N satisfies E[lj2jZ„|] < oo and E[ln„Z„] = 
for all n G N. Moreover, the definition of Z : N x il — > M and assumption ([T5| ensure 

ln„V{Yn) - E[lo„ V(r„) I (Yk)ke{o.j,...,n-i}] + ^« 

= lo„E[y(r„) I (n)fce{o,i,...,n-i}] + lo„^« < e''(*"-*"-iV(y„_i) + lo„Z„ 

P-a.s. for all n G N. An application of Proposition 12.11 thus proves the second inequality in (|20l) . Next observe 
that Holder's inequality implies 

E[X] < ni^X] + {¥[inr]f-'^'''>\\X\\L.inm (25) 



for all ft G T, p E [l,oo] and all J^/;B([0, oo))-measurable mappings X : fl ^ [0,oo). Combining (PU]) and ([^ 
finally results in 



E[^(>;0] < E[ln„ViYn)] + ||nr„)|Up(0;R) (P[(f^„)1) 



(i-i/p) 



< e''(*"-*«)E[F(yo)] + r(y»)HLnO;K)[(^E ^(^,^^_^,) j^[^(^°)] 



(1-1/p) (26) 



for all n G No, p G [l,oo] and all £/i3([0, oo))-measurable functions V : E ^i' [0,oo) with V{x) < V{x) for all 
X G E. The proof of Corollarv l2.2l is thus completed. D 



Corollarv l2.2l in particular, proves estimates on the quantities supj,g{o,i,...,n} E[lOfc V^(yfe)] for n G N (see the 
first inequality in ^U\i). Here fin C fl, n G Nq, are typically large subevents of the probability space {fl,J^,F) 
and V: E ~> [0, oo) is an appropriate Lyapunov-type function (see Corollarv 12.21 for details). Under suitable 
additional assumptions, one can even obtain an estimate on the larger quantities E[supj.gr-o^i^...^„} InfcV'(^fc)] 
for n G N. This is subject of the next corollary. 



Corollary 2.3. Let p G M, p G [l,oo), let (51, J^, P) be a probability space, let {E,£) be a measurable space, 
let i„ G M, n G No, be a non- decreasing sequence, let (^: [0,cx3) — > (0,oo], v. N -^ [0, c») be functions, let 
V: E ^ [0, oo) be an £/B{[0,oo)) -measurable function and let Y : Nq x Vl —^ E, Z: N x 57 — >■ R 6e stochastic 
processes such that the process X]fc=i ~^^k^k, ri G N, is a martingale and such that 



la„T^(r„) < e''(*"-*"-iV(r„_i) + ln„Z„ 



lOn^nlU 



P(0;R) 



< V-n 



sup lo,e''(*"-*'=V(rfe) 

fee{0,l,...,n-l} 



LP(n;R) 



for allneN where 5]„ := r\'kZo{V{Yk) < C(ifc+i - tk)} e T for all n G Nq. Then 



sup In^e 
/£e{o,i, ■■-,"} 



-pt* 



y{yk) 



LviU-f 



<2||V'(l?))|liP(n;R)exp Xp 



El 

Lfc=l 



Vk\ 



Ph 



for all n G No . 



(27) 
(28) 



(29) 



Proof of Corollarv \2.3\ Inequality PTI) . the triangle inequality, the estimate (a + 6)^ < 2a^ + 26^ for all a, 6 G 
and the definition © of Xp G [0, oo), p G [1, oo), applied to the martingale X^ILi '^vtk^^''^'' Zk, n G No, yield 



sup ln,e- P'^ViYk, 

fee{0,l,...,n} 



< 2 e" 



-pto 



^(i"o) 



< 2 \\e-P^°V{Yo 



lLP(n;R) 



lLP(n:] 



<2||e-''*«ni^o)|L,(^^«) 



LP(n;I 



sup 

fce{l,2,...,n} 



Eif^^^ 



-ptl • 



1=1 



LP(0;B 



(30) 



2xp5;]||la,e-''*'=Zfc 



lLP(f2;I 



fc=l 
n-1 



2XpE 



i^fc+i| 



fe=0 



sup 
ie{o,i,---,fe} 



lo^e-"*' 



^(>^0 



LP(0;R) 



for all n G No where the last inequality follows from assumption (|28p . Consequently, Gronwall's lemma for 
discrete time yields 



sup la,e^''*^F(rfe; 

kl^{0,l,...,n} 



<2\\e 



-pta 



Lp{n-fi 



^(^o)|Lp(j,.R)exp 



/ ri-l 

2xpEI'^'^+iI 



k=0 



for all n G No- This finishes the proof of Corollary [ 

An application of Corollary 12 . 31 can be found in Lemma [2.281 below. 



D 



2.1.2 Semi moment bounds 

In the case of nonlinear SDEs, it has been shown in [35] that the Euler-Maruyama approximations often fail 
to satisfy moment bounds although the exact solution of the SDE does. Nonetheless, the Euler-Maruyama 
approximations often satisfy suitable moment bounds restricted to events whose probabilities converge to one 
sufficiently fast; see Corollary 4.4 and Lemma 4.5 in _33 . This is one motivation for the next definition. 

Definition 2.4 (Semi boundedness) . Let a G (0, oo], let L C M. be a subset o/K, let {E,£) be a measurable space, 
let {Q,J-,¥) be a probability space and let V : E ^- [0, oo) be an £/B{[0,oo)) -measurable mapping. A sequence 
Y^ : Lxfl~?'E,N£N, of stochastic processes is then said to be a-semi V -bounded (with respect to V) if there 
exists a sequence ft^ G o'n{{Y/^)tei) — o'niY^) C J-, N E N, of events such that 



limsup supE[ln„l^(i;^)] +N"-¥[{nNy] < oo. 



(31) 



Moreover, a sequence Y^ : /x 57 — > _E, iV G N, of stochastic processes is said to be 0-semi V -bounded (with respect 
to Pj if there exists a sequence Vti,i G aii{Y'^), N G N, of events such that lim sup^_j.oQ sup^gj ]E[ls7„y(y("'^)] < 
oo and \\v[\i^^ao^\{S^N)\ — 0. 



We now present some remarks concerning Definition 12.41 First, note that the concept of semi boundedness 
in the sense of Definition 12.41 is a property of the probabihty measures associated to the stochastic processes. 
More precisely, in the setting of Definition 12.41 a sequence of stochastic processes Y^ : I x ft ^ E, N G N, is 
a-semi V^-bounded if and only if there exists a sequence An G S^'' , -/V G N, of sets such that 

limsupi sup / lAr,{x)V{Trt{x))PY«{dx)+N°' -WYNliANy]} <oo (32) 

Af->oo yte[0,T]JE^' J 

where Py« : £'>^' -> [0, 1], A^ G N, with Pyiv [A] = P[y^ e A] for all A e S"^^ and aU A^ e N are the probability 
measures associated to Y^ : I x il -^ E, N E N, and where vrt : E^' -^ E, t (1 I, with Trt{x) — x{t) for all 
X G E^^ and alH G / are projections from E^^ to E. Semi boundedness in the sense of Definition 12.41 is thus 
a property of the sequence Pyw : £^^ — >• [0, 1], A^ G N, of probability measures on the measurable path space 
{E^-' , £®^) . Moreover, observe that in the case where the index set / appearing in Definition 12 .41 consists of only 
one real number to G K, i.e., / = {to}, the sequence Y^ : I x H. ^ E, N G N, oi stochastic processes reduces 
to a sequence Z^: H. —>■ E, N E N, oi random variables with Zn{uj) = Yj^ (uj) for all w G il and all TV G N. 
Furthermore, observe that in the case a = cx) in Definition 12. 4[ condition (PT|) is equivalent to the condition 
limsupjv_j.oo sup(g/E[T^(y(^)] < oo. Next we would like to add some further comments to Definition 12.41 For 
this we first note the following well-known characterization of convergence in probability (see, e.g.. Exercise 6.2.1 
(i) in Klenke [40] and Remark 9 in [38]). 

Lemma 2.5 (Convergence in probability). Let (r2,J^, P) be a probability space, let (EjcIe) be a separable metric 
space and let X : D, ^ E and Y/v : fl ^ E, N Cz N, be J- /B{E) -measurable mappings. Then the following three 
assertions 

(i) the sequence Y^ , iV G N, converges to X in probability, i.e., limAr_j.oo P[c^£;(Ar, Yat) > e] = for all 
e G (0,oo), 

(a) it holds that l[mN^ooP[dEiX,YN) < l] = 1 and limTv^oo E[l{rfB(x,yjv)<i}'^B(^' ^A^)] == 0' 
(Hi) there exists a sequence fi^r G J^, TV G N, with liinjv_i.oo ]E[ls7„d£;(X, Y/v)] — and limjv-j-o P [^w] = 1 
are equivalent and each of these assertions implies 

\imN^ooV[dE{X,YN) > l] =0 and sup^grjE[l{rf^(x,Y„)<i} Me(2;, Ytv)!''] < oo (33) 

for all X E E with E[ \dE{x, X)\^ ] < cxi and all p G (0, cxd). 

Proof of LemmalKR Note that if limN^oof[dE{X, Yn) > e] = for all s G (0, oo), then liniAr^ooP[rf£:(-'f , Yn) < 
l] = 1 and limjv_j.ooE[l{(i^(x.y„)<i}rf_E(-'^, ^Af)] = due to Lebesgue's theorem of dominated convergence. This 
shows that (i) implies (ii). Next observe that (ii), clearly, implies (iii). Moreover, if there exists a sequence Jljv G 
J", iVGN, with limAr^ooE[ln„d£;(X,F^)] = and limAr^oP[^Jv] = 1, then limjv^ooP[ln«c«i=;(^,>^Ar) > e] = 
for all £ e (0, oo) and therefore 



lim ¥\dEiX,YN) >e]< lim Pff^^v n {^^(A:, Fjv) > e}] + lim F\{nNY 

Af->oo Af-i-oo JV-i-oo 

= lim ¥[ln^dE{X,YN)>e] =0 



(34) 



for all e G (0,00 ). This proves that (iii) implies (i). Finally, observe that 

||l{dE(x,y„)<i}t^s(2;,>V)||iP(s^.R) 

< ||l{dE(x,Y-„)<i}'^B(2^7-'^)|Lp(n;K) + ||l{<iB(x,y«)<i}C?£(-'^'^Jv)||iP(o.R) (35) 

for all a; G -B, A^ G N and all p e (0, oo). The proof of Lemma I^TSl is thus completed. D 

Let us now study under which conditions Euler-Maruyama approximations are a-semi ^-bounded with 
a G [0,oo] and y : M'' — ?► [0,oo) appropriate and d G N. First, note that convergence in probability of the Euler- 
Maruyama approximations has been established in the literature for a large class of possibly highly nonlinear 
SDEs (see, e.g., Krylov [JS], Gyongy & Krylov [35], Gyongy [3T] and Jentzen, Kloeden & Neuenkirch iM ). For 
these SDEs, one can thus apply Lemma [2.51 to obtain the existence of a sequence of events whose probabilities 
converge to one and on which moments of the Euler approximations are bounded due to p3p in Lemma 12.51 
This is, however, not sufficient to establish a-semi ||-|| ''-boundedness of the Euler-Maruyama approximations 
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with a G [0, oo] and p S (0, cxi) since the events in p3p do, in general, not satisfy the required measurability con- 
dition in Definition 12.41 Moreover, we are mainly interested in a-semi |j-|j^-boundedness of the Euler-Mayuama 
approximations with a,p S (0, oo). Lenima l2.5l onlv shows that the complement of the events on the right side of 
([55)1 converge to zero but gives no information on the rate of convergence of the probabilities of these events. So, 
in general, a-semi ||-||^-boundedness with a,p S (0,oo) cannot be inferred from convergence in probability. Here 
we employ the theory of Subsection 12.1.11 to obtain semi boundedness for the Euler-Maruyama approximations 
(see Theorem 12.131 and Corollarv l2.9l below'). In particular, Corollarv 12 . 2 1 immediatelv implies the next corollary. 

Corollary 2.6 (Semi moment bounds). Let (r2,J^, P) he a probability space, let {E,£) be a measurable space, 
let p,T Cz (0, oo), a G (1, oo], q G (1, oo], A^o €^,letV:E-~> [0, oo) be an £/B{[0, od)) -measurable mapping and 
let Y^ : {0, 1, . . . , N} x il — > E', A^ G N, be a sequence of stochastic processes satisfying 

l{V(y«)<| « p} ■^ViY^'i,) I (n^)fee{0,l,...,n}] < e"^ ' ^(^rf ) (36) 

P-a.s. for aline {0,1,..., N- I}, N G {Nq, Nq + 1, . . .} and \imaupj^^^E[V{YQ^)] < oo. Then the sequence 
Y^ : {0, 1, . . . , N} X n —?' E, N e N, of stochastic processes is {a — l)-semi V -bounded with respect to P. 

Proof of Corollarv \2.6\ Corollarv 12.21 above with the truncation function C,: [0, oo) — > (0, oo] given by C(0) = 
and C(i) — t^" for all t G (0, oo) and with the sequence i„ G M, n G No, given by t„ ~ nT/N for all n G Nq 
implies 

E[lo„y(r,r)] < e^ • E[y(ro^)] < e"^ • E[V{Y,'')] and (37) 



pii^Ny 



< 



\k=0 



Y: e^ -^ViY,^)] < i- . iV . e"^ • ^[^(^0^)] (38) 



for aU 71 G {0, 1, ... , N} and aU N G {Na,No + 1, . • . } where nN := n^=o{V{Y,^^) < [T/Nf] G J" for all 
A^ G N. Combining ([57)) . ([55]) and the assumption limsup^^^^ E[y(yg^)] < oo implies 

limsup ( sup E[l,-i^V{Y^)\ +iV("-i) •P[(r!w)^] ) < oo. (39) 

Af^oo \ ne{0,l,...,W} / 



The proof of Corollarv 12.61 is thus completed. D 

2.1.3 Moment bounds 



Corollarv l2.61 in particular, establishes moment bounds restricted to the complements of rare events for sequences 
of stochastic processes. In some situations, the sequence of stochastic processes fulfills an additional growth 
bound assumption (see (j4ip and (|46)) below for details) which can be used to prove moment bounds on the 
full probability space. This is subject of the next proposition. The main idea of this proposition is a certain 
bootstrap argument which exploits inequality (pi]) in Corollarv l2.2l (see also estimate (US])). 

Proposition 2.7 (Moment bounds). Let (17, J^;P) be a probability space, let {E,£) be a measurable space, let 
T G (0,oo), p G [0,00), a G (l,oo], A^o G N, let V,V : E -^ [0,oo) be £ /B {[0, 00)) -measurable mappings with 
V{x) < V{x) for all X e E and let Y^ : {0, 1, . . . , A^} x fi — > ii', A^ G N, be a sequence of stochastic processes 
satisfying 

l{y(y«)<| « PI • E[l/(y„^i) I (l/)fce{0,l,...,n}] < 6^ ' V {Y^ ) (40) 

P-a.s. for allne{Q,l,...,N -1} and all N G {A^o, A^o + 1, ■ • ■ }■ Then 

limsup sup E[y(y„^)] (41) 

W-i-oo ne{Q.l,....N} 

< e''^ fl + lim sup E [1/(1/)]) fl + T^^^-i/p) limsup Ui-")^-!/^) sup \\V{Y^)\\L,(n:W,) 

\ N^oo J \ N^oo L 0<n<Af 

for all p G [1, oo]. 

Proof of Proposition \2. 7[ With the truncation function C,: [0, oo) — ?> (0,oo] given by C(0) = and C,{t) = i~" 
for all t G (0, oo) and with the sequence i„ G K, n G No, given by i„ = nT/N for all n G No we get from 
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inequality (PT|) in Corollary 12.21 that 

^ViY,^)] < e'^ElViY,^)] + ||nr„^)|U.(a;i 



«-i ^^^\ i(i-i/p) 



fc=0 



^e^ E[n>^o^)] 



(i-i/p) (42) 



< e"^ (1 + E[ViY,^)]) (l + I T p(i-Vp) ^(1-1/,) ||y(yJV)||^^^^^^^) 
for all n e {0, 1, . . . , A^} and all N eN. This completes the proof of Proposition 12. 71 D 

2.1.4 One-step approximation schemes 

Let T e (0, oo), let {il, T , {Tt)t^\Qj:\,^^ be a filtered probability space and let {E, £) be a measurable space. We 
are interested in general one-step approximation schemes for SDEs with state space E. The driving noise of the 
SDE could be, e.g., a standard Brownian motion, a fractional Brownian motion or a Levy process. In such a 
general situation, a sequence of approximation processes of the solution process with uniform time discretization 
is often given by a measurable mapping \I>: i? x [0,r]^ x Jl — ^ E as follows. Let ^: fi ^ _E be a measurable 
mapping. Define stochastic processes Y^ : {0,1,..., A^} x fi ^ i?, A^ e N, through Yq = ^ and 

>;f+i = *(>;f,^,^^^) (43) 

for all n e {0, 1, . . . , A^— 1} and all A^ G N. For example, in the setting of the introduction, the Euler-Maruyama 
scheme ([3]) is given by the one-step function *(x, s,t) — fi{x){t~ s) + (j{x)(Wt — Wg) for all a; G M'', s,t £ [0,T]. 
The general approach of Subsections 12.1.114^71751 can be used to study moment bounds for stochastic processes 
defined as in (j43l) . 

In Subsections 12.21 and 12 . 31 below, we focus on finite-dimensional SDEs driven by standard Brownian motions. 
Due to the Markov property of the Brownian motion, an approximation scheme can in this case be specified 
by a function of the current position, of the time increment and of the increment of a Brownian motion. More 
precisely, let d,m €N, € (0, T] and let $ : M'* x [0, 9] x R™ -S> K"^ be a Borel measurable function. Using this 
function and a uniform time discretization, we define a family of stochastic processes Y^ : {0,1,..., A^} xft -^ R**, 
A^ e N n [f , cx)), through F(f = ^ and 

>;^i = *& {Yn , w ' Wi^ii^ - W^ ) (44) 

for all n G {0, 1, . . . , A^ — 1} and all A^ G N n [j-, oo). In this notation, the Euler-Maruyama scheme ([3]) is given 
by ^{x,t,y) — p,{x)t + a{x)y for all {x,t,y) G M'' x [0,6'] x R™. Moreover, in this setting, condition P^ in 
Corollarv 12.21 on the approximation processes Y^ : {0, 1, . . . , A^} x J7 — >• R'', A^ G N, follows from the following 
condition (see (gS])) on the one-step function $: R"^ x [0,0] x R'" -^ R^. 

Definition 2.8 (Semi stability with respect to Brownian motion). Let 9 G (0, oo), a G (0, oo], d, ?ti G N and let 
y : R"^ — > [0, oo) 6e a Borel measurable function. A Borel measurable function $ : R'' x [i),9] x R™ — > R"^ is then 
said to be a-semi V -stable with respect to Brownian motion if there exists a real number p G R such that 

E[V{^{x, t, Wt))] < e''* • V{x) (45) 

for all {x,t) G {(y, s) G R*^ x (0, 6*] : a = oo or V{y) < s^"} where W : [0,9] x ^ -^ R™ is an arbitrary standard 
Brownian motion on a probability space (f2, T, P). In addition, a Borel measurable function $ : R'' x [0, 6] x R™ — > 
R'* is simply said to be V-stable with respect to Brownian motion if it is oo-semi V -stable with respect to Brownian 
motion. 

Let us add some remarks to Definition 12.81 Inequalities of the form (I45p with a = oo have been frequently 
used in the literature; see, e.g.. Assumption 2.2 in Mattingly, Stuart & Higham [51] and Section 3.1 in Schurz [55] . 
In particular, inequality (16) in Schurz [55] (see also Schurz [55]) defines a numerical approximation of a similar 
form as (|43]) to be (weakly) F-stable if inequality (|45]) holds with a = oo where V : W'' -^ [0,oo) is Borel 
measurable. 

In the case a G (0, oo), inequality (jiS]) in Definition!^ is restricted to the subset {y G R'^ : V{y) < t^"} C R"* 
which increases to the fuh state space R'* as the time step-size t G (0, T] decreases to zero. Roughly speaking, 
the parameter a G (0,oo) describes the speed how fast the subsets {y G R'^: V{y) < i~"}, t G (0,T], increase to 
the full state space R'^. 

Next note that if G (0,oo), /3 G (0,oo], d, to G N, if V^: R'' — >■ [l,oo) is a Borel measurable function and if 
$: R'^ X [0,0] X R™ — >• R'^ is ^-semi ^/-stable with respect to Brownian motion, then $ is also a-semi F-stable 
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with respect to Brownian motion for every a S (0, /3]. However, the converse is not true in generaL In particular, 
Theorem 12.131 and Corollary 12 . 1 71 show in many situations that the Euler-Maruyama scheme is a-semi V^-stable 
with respect to Brownian motion for some a € (0,cx)) but not V-stable with respect to Brownian motion if the 
coefficients of the underlying SDE grow superlinearly. 

It follows immediately from the Markov property of the Brownian motion that the approximation pro- 
cesses (j44l) associated to an a-semi F-stable one-step function satisfy condition (p6| where a e (0, oo) and 
y : R'' — >■ [0, oo) are appropriate and where d e N. The next corollary collects consequences of this observation. 

Corollary 2.9 (Semi moment bounds and moment bounds based on semi stability with respect to Brownian 
motion). Let rf, m G N, T G (0,oo), 9 g (0,r], a g (l,oo], p g [l,oo], let V,V : R"^ -^ [0, oo) be Borel 
measurable mappings with V{x) < V{x) for all x g M.'^, let (O, J^, {J^t)te[o,T]i^) ^6 a filtered probability space, let 
W: [0,T]xn~^ M™ be a standard {Tt)telo,T]- Brownian motion, let<^: R'^ x [0, 9] x M" -^ R'^ be a-semi V -stable 
with respect to Brownian motion and let Y^ : {0, 1, . . . , N} x J7 — > R**, N G N, be a sequence of {J^t)te[o.T] -adapted 
stochastic processes satisfying liTasup]^^^K[V{Yf^)] < oo and l^i = $(F^, j^,W(^n+i)T/N — M^tiT/a') for all 
n G {0, 1, . . . , A^ - 1} and all N £ N n [T/9, oo). Then the stochastic processes Y^ , N e N, are {a - l)-semi 
V -bounded with respect to P. Moreover, if 

hmsup (7v(i-)(i-i/rt . sup„g{o,i,...,iV} \\ViY^)\\Lnn-M)) < oo (46) 

A''— ^oo 

in addition to the above assumptions, then the stochastic processes Y^ , N £ N, are also co-semi V -bounded with 
respect to P, i.e., limsupjv_).oo sup„g{o,i,---,JV} ^[^(^ri^)] < o°- 

Corollary 12.91 is an immediate consequence of Corollary 12.61 and Proposition 12.71 Below in Subsections 12.21 
and 12.31 we will study both explicit and implicit one-step numerical approximation processes of the form (j44p . 
For these approximation processes we will then give sufficient conditions which ensure that the function $ : R"^ x 
[0, 9] x R™ — J> R'^ appearing in pi)) is a-semi ^-stable with respect to Brownian motion so that PO]) is fulfilled 
and the abstract results developed in Subsections l2.1.1Ii^T01 can thus be applied. 

To the best of our knowledge, the idea to restrict the Euler-Maruyama approximations to large subevents of 
the probability space which increase to the full probability space as the time discretization step-size decreases 
to appeared first in [33] (see Section 4 in [33] and also [34] for details). The idea to restrict a Lyapunov- 
type condition of the form (|45|) on the one-step function of (Metropolis-adjusted) Euler-Maruyama schemes to 
subsets of the state space which increase to the full state space as the time discretization step-size decreases to 
appeared first in Bou-Rabee & Hairer [S] in the setting of the Langevin equation (see Lemma 3.5 and Section 5 
in [F for details). Finally, a bootstrap argument similar as in ([^T|) and Proposition 12 . 71 appeared first in [33] for 
a class of drift-tamed Euler-Maruyama approximations (see the proof of Lemma 3.9 in [33] for details). 

2.2 Explicit approximation schemes 

This subsection investigates stability properties and moment bounds of explicit approximation schemes. We 
begin with the Euler-Maruyama scheme in Subsection 12.2.1] and then analyze further approximation methods in 
Subsections I^X^ and 1^X31 below. 

2.2.1 Semi stability for the Euler-Maruyama scheme 

The main result of this subsection. Theorem 12.131 below, gives sufficient conditions for the Euler-Maruyama 
scheme to be a-semi ^/-stable with respect to Brownian motion with a g (0,oo), V : R'' -^ [0,oo) appropriate 
and d g N. For proving this result, we first present three auxiliary results (Lemmas I2.10I[2TT2| . The first one 
(Lemma I2.10p establishes a-semi V^-stability with respect to Brownian motion of the Euler-Maruyama scheme 
with a g (0,oo), V : R'^ — 7> [0,oo) appropriate and d g N under a general abstract condition (see inequality (|Tf)) 
below for details). 

Lemma 2.10 (An abstract condition for semi F-stabilitiy with respect to Brownian motion of the Euler- 
Maruyama scheme) . Leta,T (E (0,oo), d, m g N, let n: R'^ -^ R'^, a = (CT,j)jg{i,2,...,d}je{i,2,...,m} : R'^ ^ R'^''"' 
be Borel measurable functions, p,p g R, let (r2,J^, P) be a probability space, let W: [0,T] x ft — > M™ be 
a standard Brownian motion and let V : R'^ — > [0, oo) be a twice continuously differentiable function with 
iG^.aV){x)<p-V{x) and 

E[{gf,^„V){x + ^iix)-t + a{x)Wt,x) - {g^^„V)ix)^ <p-Vix) (47) 

for all {x,t) g {{y,s) eR'^ x (0,T]: V{y) < s""}. Then 

E[V{x -f ii{x) ■ t + a{x)Wt)] < e^P+P^' ■ V{x) (48) 
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for all {x, t) G {(y, s) G R** x (0, T] : V{y) < s "}. In particular, the Euler-Maruyama scheme is a-semi V -stable 
with respect to Brownian motion. 

Proof of Lemma \2.1(A Ito's formula shows 



E 



V{x + fi{x) ■t + a{x)Wt) 
V{x)- 



ds 



{G,,,aV) {x + ^i{x)-t + a{x)Wt , X 

= Vix) + {gV)ix)-t+ f E\{g^,„V){x + iiix) ■t + aix)Wt,x) - {g^^,v)ix) 
Jo '- 

<Vix)+p-t-Vix)+ [ E\{g^,,V){x + fi{x)-t + aix)Wt,x)-{g^,aV)ix) 
Jo '- 

for all X G M'^ and all i G [0,r] and inequality (H7|) therefore implies 

E V{x + n{x) ■ t + a{x)Wt) < V{x) +p-t- V{x) +p-t- V{x) 

= V{x) {l + p-t + p-t)< e^P+P^' ■ V{x) 

for all (a;, t) G {{y, s) eR'' x (0, T] : V{y) < s""}. The proof of Lemma EHH is thus completed. 



(49) 



ds 
ds 



(50) 
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In many situations, it is tedious to verify (j47p . We therefore give more concrete conditions for a-semi Te- 
stability with respect to Brownian motion of the Euler-Maruyama scheme with a G (0,oo), V: R'' —>■ [0, oo) 
appropriate and d G N in Theorem 12.131 below. For establishing this theorem, the next simple lemma is used. 

Lemma 2.11. Let T G (0,oo),c G [0,oo),p G [l,oo) he real numbers and let y: [0,T] -^ M. be an absolute 
continuous function with y'{t) < c|y(i)| for \\[QT]-a'lmost all t G [0,r]. Then 



y{t)< |y(o)l 



|1/P I ct 



< 2^'P- 



|y(o)l 



y{t) < |2/(0)| + ct [ |j/(0)|^/^ + cty'' 'U \y{0)\ + 2(^-1) [ct |y(0)|(^-^/^) + Ict^ 
for all t G [0, T] . 
Proof of Lemma \2. 11\ The assumption y'{t) < c\y{t)\ '^' for A|[o j-j-almost all t G [0,T] implies 

'y'{s){yis)f/P-'Us<c{t-to) 



(51) 
(52) 



(53) 

for all to G {T{t),t] and aU t G [0,T] with y{t) > where the function r: [0,T] -^ [0,T] is defined through 

r(t) ~ maxfjO} n {s G [0, t] : y{s) = O}] (54) 

for all i G [0,r]. Estimate ^^ then gives 



{y{t)f'<{y{to)y^' + '-^^ 

for all io e iT{t),t] and aU t G [0,r] with y{t) > 0. This yields 

p - I I p 



{y{t)r'<\yir{t))\ 



(55) 



(56) 



for all t G [0,T] with y{t) > 0. Inequality ([5^ implies ([5T|) . In addition, note that combining ([5T|) and the 
inequality 



(a: -I- ?/)''=: x'' -I- / r (a; + sy) y ds < x^ + r y {x + y) 



for all x,y G [0, oo), r G [1, oo) shows ([5^ . The proof of Lemma [2. Ill is thus completed. 
An immediate consequence of Lemma l2.11l are the following estimates. 



(57) 
D 
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Lemma 2.12. Let c,p S [l,oo) be real numbers and let V e C^{R'^,M) with ||(VF)(a;)|| < c\V{x)\^'^ ^^''^ 
X^d -almost all x G W^ . Then 

V{x + y) < |y(x)| +cP2(P-i) (||y|| \Vix)\^'-'/''^ + Hz/f 



for 

(58) 
(59) 



for all x,y € R"^. 

Proof of Lemma \KM The assumption ||(Vy)(x)|| < c |F(x)|^^"^/^^ for ARd-almost all x ^W^ implies 

^ {V{x + ty)) = V'{x + ty)y< \\VV{x + ty)\\ \\y\\ < c |T/(x + ty)\^^-^^P^ \\y\\ 

for Ar- almost alH G M and all a;, y G R"*. Inequality ([5^ in Lemma [2.111 hence gives ([5^ and 

V{x + y) < 2(P-i)(|F(x)| + |^[) < c^2(^-i)(|l/(x)| + \\yf) 

for all x, y G M'^. The proof of Lemma [2. 121 is thus completed. 



(60) 



(61) 
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We are now ready to present the promised theorem which shows a-semi V-stability with respect to Brownian 
motion of the Euler-Maruyama scheme with a G (0, oo), T^: R"^ — >■ [0, oo) appropriate and d G N. It makes use 
of Lemma 12.101 and Lemma 12.121 above. 

Theorem 2.13 (Semi ^-stability with respect to Brownian motion for the Euler-Maruyama scheme). Let 
T G (0,oo), d, 771 G N, p G [3,oo), c, 70,71 G [0, 00) be real numbers with 70+71 > 0, let /i: R'' — > R'', 



hdxm 



be Borel measurable functions and let V G Cp(R , [1, 00)) with 



{g,..aV){x)<c-V{x), \\fi{x)\\<c\V{x)\i^\, ||a(x)||i(R™,R.)<c|F(a;)|i^^J (62) 

for all X G R''. Then the Euler-Maruyama scheme R'^ x [0,T] x R™ 3 {x,t,y) 1-^ x + ^.{x)t + a{x)y G R"* is 
plili + 2(70 V 7i))-semi V-stable with respect to Brownian motion. 



Proof of TheoremlMM Throughout this proof, let {n,T,P) be a probability space and let W : [0,r] x 17 ^ R™ 
be a standard Brownian motion. We will prove semi T^-stability with respect to Brownian motion for the 
Euler-Maruyama scheme by applying Lemma 12.101 Our aim is thus to verify (|47l) . For this note that 



E 



\{g,.,aV){x + ll{x) ■ t + <j{x)Wt,x) ~ {g^,,,v){x)\ 



< E 



\V' {x + n{x) ■ t + a{x)Wt) -~V'{x)\ 



LiW^,] 



lM(a;)|| 



•E 



/ m 

\V"{x + fi{x) ■ t + a{x)Wt) - ^^"(a;)|L(.,(R.,R)] E 11^^(^)11 



(63) 



\k=l 



<m{c+lf j^E[||F«(a; + M(2;)-i + cr(a;)W"t) ~V^'\x)\\^^^ 



\Vix)\ 



[Z,^]- 



for all X G M'' and all t G [0,r]. Next observe that the assumption V G Cp{R'^, [l,oo)) and estimate (|55|) imply 
the existence of a real number c G [c + 1, 00) such that 



y«(2;)-yW(x) 



< 



V'-'+^\x + r{y ^ x)) 



L(> + i)( 



\y — x\\ dr 



L(*)(R<',R) Jo 

<ci I \V[x + r{y~ x))\' " 'drj|jy-x|| 

< |c|(^+^) 2^ {\V{x)\'^^i^ + \\y ~ xt—'^) \\y - x\\ 
= \c\^'+'hp[\V{x)\''^^\\y-x\\+\\y~xt-''>) 



(64) 
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for all x,y eR'^ and alH G {1, 2}. Putting ((Ml) into ^^ then results in 



E [I {G,.,aV) {X + tl{x) ■ t + <j{x)Wt,x) - {Qy.,ay){x) 

<to|c|'^+^)2'' [^E[||/x(a;)-i + CT(x)VKi||] \V{x)\ 



[ ""-^^'7 



(65) 



+ TO C 



(p+3) 



2W^E[||Ai(:E)-t + c7(x)Wt||(^-^^] IV^C 
for all a; G M'' and all t G [0, T\. In the next step, we put the estimate 

E 



^(i-i)+' 



Pi{x) ■ t + <j{x)Wtf\ < (m + ly ||^(:r)|r ■t'' + J2 W'^^ixW ' ^ 

V fc=i 

< c"" (to + r + l)(2'-+i) l' |T/(x)|[^^] f + |y(x)|[^^^^^] ti) 



W, 



(fc)| 



(66) 



c'' {m + r + 1) 



(2r+l) 






tl 



for all a; G M'', t G [0, T], r G [0, oo) into ([53]) to obtain 



E 



\{g^,aV){x + fi{x)-t + <j{x)Wt,x)-{g^^„V){x)\ 

't-(.-i)+"'(,-i)/J 



<(TO+p)^ncr'^+"2P Yl 1^(^)1 



"St 



viJ=i 



(67) 






[ "'"-^'^""/'""-^''V l] ^iE_2l^ 



for all X G M'' and all t G [0, T]. Lemma r2 . 1 01 therefore shows that the Euler-Maruyama scheme is a-semi V-stable 
with respect to Brownian motion, i.e., there exists a real number p G M such that E[y(a; + fi{x) ■ t + a{x)Wt)] < 
eP* ■ V{x) for all {x, t) G {{y, s) G M'* x (0, T] : V{y) < s""} where 



a := mm | mm^ 

ije{i:2} 



p-ip-i) 



J • (7(.-i) + ^) 'j- • (7(.-i) + ^^^^^) 



Finally, note that 



a = mm min< 
i,ie{i,2} ' 



J- 7(»-i) + 



P 



7(3-1) 



fi'-i-) I 7(3-1) 



J ■ 



-^) ^ J 



(69) 



^,je{i,2} \^ (j . 7(,_i) + 7(^_i)) y maxje{i,2} (j ■ niax(7o,7i) + 7o--i)) 
P 



71 +2max(7o,7i)' 
This completes the proof of Theorem 12.131 
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If y: M'' ^- [l,oo) is a Lyapunov-type function in the sense of Theorem 12.131 with d G N, then the function 
M'' 9 a; 1-^ |y(a:;)|' G [l,oo) with q G (0, oo) appropriate is in many situations a Lyapunov-type function too. 
This is subject of Corollary 12.151 below. The next result is a simple lemma which will be used in the proof of 
Corollary [2T51 

Lemma 2.14 (Dcriyatiyes of potencies). Let q G [l.oo), p G [2,c»), d G N and let V G Cp{M.'^, [l,oo)). Then 
the function V-.M!^ ^ [l,oo) given by V{x) = {V{x)f for all x &W^ satisfies V G Cl^^W^, [l,oo)). 
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Proof of Lemma \2.14\ Note that 

V'{x){v^) = q (l/(a;))(«-'^ V'{x){v^), (70) 

V"{x){v^,V2) ^q{q-l) {V{x))^''-^^ V'{x){v^) V {x){v2) + q {Vix))^"-^^ V"{x){vuV2) (71) 

for all X, wi, W2 G R'' and 

f(3)(^)(„^,„2, „3) = 9 (5 - 1) (g - 2) {V{x)f^-^^ V'{x){v^) V'{x){v2) V'{x){v:i) 
+ q{q-l) {Vix))^"-^^ V"{x){v,,V2) V'{x){v:i) 

+ q{q-l) {Vix))'^"-^^ V"{x){vi,v^) V'{x){v2) (72) 

+ q{q-l) {V{x)i''-^'> V'{x){v^) V"{x){v2,v^) 
+ q{V{x))'^^-^^V'^^\x){v,,V2,V:i) 

for all wi, W2, v-i e M'* and all a; e {y e R"^ : y" differentiable in y}. The assumptfon V e C^iW^, [1, 00)) therefore 
shows that V G Cpg(M'^, [l,oo)). The proof of Lemma [2.14l is thus completed. D 

Corollary 2.15 (Potencies of the Lyapunov-type function). Let T g (0,oo), d, jti G N, p G [3, 00), q G [l,oo), 
c, 70; 7i G [0, 00) he real numbers with 70 + 71 > 0, let jjl : W^ — > W^ , a : M'' — > M''^™ he Borel measurable functions 
and let V eCl{W^, [I, 00)) with \\ii{x)\\ <c\V{x)\\-^\ lk(a;)llL(R™,B<i) <c|y(x)|[^] and 



{q-l)\\ V'{x)cj{x)\ \lsi:^„ 
2 • F(a;) 



(ep,.^)(:^) + ^^ ^__Z_A_^^i£l£:i£l < c . y(x) (73) 



/or aH X G M''. T/ien i/ie Euler-Maruyama scheme R"^ x [0,T] x M™ 9 {x,t,y) 1-^ x + ^{x)t + a{x)y G M"* is 
pq/(pii + 2(70 V 71)) -semi |V^|''-staWe wit/i respect to Brownian motion. 

Proof of Corollary \MJ5\ First, define the function f: R'^ -> [l,oo) by V{x) = (T^(x))'' for all x G M''. Then note 
that this definition ensures 

\\fiix)\\<c\V{x)\i^^\ and Mx)]]^^^^^^.) < c\Vix)\i^^\ (74) 

for all a; G R'^. In addition, observe that 

t T\ ( ('^ ~ 1) ll^'(2;)cr(a;)|| rro^p™ luA 

{Q,Ay'-)){^)-r{V{x)t-^^ Uq,,.V){x)^^- \.y\^^ ^ M (75) 

for all X G R'* and all r G (0,00). Therefore, we obtain 

(e,,.F)(x) . , (V^(.))'-^' ((a.,.T/)(x) + ^'''^''"'2%"^''"^^""''"^ ) (76) 

for all X G R'^ and ([75)1 hence gives 

(e^,aV^)(x)<g-c-1/(x) (77) 

for all X G R*^. Combining (ITU) . ([77]) . Lemma [2.141 and Theorem 12.131 then shows that the Euler-Maruyama 
scheme is — , ^f'' ^ — y-semi F-stable with respect to Brownian motion. The proof of Corollary 12.151 is thus 
completed. D 

Note that in ([75)1 the norm in the Hilbert space of Hilbert- Schmidt operators from R™ to R is used where m G 
N. The definition of that norm and more details on Hilbert-Schmidt operators can, e.g., be found in Appendix B 
in Prevot & Rockner [5D]. In the next step, Theorem 12.131 is illustrated by a simple example. More precisely, 
the next corollary considers the special Lyapunov-type function T^: R'* — >■ [l,oo) given by V{x) = 1-1- ||x||^ for 
all X G R'^ with p G [3, 00) and d G N. 

Corollary 2.16 (A special polynomial like Lyapunov-type function). Let T G (0,oo), 0,70,71 G [0,oo), p G 
[3,00), d,m £N be real numbers with 70 + 71 > 0, let ^: R'' — > R'', a: R'' — > M''^'" be Borel measurable functions 
with 

(x, Ai(x)) + (£^||a(x)||^5(„„ R.) < c (1 + ||xf ) , (78) 

||a*(x)||<c(1 + ||x||[^«+i1) and ||a(x)|U(M„.,K.) < c (l + ||x||[^] ) (79) 

for all X G R''. Then the Euler-Maruyama scheme R'^ x [0,T] x R™ 9 {x,t,y) ^ x + ^j,{x)t + a{x)y G R'' is 
p/{li + 2(70 V 7i))-semi (1 + \\x\Y')^^^d- stable with respect to Brownian motion. 
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Proof of Corollary\ME Let V : W^ ^ [l,oo) be given by V{x) = 1 + \\xf for all x g W^. Then note that 
inequality ([75]) implies 

(^M,.l^)(^) < pII^II^""'^ (^,A^(^)) + £^ ll^ll^""'^ lk(^)ll?.s(E"^.«^) < '^■P-C-V{x) (80) 

for all a: e M''. Next observe that 

Mx)\\ < c (1 + ||a;|| [^"+'1 ) < 2c |F(a;)| t^] , (81) 

lk(2^)llL(R™,R<')<c(l + ||a;||[^])<2c|y(x)|[^] (82) 

for all X e R''. Combining (gOll-dHSl), the fact V G C^(IR'^, [1, oo)) and Theorem \n^ hence shows that the 
Euler-Maruyama scheme is p/(7i + 2(70 V 7i))-semi T^-stable with respect to Brownian motion. The proof of 
Corollarv 12.161 is thus completed. D 

Theorem l2.13[ Corollarv l2. 151 and Corollarv 12 . 161 give sufficient conditions for the Euler-Maruyama scheme to 
be a-semi ^/-stable with respect to Brownian motion with a £ (0,oo), V: M.'^' — > [0, 00) appropriate and d G N. 
One may ask whether the Euler-Maruyama scheme is also ^-stable with respect to Brownian motion. The 
next result, which is a corollary of Theorem 2.1 in |35j (which generalizes Theorem 2.1 in [36)). disproves this 
statement for one-dimensional SDEs in which at least one of the coefficients /i and a grows more than linearly. 
There is thus a large class of SDEs in which the Euler-Maruyama scheme is a-semi V^-stable but not F-stable 
with respect to Brownian motion. 

Corollary 2.17 (Disprove of V^-stability with respect to Brownian motion for the Euler-Maruyama scheme). 
Let T, p £ (0, 00), a, c G (1, 00), let fi^a: M — > K he Borel measurable functions with \fJ.{x)\ + |cr(a;)| > -^^ for all 
a; G M with \x\ > c and with supj,g]g |cr(a;)| 7^ 0. Then there exists no Borel measurable function y: M — > [0, cx)) 
which fulfills limsup^x q supj.gjj (i+v(x)) "^ °^ '^'^'^ 

E[V{x + fi{x) ■ t + a{x)Wt)] < eP* ■ V{x) (83) 

for all {x,t) G M X (0,T] where W^: [0, T] x O — > M is an arbitrary standard Brownian motion on a probability 
space {n,J',F). 

Proof of Corollary \2.17\ Suppose that F: K — > [0, 00) is a Borel measurable function which fulfills (|83| and 
limsup^N^o sup^gg (iiy/j.-)\ < 00. Then there exists a real number r G (0, 00) such that \x\^ < -^ (1 + V{x)) for 
all X G K. Next observe that the assumption sup^.^^ |o'(a;)| 7^ implies that there exists a real number xq G M 
with (t(xo) 7^ 0. Now define a sequence Y^ : {0, 1, . . . , N} x J7 ^- R, A^ G N, of stochastic processes by Fq^ = xq 
and r„^i = Y^ + ii{Y^)^ + a{Y^){W^^n+i)T/N - W^t/n) for all n G {0, 1, . . . , TV - 1}, TV G N. Corollary [U 
then implies 



1 <ifl + limsupE[V^(rj^)]) <\.{l + eP'^ ■V{xo))< 

\ N^oo J 

This contradicts to Theorem 2.1 in 35j. The proof of Corollarv 12. 171 is thus completed. D 



limsupE[||rj(rirj < i ( l + limsupE[V^(rj(,^)J )< 1(1 + 6"^ • V^(xo)) < 00. (84) 



2.2.2 Semi stability for tamed schemes 

In the previous subsection, semi V^-stability with respect to Brownian motion for the Euler-Maruyama scheme 
has been analyzed. In this subsection, semi TZ-stability with respect to Brownian motion for appropriately 
modified Euler-type methods is investigated. We begin with a general abstract result which shows that if two 
numerical schemes are close in some sense (see inequality ((85)) below for details) and if one of the two numerical 
schemes is a-semi 1^-stable with respect to Brownian motion with a G (0, 00), V : W^ — > [0,oo) appropriate and 
d G N, then the other scheme is a-semi l/-stable with respect to Brownian motion too. 

Lemma 2.18 (A comparison principle for semi 1^-stability with respect to Brownian motion). Let a,T £ (0, 00), 
c G [0,00), p G [1,00), d,TO G N, let (il, J",P) be a probability space, let W : [0,T] x fi ^- R™ be a standard 
Brownian motion, let $, $: R'^ x [0,r] x R" -)■ R'^ be Borel measurable functions and let V G Cp(R'', [0,oo)) be 
such that $ : W^ x [0, T] x R™ — >■ R'^ is a-semi V -stable with respect to Brownian motion and such that 

{E\^\^{x,t,Wt)-Hx,t,Wt)\\'']y'' <c-t-\V{x)\^''' (85) 

for all {x,t) G {{y,s) e R'^ x {0,T]: V{y) < s""}. Then l> : R'' x [0,T] x R™ -^ R'' is also a-semi V-stable with 
respect to Brownian motion. 
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Proof of Lemma \2.18l Inequality (j59p in Lemma 12.121 implies the existence of a real number c G [c + 1, oo) such 
that 



E 



V{^x,t,Wt)) 



E 



< E 



V{^x,t,Wt))'^ +2P\cfE[\V{<^{x,t,Wt))\^\mx,t,Wt)-^x,t,Wt)\ 
+ 2P\c\PE[\mx,t,Wt)-Hx,t,Wt)\\'' 
and Holder's inequality hence gives 



(86) 



E 



V{^x,t,Wt)) <E Vmx,t,Wt)) 



2P\ 



\c\''(E[v{<^{x,t,Wt)) 



(p-i) 



E 



mx,t,wt) - Hx,t,wt)f 



(87) 



+ 2P\c\PE \\<i>{x,t,Wt)~<i>ix,t,Wt)\ 



for all X e M''. Moreover, by assumption there exists a real number p e [0, oo) such that 



E 



V{^{x,t,Wt)) <eP^-V{x) 



for all {x,t) e {(y, s) G 



E 



(0,T] : V{y) < s""}. Putting §E\i and dHE]) into dHZ]) then results in 



V{^ix, t, Wt)) < e"* • V{x) + 2P \c\P (e^* -c-t- V{x) + {cf ■ t^ ■ V{x)) 

< (e"* +t{T+ If e"^ 2(P+i) \cfP ) • V{x) 

< exp( [p + (T + If e''^2(P+i) |c|'^] • t) • V{x) 

for all (a;,i) G {(y,s) G K'' x (0,r]: F(y) < s""}. The proof of Lemma [^T^ is thus completed. 



D 



A direct consequence of Lemma 12.181 is the next corollary. It proves a-semi ]/-stability with respect to 
Brownian motion with a G (0,oo), ]/: M'^ — > [0,oo) appropriate and d G N for a class of suitably "tamed" 
numerical methods. 

Corollary 2.19 (Semi V-stability with respect to Brownian motion for an increment-taming principle). Let 

a,T & (0,oo), I3,ce [0,oo), p G [l,oo), d,m G N, let {n,T,¥) be a probability space, let W : [0,T] x fl -)■ R"" 
be a standard Brownian motion, let V G Cp(R'*, [0,oo)), let <I>: K'' x [0,r] x M™ — )• K'' &e a Borel measurable 
function which is a-semi V-stable with respect to Brownian motion and assume that 



E 



\\<^{x, t,Wt)- x\M < c ■ tP(i-^) • V{x) 



(89) 



for all {x,t) G {{y,s) G R"* x (0,T]: V{y) < s""}. Then the function R'* x [0,r] x M™ 9 {x,t,y) ^ x + 
,, \^„'^,Z^\ FT G R'^ is a-semi V-stable with respect to Brownian motion. 



Proof of Corollary \2.19\ We show Corollary 12.191 by applying Lemma 12.181 We thus need to verify inequal- 
ity dSSD- For this let $ : M'* x [0, T] x R™ ^ R'' be given by $(i, x,y) ^ x + ^ax(i!^l[i'(l7,^)-.||) for aU x G W\ 
i G [0, T], J/ G R" and note that 



\'^(x,t,Wt)-'^{xX'^t^ 



^\\^{x,t,Wt)-x\\-- 
<tf'\\^{xXWt)-xf 



Lax(l,i''||$(a;,i,y) -a;||) - 1 
max(l,i/3||$(x,i,2/)-x||) 



(90) 



for ah X G R'^, i G [0,T] and all y G R™. Estimate ([HS]) therefore shows dSi]) and Lemma ETT51 hence completes 
the proof of Corollary [?J9l D 



The next result is an immediate consequence of Corollary [2T9] and Theorem l2.13l It shows a-semi ^-stability 
with respect to Brownian motion with a G (0, oo), y : R*^ — >■ [0, oo) appropriate and d G N for a suitable "tamed" 
Euler-Maruyama scheme. 
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Corollary 2.20 (Semi V^-stability with respect to Brownian motion for an increment-tamed Euler-Maruyama 
scheme). Let T G (0,00), d,m e N, p S [3,cx)), 0,70,71 S [0,oo), let /i: M"* ^ M'^, a-.R'^^ K'*^'" be Borel 



measurable functions and let V € Cp{M.'^, [1, 00)) with 



forallxe R". Then the funet^onR'ix[0,T]xR"^ 3 {x,t,y) ^ x+^^^^^^^igg^glfj^ e R'' is ,^+,.„4,„,^^,i/,) - 
semi V -stable with respect to Brownian motion. 

Proof of Corollary \2.2(A We prove Corollary 12.201 by using Theorem 12. 131 and Corollary 12.191 and need to verify 
condition dHU). Note that the assumptions ||/x(x)|| < c\V{x)\^'^'>+'^^/p and ||cr(a;)||L(R™,Rd) < c\V{x)p^+'^^/^^P^ 
for all a; e M'' imply 

¥.\Mx) ■ t + a(x)Wt\\^''] <(m+l)('P-^) E[|W^i(')|2p] (mx)]]^" t'P ^y^\\ak{x)\f'' tA 
^ ^ V fc=i / (91) 

< c^P (m + if^ V.[\W[^\^P] {\V{x)\^^'"'+^'^ t^P + \Vix)\'^'"+^^ iP) V{x) 

for all a; e M'' and all t G [0,r]. Applying Theorem 12. 131 and Corollary 12.191 therefore shows that $ is a-semi 
y-stable with respect to Brownian motion with 

a:==minf-^^,^^, -^ -V(0,oo). (92) 

V27o + l'7i + l 7i+2(7oV7i)y ^ ' ^ ' 



Next note that 



P 



^(70 + ^, 71 + 1,71 + 2(70 V 71)) 

p p y^"^) 



max(7o + 5, 71 + 2 max(7o, 71, 5)) 71 + 2 max(7o, 71, 5) 

The proof of Corollary 12.201 is thus completed. D 

Note that, under the assumptions of Theorem 12.131 the Euler-Maruyama scheme is p/(7i -I- 2max(7o,7i))- 
semi y-stable with respect to Brownian motion but the appropriate tamed numerical method in Corollary l2.20l 
is p/(7i -|- 2max(7o,7i, l/2))-semi V-stable with respect to Brownian motion. 

2.2.3 Moment bounds for an increment-tamed Euler-Maruyama scheme 

In Subsections 12. 2.1] and 12. 2."^ above, we established under suitable assumptions that both the Euler-Maruyama 
scheme (see Theorem 12. 131 in Subsection 12. 2. ip as well as appropriate tamed Euler methods (see Corollary l2.20l 
in Subsection 12.2.2^ are a-semi IZ-stable with respect to Brownian motion with a G (l,oo), V : R"^ ^^ [0,00) 
appropriate and d G N. Corollary 12.61 can then be applied to show that these approximations have uniformly 
bounded moments restricted to events whose probabilities converge to one with convergence order (a— 1). In the 
case of appropriately tamed numerical methods (see Subsection 12.2.21 for a few simple examples and Section 13.61 
below for more examples), it can even be shown that the approximations have uniformly bounded moments 
without restricting to a sequence of events whose probabilities converge to one sufhciently fast. In the next 
result this is illustrated in the case of the increment-tamed Euler-Maruyama scheme from Corollary 12.201 Its 
proof is based on an application of Corollarv l2. 91 above. 



Corollary 2.21 (Moment bounds for an increment-tamed Euler-Maruyama scheme). Let T G (0, 00), d, m G N, 
p G [3,00), c,7o,7i G [0,00), let ij:R'^-)-R'^,a:R'^-)- R'*^™ be Borel measurable functions, let (1^, J",P) be a 
probability space with a filtration {J^t)t<£lQ,T], ^^t W: [0,T] x Jl ^ M™ be a standard {Tt)te[o,T]- Brownian motion 
and let ^: fJ -J^ R** be an J^o/B{R'^) -measurable mapping. Moreover, let V G Cp{R'^, [l,oo)) with E[V{^)] < 00 
and with 

ig,,aV)ix)<c-Vix), \\ti{x)\\<c\V{x)\V-^], ||a(x)L(R..,R.) <c|F(a;)|[^] 

for all x G R''. Furthermore, let Y^ : [0,T] x 57 — > R'', iV G N, be a sequence of stochastic processes given by 



F(f = ^ and 



M(?^ ) ^ + <^(^D (Wi,.-,i,T - Wnl 



' N 

max(l, ^MY^)^ + a{Y^){W^^^^ - W^)\\) 



^t -^BT + {—-n)- ——jr^rTr-ToWTT^—TcyWTTTTr 777 — 7TTT v94j 



for all te (if,^^^^], n£ {0, 1, . . . , A^- 1}, iV G N. Then sup ,^^j^ sup te[t)^T]m\Yt^\\''] < 00 for all q e [0,oo) 
which satisfy q < 2^i+4max^(7o,7i,i/2) ~ ^ a"*^ sup^GK-^ lkllV^(a;) < 00. 
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Proof of Corollarv \2.21\ Let q E (0,oo) be a real number which satisfies q < 



p 1 r,„rl 

27i+4max(7o,7l,l/2) 2 

sup^gRd ||a;|j9/l/(a:) < oo. Next observe that Corollary l^?^ shows that the function $: M'' x [0,T] x M™ -^ W^ 
given by $(x, t,y)^x+ ^,^(i!tfu(t)t+a'(x)vl\) ^r aU x e R'^, t G [0,T], y G K™ is a-scmi IZ-stablc with respect 

7i+2max(7o,7i.i/2) ^ (0'°«)- Morcover, note that 



to Brownian motion with a — — 



N-l 



Ir.^ll < 



iieii + E 



<iieii + 



n=0 



max( 1, 






(95) 



for all i e [0, T] and all A^ G N and therefore 



sup sup (n^^'^WWY/^W'^I 
NeNte[a,T] ^ 



L°°(0;B 



< (X). 



The inequality 



hence shows 



l-a= -2 



a 1 
2" " 2 



= -2 



sup sup 

NeNte[o,T] 

Corollarv l2.9l can thus be applied to give 



271 +4max(7o,7i,i) 



< -29 



< 00. 



lim sup sup E 

N^oo ne{0,l,...,N} 



Combining this and (|96p finally implies 



sup sup E 

AfeNte[o,T] 



Yt 



Y, 



NUI 



N ll"? 



< 00. 



< 00. 



The proof of Corollary 12.211 is thus completed. 



(96) 

(97) 
(98) 

(99) 

(100) 
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Corollary 2.22 (Potencies of the Lyapunov-type function). Let T £ (0,oo), d,m E N, p E [3, oo), q E [l,oo), 



:;,7o,7i £ [0, oo), let fj.: 



]>dxr 



be Borel measurable functions and let V E CpiM."^, [l,oo)) 



with \\^l{x)\\ < c|F(x)|('^«+i)/P, |k(a;)||i(R™^Rd) < c |y(a;)|('^i+2)/(2p) and 

{Qt,,aV){x) + Y^vTx) <c-V{x) (101) 

for all X E M''. Moreover, let (J7, J^, P) be a probability space with a filtration {J't)te[a.T], l^t W : [0,T] x J7 — > 
M™ be a standard {J-t)tGlo,T]-Brownian motion, let £^: ft —>■ R'' be J- q/B{W^) -measurable with E[V^(^)] < oo 
and let Y^ : [0,T] x 51 — >■ W^ , N E N, be a sequence of stochastic processes given by Y(J^ = ^ and (j94p 
Then sup^gpj supfgrpTi E[||Fj^||''] < oo for all r E [0,oo) which satisfy r < -^ ^ 

SUP;:„gR<i \\x\Y/yi.x) < OO- 



+4max(7o,7l,l/2) 2 '^"'^ 



Proof of Corollary \K. 



Define the function V : 



[1, oo) by V{x) = (^(2;))' for all x eR'^ and observe that 



yix)\\ <c|f(a;)|[^] 



and 



\aix)\\ 



Tl+21 



i) < c\V{x)\^^^ 



(102) 



for aU x E M"*. Lemma [2Tl implies that V E C^q(R'', [1, oo)) and inequality dZZ]) yields that Gt,,aV < qcV. 
Therefore, Corollarv 12.211 shows that sup^gpjSup(g[o,T] ^[ll^^ll'^] < ^c for all r E [0,oo) which satisfy r < 
27i+4maxf7o.7i,i/2) ~ 5 ^^^ sup^gRd || a;|| ''/V" (x) < OO. The proof of Corollary H^ is thus completed. D 



We now illustrate the moment bounds of Corollarv 12.211 and of Corollarv l2.22l bv two simple corollaries. 



Corollary 2.23. Let c,T E (0, oo), d,m eN, let ^: 



pdx? 



be Borel measurable functions 



and let V EU 



pe(o,oo)C'p( 



[1, oo)) with limsup^N^o supj.gRd yhP) < °o '^^'^ ™*^^ 



sup 

xeTSL'' 



{Gi.,aV){x) , r\\V'{x)a{x)\\ 



V{x) 



\VixW 



< oo, sup 

xeM'' 



\^i{x)\\ + ||cr(a;)||L(R" 



(1 + 11^11^) 



< oo 
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for all r G [0,oo). Moreover, let (il, J^,P) be a probability space with a filtration {J^t)te[o.T]j ^^i ^ '■ [0,T'] x 
n — > R™ be a standard {J^t)t£[o,T\-Bro'wnian motion, let £,: 51 — > M.'^ be an J-q/B{W^) -measurable mapping with 
E[y(^)] < oo and let Y^ : [0,T] x V, —>■ R'*, N G N, be a sequence of stochastic processes given by Y(^ = ^ and 
((Ml) . Then sup^g^sup^gfo^j.] E[||yt^||«] < oo for all q e [0,oo). 

Corollary 12.231 is an immediate consequence of Corollary 12.221 Note also for every d E N that a function 
T^: R"* — > [l,oo) is in Upg(o_oo)C'p(K'^, [l,oo)) if and only if it is twice differentiable with a locally Lipschitz 
continuous second derivative and if there exists a real number c £ (0,oo) such that X]i=i 1!^ (2^)IIl(»)(e<',r) ^ 
c |F(a;)|[^"^''^l for Agd-almost all x € R"*. The next corollary of Corollarv l2.21l is the counterpart to Corollarv l2.16l 

Corollary 2.24 (A special polynomial like Lyapunov-type function). Let T E (0, oo), c, 70,71 E [0, 00), p E 
[3, 00), d,m eN, let fi:R'^^ R'^, a: R'^ ^ jjdxm ^^ Borel measurable functions with 

{x, ^^{x)) + l£^||a(x)||^5(j,„ R,) < c (l + \\xf) , (103) 

Mx)\\<c{l + \\x\\^'"'+'^) and \Hx)hiM-.m^)<c{l + \\xf-^^) (104) 

forallxEW^. Moreover, let {Vt,J- ,W') be a probability space with a filtration {J^t)te[o,T]! letW: [0,T] xf2 — > R'" be 
a standard {Tt)te[a.T]-Brownian motion, let ^: fi — > R'' be an J-q/B{M.'^) -measurable mapping with K[V{£^)] < 00 
and let Y : [0, T] x fl — )• R , N E N, be a sequence of stochastic processes given by Y^ — ^ and (|94p . Then we 
obtain sup^gp, sup^gfo j.] E[||Ft^||9] < 00 for all q E [0,oo) which satisfy q < 2^-1+4 max'(7o,7i, 1/2) " ^■ 

Proof of Corollary\K^ Let V': R"* ^ [1, 00) be given by V{x) = 1 + ||a;||^ for all x E R'^. Then note, as in the 
proof of Corollarv l2.16l that 

r 70+1 1 [21+11 

ig^,aV)ix)<2-p-c-Vix), \\^iix)\\<2c\V{x)\^^\, Mx)\\^^^^^^,^<2c\V{x)\^^^\ (105) 

for all X E R'^. Combining this, the fact V E C^(R'', [1, 00)) and Corollary H^ then shows that 

sup sup E[||rt^||«] <oo (106) 

NeNte[o,T] 

for all g E [0,oo) which satisfy g < - — ^^^ . IT~5 and sup^^ggd ||a;|l'/y(a;) < 00. The proof of Corollarv l2.24l 

is thus completed. D 

2.3 Implicit approximation schemes 

In this subsection, stability properties and moment bounds for implicit approximation schemes are analyzed. 
The main results of this subsection are Corollarv l2.27l for the fully drift-implicit Euler scheme and Lemma [2.281 
for partially drift-implicit approximation schemes. 

2.3.1 Fully drift-implicit approximation schemes 

Corollary 12.271 below proves uniform bounds on the g-th moments of fully drift-implicit Euler approximations 
for a class of SDKs with globally one-sided Lipschitz continuous drift coefficients where q E [0, 00). This result 
generalizes, in the case of the fully drift-implicit Euler scheme, Theorem 3.6 of Mao & Szpurch [55] which 
establishes uniform bounds on the second moments of the numerical approximation processes. First we prove 
two auxiliary lemmas (Lemma 12.251 and Lemma 12.261) . The first lemma is a slight generalization of Lemma 3.2 
of Mao & Szpruch 52J. 

Lemma 2.25. Let d E N, c E (0, 00) and let fi: R"^ ^ R"^ be a function with {x,fi{x)) < c (l + ||x|p) for all 
xER'K Then 

l + \\x- fi{x)sf < e^"^*-") (1 + Ik - Kx)tf) (107) 

for all X e R'' and all s,t E [0, ^] with s < t. 

Proof of Lemma \2.25\ Throughout this proof, let F : [0,oo) x R'' — > R'^ be a function defined through Ft{x) := 
X — fj.{x)t for all {t, x) E [0, 00) x R*^. The assumption {x, n{x)) < c (l + ||a;|p) for all a; e R*^ then implies 

1 + ||F,(a:)f = 1 + \\Ft{x) + {t - s)^i{x)\\' 

= 1 + \\Ftix)f + 2{x ~ tii{x), {t - s)ii{x)) + {t- sf y{x)f 
= 1 + \mx)f + 2{t s){x,^,{x)) -{t + s){t s) Mx)\\' 
<l + \\F,{x)f-^2c{t-s){l-^\\xr) 
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for all a; G K*^ and all s, i G [0, oo) with s < t. The special case s = in (|108p shows 



l + ll^ll^^ il-2tc) 
for all a; G K'^ and alH G [0, ^). Finally, we apply pTO)) to pII5)) and arrive at 

1 + \\Fsix)f < 1 + \mx)f + 2c (t - .) ^^+™^^ 
= (l + ||F,(a:)||^) (l + f^) < (l + ||F,(x)||^) 

for all a; G M'' and all s, i G [0, ^) with s < t. 



(109) 



g(l-2ct) 



(110) 



n 



Lemma 2.26 (Stability of the fully drift-implicit Euler scheme). Let d,m G N, c £ (0,oo), p G [2,oo) and let 



/i: M'* -^ R'^, cr: M'' -> R'^^™ be functions with 



r:,fi{x)) 



(p-i) 



2 i|C^(a^)liHS(R'",R<') ^c(1 + ||.t|| 



for all X G M'' . T/ien there exists a real number p G R such that 



E 



{l + \\x + a{x)Wtfy^ < e"* (l + \\x - ^^(2^)^11")' 



(111) 



(112) 



for all X G R'^, t G [0, ^] and all q G [0, |] where W : [0, cxd) x fJ — > R™ is an arbitrary standard Brownian motion 
on a probability space (f2, J-", P). 

Proof of Lemma \2.2b\ Throughout this proof, let F : [0,cx)) x R'' — ^ R'^ be a function defined through Ft{x) := 
X - n{x)t for aU (i, x) G [0, oo) x R'' and let e^' := (1, 0, . . . , 0), . . . , e^^ := (0, . . . , 0, 1) G R™ be the canonical 
basis of R™. Ito's lemma yields 



E 



l + ||a; + cr(a;)T^t|| 



= (1 + ||.t|H +q\W{x'"^ 



HSl 



{l + \\x + a{x)W.fy''''> 
q{2q-2)J2 r^{^ + \\x + 'y{x)W,fY''-''' 



ds 



\{x + a{x)W,,aix)eTW 



fc=i 



<[l + \\xfy +qi2q-l)\\aix)\'' 



ds 



(113) 



_f/S(K'",R'') 



E 



l + ||x + a(a;)M^,||- 



(9-1)' 



ds 



for all t G [0, oo), x G R"* and all q G [1, oo). In addition, the fundamental theorem of calculus implies 

(l + \\Ft{x)\\'y = (l + ||x||')' -2qJ (l + ||F,(a;)||')^''"'^ (F,(x), ^.(x)) ds (114) 

for all i G [0, oo), X G R"* and aU q G [1, oo). Combining (|113p and (I114p gives 



< 



I + \\x + a{x)Wt\ 
(l + \mx)\\'y +2q ni + mix)f) ''^'' {F,ix),^,ix)) 



(«-i) 



+ q{2q-l)\\aix)\\ 



HS( 



E 



I + \\x + a{x)Ws\\ 



ds 



(115) 



ds 



for all t G [0, oo), a; G R'' and all q G [l,oo). Roughly speaking, we now use (|115p to prove (|112p by induction 
on g G [0, |]. More precisely, let k: [0, |] — >■ [0,oo) be a function defined recursively through K{q) := 6cq for all 

g G [0, 1] and through K{q) := 7p^ (c + 1)^ exp { '^%~^'' + p) for aU g G (n, n + 1] n [0, |] and aU n G N. We then 
prove 

E[(l + ||x + a(a;)m||')'] <e«('')*(l + ||^t(x)||')'' (116) 
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for all t G [0, j^], X G M.'^ and all q E [n, n + 1] n [0, |] by induction on n £ Nq. For the case n = and g = 1, we 
apply assumption (jllip and Lemma 12.251 to the right-hand side of (|115l) and get 



(117) 



E 1 + \\x + a{x)Wt\\ 

< 1 + ll^t(a^)ll' + J {^{x,fiix)) + |k(a:)||^s(R„ R.) - 2s||/i(a:)f ) ds 

< 1 + ||Ft(a;)||' +2c f (l + ||a;f ) ds < 1 + \\Ftix)f + 2tce^'' (l + \\Ftix)\\ 

< (1 + 2ect) (l + ||^*(:r)||') < e^-* (l + \\Ft{x)f) < e'^«* (l + ||F,(a;)f ) 

for all t S [0, j^] and all x E M.'^. Next observe for every q E (0, 1) that the function [0, oo) 3 z i-^ z'' E [0, oo) is 
concave. Hence, Jensen's inequality and (J117p imply 



E 
< 



l + \\x + (7{x)Wt\\ 



< E 



l + ||a; + (7(2:)W^t||' 



(e'^W* (l + ||F,(a:)||^))' = e'^^^)'* (l + |iF,(x)f )' 



(118) 



for all t E [0, j^], X E M.'^ and all q E [0, 1]. This proves (|116l) in the base case n = 0. For the induction step 
n — )■ n + 1, apply the induction hypothesis on the right-hand side of (jllSp to obtain 



E 



[l + \\x + a{x)W,fy]-{l + \mx)f) 



< 



,(9-1) r 



(l + \\F,ix)f) [2q {F,ix),^i{x)) + q {2q - 1) e^^-'''^ \\aix)\\ls^^^^^,^ 

2qJ {l + \\Fs{x[ 



(9-1) 



,(9-1) 



(29-1) 



^x,^j.{x)) H 2 I1'''(^)IIhs(r™,r<*) 



+ 2q (l + \\F,{x)f) ' [iH^ (e«('-i)« - l) \\a{x)\\l,^^^^^.^ - . ||^(x)|| 



ds 



ds 



and assumption (jllip and the fact that the function (0, oo) 3 x i-^ — E (0, oo) is increasing hence show 



E 



{l + \\x + a{x)Wtfy]-{l + \\Ft{x)f) 



<2qc / l + ||F,(a:) 



(9-1) 




/■* / ,\(9-l) 

2qJ [l + \\F,ix)f 



i + \\x\nds 



i1-h] 



--^\\a{x)\'' 



HS{Ti 



IM-^)II 



ds 



(119) 



<2qc (l + ||F,(a;)f) ' {l + \\x\\^) ds 



2q / l + ||F,(x) 



(9-1) 



2cpexp(^^)||a(a;)|l 



2 
HSC 



M^W 



ds 



for all t E [0, ^], a; e M'' and aU q E [n -I- 1 , ri -I- 2] n [0, |]. Next observe that Young's inequality and again 
assumption (jllll) give 



r\W{x)\\HSiR'-M'i)- llM(a^)ll 



2r 



(p-1) 



f:,^^(^^))+(£_ll\l^^:,)\\ 



HSiR'^.R'^) 



^ 2.C (1 + lNP) ^ 2. NIHMx)|| „ 1,^(^)1,2^ 



< 



2rc r^ 

(P - 1) ^ (p - 1)' 



b-1) 



2rc(l + ||x|p) r^ljxf 
(p-1) ^(p-1)' 



(120) 



(l + llxf) < (2rc + r2)(l + ||.x||2) 
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for all a; e R'^ and all r e [0, oo). Combining (I119p and (|120l) implies 



E 



l + ||a; + a(a;)W^,||')'] - {l + \\Ft{x)fy 



(9-1) 



(l + \\xf)ds 



<2qc f (l + ||F,(a;) 
+ 2qj' (l + \\Fs{x)fy'''^ ■ s ■ {4c^p + 4cV) exp(^^) (l + \\x\\') ds 
and Lemma [2.251 hence gives 



(121) 



E 



l + \\x + a{x)Wt\\ 



< 

< 
< 



,(9-1) 



(122) 



2qc{l + 6cp'e^p[^<^)) j\l + \\F,ix)fy' {l+\\xr)ds 
pc(l + 6cp2exp(^^)) f\iMt-sK,-r)+4ct) (^i + \\Ft{x)fy ds 
7ip3 (c+ l)^xp(%i) +p) (l + ||F,(x)f )' 
for all t e[0,j^], X eR'^ and all q e [n + 1, n + 2] n [0, |]. The estimate 1 + r < e'' for all r e M therefore shows 



I + \\X + (7{x)Wt\\ 



+ t7p' (c+ l)2exp(^^ +p)] (l + \\Ft{x)f) ■ 
<e*-(9)(i + ||Ft(a;)f)' 



for alH e [0, ^], x e M'' and all q G [n + 1, n + 2] n [0, |]. This finishes the induction step and inequality (J116p 
thus holds for all n E Nq. In particular, we get from inequality (|116l) that 



E 



l + \\x + a{x)Wt\\'^j < expft • sup,,g[o £] K(r)j (l + \\x - fi{x)t\\^] 



(123) 



for all t E [0, j^], X E M.'^ and all q E [0, |]. This and the estimate sup^^rp pi K(r) < cx) then complete the proof 
of Lemma [123 ' □ 



Now we apply Lemma 12.261 and Corollarv l2.2l to obtain moment bounds for fully drift-implicit Euler approx- 
imations. 

Corollary 2.27. Let d,m E N, c,T E {0,oo), p E [2, oo), let (^Q,J^,{J^t)te[a,T]jP) be a filtered probability space, 
let W: [0,r] X O — > M™ be an {Tt)te[o.T]-Brownian motion, let £_: ft ^ M.'^ be an Ft^ / B{K^)-measurable function 
and let PL : W^ -^W^, a:W^ ^ R''^™ be Borel measurable functions with E[||^(OP] < oo and 



{x - y, fi{x) - n{y)) < c ||a; - y|| 



(^,^(^)) + (£_l)||^(^)|j 



HS(R'",R'') 



<c l + llxll 



(124) 
(125) 



for all x,y E R'^. Then there exists a unique family Y^ ; {0, 1, . . . , N} x il ^ R'^, N E NO {cT, oo), of stochastic 
processes satisfying Yfl^ — ^ and 

>;^l = y,^ + Kyn+l)w + <y{yn)iWin+l)T/N - W„t/n) (126) 

for all n E {0, 1, . . . , A^ — 1} and all N E N D {cT, oo) and there exists a real number p E (0, oo) such that 



lim sup sup E 

N^oa ne{0,l.....,N} 



N\\2^<l 



{1+ lie in 



< e"^ • E 



{l + lieff 



(127) 



for all qE [0,f]. 



Proof of Corollary \2.27\ Throughout this proof, let F: [0, oo) 
x — ii{x)t for all {t,x) E [0,oo) x R*^. Next observe that 



l"^ be a function defined through Ft{x) 



{x ~ y, {fi{x)t - x) - {fi{y)t - y)) < {ct - 1) \\x - yf 



(128) 
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for all x,y E M.'^ and all t E [0,oo). This inequality ensures the unique existence of the stochastic processes 
Y^: {0,l,...,N} xn^R'^, N eNn {cT, cx)). In the next step, note that 



\FT/N{Yn+l] 



K + <^iyr^){Win+l)T/N " W'nT/Tv) || 



(129) 



for all n £ {0, 1, . . . , A^ — 1} and all A^ G N n {cT, oo). Lemma [2.261 hence implies the existence of a real number 
p e M such that 

e[{1 + \\Fr/N{Y^\^)rr I Y^^] < exp(f ) • {l + ||fT/A'(r„^)|p}^ (130) 

P-a.s. for all n e {0, 1, . . . , A^ - 1}, A^ e N n [4cT, oo) and all q e [0, §]. Next fix a real number q £ [0, f ] and 
we now prove p?7)) for this q e [0, f ]. If E[||^|p9] = oo, then (fT?7)) is trivial. So we assume E[||^||29] < oo for 
the rest of this proof. Hence, we obtain that E[||^|p'' + ||/^(^)|P'] < oo. Now we apply Corollary 12.21 with the 
Lyapunov-type function F : M"^ — >• [0, oo) given by V{x) = {l + ||-FT/iv(2:)|P}* for aU x e M'^, with the truncation 
function ( : [0, oo) — >■ (0, oo] given by ({t) = oo for all t G [0, oo) and with the sequence i„ G K, n € No, given by 
tn = rtun{nT/N, T) for all n e No to obtain 



sup E 
T.e{o,i N} 



{i+\\FT,N{Y^)r] 



< e''^ • E 



219 



{1 + II^T/^(0II'} 



(131) 



for all n £ {0, 1, ... , A^} and all A^ G N n [4cT, oo). Lemma [2.251 and the dominated convergence theorem hence 
give 



lim sup sup I 

N^oa ne{Q,l.....N} 

= e^"^ • E 



{i + r,riP} 



21 9 



N 



lim {l + \\FT,Nm^] 



< lim 



= 6"^ 



/ 4cT 

(e « 



=PT 



E 



2-19 



{i + \\FT/Nmn 



219 



{1 + iieir} 



and this completes the proof of Corollary 12. 271 



(132) 
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2.3.2 Partially drift-implicit approximation schemes 

In Subsection 12.3. II above, moment bounds for the fully drift-implicit Euler schemes have been established. This 
subsection concentrates on partially drift-implicit schemes. 

Lemma 2.28 (Partially drift-implicit schemes for SDEs with at most linearly growing diffusion coefficients). Let 
h £ (0, oo), c G [1, oo) be real numbers with he < j, let (O, T, P) be a probability space, let W : [0, oo) x fi — > M™ 

— > ]R<^^™ be Borel measurable functions with 



be a standard Brownian motion, let f. M. x R'^ ^ M."- , a : M. - 

(y, ip{x,y)) <c{2+\\x\\^ + \\y\\^) , lk(x)||^5(„„ „,) < c (l + \\x\\^) 

for all X, y £ R^ and let Y : Nq x fl ^ R^ be a stochastic process with 

Yn+l = Yn + ^(r„, r„+l) h + a{Yn) {W^n+l}h " Wnh) 

for all n G No . Then 



sup IJYfcl 
fce{o,i,...,n} 



Lp{n-f 



< 2 (l + l|i"o|liP(n;R<i)) exp((4-H/xp/2)c'n/i) 



(133) 
(134) 

(135) 



for all p £ [4, oo) and all ri G No- 

Proof of Lemma \2.28[ First of all, let p £ [4, oo) be arbitrary. If E[||yoP] = oo, then inequality (|135[) is trivial. 
Thus we assume E[||yoP] < oo for the rest of this proof. In the sequel, we will show (I135p by an application of 
Corollary 12.31 For this we define a stochastic process Z : N x O — > M'* through 



Zn ■■^ ,^ ^^^. (2(r„_i,a(y„_i)(iy„;, - W^(„-i)0) 



(1 - 2hc) 

+ \\<7[Yn-i){w^h - w^(n-i)/.)f - h lk(y„_i)|| 



HS( 



(136) 



for all n G N. Then observe that Z„ is integrable and that E[Z„ | (2'fe)j,£{o,i,...,n-i}] = P-a.s. for all n G N. 
It thus remains to verify inequality ^\Ff\ and inequality ([^51) in order to complete the proof of Lemma [2. 281 For 
inequality ((27|) note that Young's inequality and assumption (I133P give 



2 \\Y, 



- 2 (r„, r„_i + a(r„-i)(w„ft - ^^(n-i)/.)) + 2/j {Yn. (^(r„-i, y„)) 

< ||y„||' + ||y„_i + (T{Yn-i){Wr,h - M^(„-i)h) f + 2/ic (2 + ||r„-i|p + ||r„||2) 



(137) 
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for all n E N. Rearranging (|137p yields 



1 + \\Yn 



< 



1 + ||r„_i + a{Yn-i){Wnh ~ M/(„-i),0 f + 2hc (1 + ||r„-i|p) 



i + ||y„_if + /i||a(y„_i^"' 



(1 - 2hc) 



(1 - 2hc) 

fi + ||r„_i||2)(i + 3/ic) 7, , „ „„, 



'^n 



(138) 



(1 - 2hc) 

for all n G N where the last inequality follows from the estimate jirff < e^^ for all x G [0, -j]. Estimate (|138p 
is inequality (l?fl) with p ~ 7c G [0,oo). In the next step, we verify inequality ([^5]). Observe that the estimate 
||X — E[X]||^g,j^.g-) < 2 ||X||^5(Q.gj for all q G [l,oo) and all J^/;B(R)-measurable mappings X: fJ — > K and 
Lemma 7.7 in Da Prato & Zabzcyk [TB] give 



\Zn\ 



LP/2(n;i 



< 



< 



< 



< 



\{Yn-l,aiY^.l)iWnh-W^n-l)h))\ 



LP/^ifl;] 



+ \\cr{Yn-l){Wnh-W^„^l)h)\ 



LP(n;R<i) 













(1/2 


-/ic) 








n/8 


KK-i, 


a(r„_ 


-i)(-)>IIl./= 


(n;ffS(R",K)) 


1 P(P-1)'1 1 
^ 2 1 


^{Yn-l)\\LP{n-HS(VL"' 


,R'')) 












{1/2- h 


c) 








pVTT 

72 


r«-i 


llk(y 


»-l) l//S(R'' 


\R'i) 


Lp/'^(n-M 


^ + chp{p - 


-1) 


1 ] \\\^ I|2 


P/2(n;R) 












(1 - 2/?,c) 








pVh 
V2 


2 


2 

- + § 


> + \\Yn-i 


||21 


LP/2(0;K) 


+ c/ip(p - 


1) 


1 -X- \W l|2 

1 + 11^^.1-1 II Lv 


/2(n;R) 



and therefore 



l^nl 



LP/2(0;1 



pyhi 



< 



1 I II V l|2|| 

^2 11-^ + ll-fn-lll |lLP/2(n;] 



(1 - 2hc) 



■chp{p-l)\\l+\\Yr,-i\ 



lLP/2(n:] 



pc 



Vh 



-l= + Vh{p-l) 



V2 



(1 - 2hc) 



(1 - 2/ic) 

ll+ll^'n-lfl 



(139) 



LP/2(n;] 



<p{p + V2~ i)cVh II 1 + ||y„_i| 



211 <cp2V2ft||i + |jy„_i||2i 



Ilp/2(0;K) - -^ v-.-ll-L I n^n-iii |iip/2(n.R) 

for all n € N. This implies inequality (l28t with i/: N — > [0,oo) given by t-n = cp^\/2h for all n e N. Now we 
apply CoroUarv 12.31 with p — 7c, with the Lyapunov-type function V: R"* — ?► [0,oo) given by V{x) = 1 + ||x|p 
for all X e W^, with the truncation function (: [0, oo) — > [0, oo] given by (^{t) — oo for all t £ [0, oo) and with the 
sequence i„ G M, n G Nq, given by i„ = nh for all n G No to obtain 



sup ||Yfe| 

feG{0,l, ■■■,"} 



< e 



7cnh 



LP/2(n;I 



sup e-^^^-Mi + ||yfe||' 

fce{o, !,■■■,"} ^ 



LP/2(n;] 



< 2e 



7cnh 



1 + ^0 



LP/2(n:] 



exp X£ 



< 2 



,fc=i 

(l + ll>'o|lL(0;K^)) exp((7 + 2p\E)c^nh 



for all n G Mo . This finishes the proof of Lemma 12.281 



(140) 
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Lemma 2.29 (A class of linear implicit schemes for one-dimensional SDEs). Let c G [0,oo), T E (0,oo), let 
(i7, J^, (J-i)tg[o.T]i IP) be a filtered probability space, let W: [0,T] x J7 — >■ M 6e a standard {J^t)telo,T]-Brownian 
motion, let S,'- ^ ^ ^ be an To/B{M.) -measurable function and let a, b, a: R — ^ R be Borel measurable functions 
with 



■{a{x)x + b(x)) + \\a{x)\'^ <c{l + x^), a{x) < c, \b{x)\'^ < c {l + x^) 



(141) 
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for all X E R. Then there exists a unique family Y^ : {0, 1, . . • , N} x il — > 
processes satisfying Y^ — ^ and 

for all n e {0, 1, . . . , N} , iV G N H {cT, cx)) and we have 



iV e N n (cT, oo), o/ stochastic 

(142) 



sup sup E 

A''GNn[2cT,oo) rie{0,l,...,Ar} 



1+ y: 



TVI 



<e(«=+2)^.E[l + |e|' 



(143) 



Proof of Lemma \2.29[ First, note that the assumption a{x) < c G [0, cx)) for all a; e R (see (J141I) ') ensures the 
unique existence of a family Y : {0, 1, . . . , N} xil —?' R, N E NO {cT, oo), of stochastic processes satisfying Yq — £^ 
and (I142P . Next observe that (I14ip implies 



{x + b{x)tf + \aix)ft 
(1 - a(x)tf 



<x 



{x + bix)t + a{x)WtY 
{l-a{x)tf 

x^ (1 - 2a{x)t) 2i {x ia{x) x + h{x)) + \ \a{x)\'] + \b{x)\' t^ 
{l-a{x)tf {l~a{x)tf 

2tc (l + x^) + c (l + a;2) t2 



(144) 



and therefore 



E 



<{l + x' 



{l-a{x)tY 



{x + b{x)t + (j{x)WtY 
(l-a{x)tf 

(2c + ct)t 



\{l-a{x)t) 



(1 -a{x)t) 
for all t e [0, -) and all x G M. Hence, we obtain 

{x + b{x)t + cf{x)Wtf 



,,,(l..^)cxpi?^±^ 



(145) 



E 



{l-a{x)tY 



{l-ctY 



< (1 + x^) e(«^+2)* 



(146) 



for all t G [0, ^] and all a; e M. Note that (jl46p shows that the linear implicit scheme M x [0, ^] x M 9 
{x, t, y) 1-^ ^ n- r w^ G M is (1 + a;^)a;GR-stable with respect to Brownian motion. Moreover, combining (jl46p 



(l-a(x)t) 

and CoroUarv 12.21 results in (J143I) and this completes the proof of Lemma [2.291 
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3 Convergence properties of approximation processes for SDEs 

With the integrability properties of Section [5] at hand, we now prove convergence in probability, strong con- 
vergence and weak convergence results for numerical approximation processes for SDEs. For this, our central 
assumption on the numerical approximation method is a certain consistency property in Definition 13.11 below . 

3.1 Setting and assumptions 

Let T e (0,oo), d, m e N, let (f2,J^, P) be a probability space with a normal filtration {Ft)t£[o,T] and let 
H^: [0,r] X ri — )■ M™ be a standard (J^t)tg[o_T]-Browman motion. Moreover, let D C M'' be an open set, let 
fi: D ^)- M.''- and a: D ^)- K''^™ be locally Lipschitz continuous functions and let X : [O.T] x fl ^>- D he an 
(J-t)tg[o,T]-a'dapted stochastic process with continuous sample paths satisfying 



Xt = Xo+ f t4X,)ds+ f c7{Xs)dWs 
Jo Jo 



(147) 



-a.s. for all t e [0, T]. Note that we assume existence of a solution process staying in the open set D. Next let 
: M'^ X [0, T] X K" ^ R'' be a Borel measurable function and let Y^ : [0, T] x fi ^ M'^, A^ e N, be a sequence 



of stochastic processes defined through Yg^ := Xq and 



Y-:=Y^ + {if-n).^Y-^,W,,^ 



Wr. 



(148) 



for all t S (^, \ ] , n £ {0, 1, . . . , A^ — 1} and all A^ e N. A central goal of this section is to give sufficient 
conditions to ensure that the stochastic processes Y^ : [0, T] x fi -> R'^, A^ G N, converge in a suitable sense to 
the solution process X: [0, T] x f) ^ £> of the SDE pi7|) . 
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3.2 Consistency 

This subsection introduces a consistency property of the increment function 0: M'' x [0.,T] x M'" -^ M'' from 
Subsection [nm which ensures that the stochastic processes F^ : [O.T] x Q ^ W^, N £ N, defined in (|148p 
converge in probabihty to the solution process X : [0, T] x J7 — )■ _D of the SDE (jl47p (see Theorem 13.31 below) . 

Definition 3.1 (Consistency of numerical methods for SDEs driven by standard Brownian motions). Let T E 
(0,oo), d,m e N, let D C R"^ be an open set and let fi: D ^ R"^ and a: D ^ jjdxm ^g functions. A Borel 
measurable function cj): M"^ x [0,r] x M™ — > M.'^ is then said to be {fi, a) -consistent with respect to Brownian 
motion if 



lim sup —1= ■ sup E 



a{x)Wt - (l)ix,t,Wt) =0 (149) 



and 



lim sup (sup y{x) - i •E[(/)(x,t,W^t)]|| ) =0 (150) 

for all compact sets K C D where W : [0, T] x J7 — > R™ is an arbitrary standard Brownian motion on a probability 
space (rj, J",P). 

Note that (|149p in Definition 13.11 assures that the expectation in (I150p is well-defined. Further consistency 
notions for numerical approximation schemes for SDEs and results for such schemes in the case of SDEs with 
globally Lipschitz continuous coefficients can, e.g., be found in Section 9.6 of Kloeden & Platen [131, in Chapter 1 
of Milstein [56] . in Beyn & Kruse [2, in Kruse |44] and in the references therein. The next lemma gives a simple 
characterization of (/i, cr)-consistency with respect to Brownian motion. Its proof is straightforward and therefore 
omitted. 

Lemma 3.2. Let T e (0,oo), d,m £ N, let D C W^ be an open set and let fj.: D ^ R'^ and a: D -^ ^d-Km 
be functions. A Borel measurable function 4>: K'' x [0,T] x M"* — >■ W^ is then {fx, a) -consistent with respect to 
Brownian motion if and only if 

limsup(^-sup,gc^E[||a(a;)M^t-0(a;,i,W^t)ll]) (151) 

and 

limsupU- sup |U(a;)t-E[0(x,t,H/t)]|| ) ==0 (152) 

for all V &N where the sets Dy C D, v eN, are given by Dy := {v & D : \\x\\ < v, dist{x, D'^) > i} for all v eN 
and where W: [0,T] x f2 — ^ K™ is an arbitrary standard Brownian motion on a probability space (17, J^, P). 

A list of numerical schemes that are (/x, (T)-consistent with respect to Brownian motion can be found in 
Section! 



3.3 Convergence in probability 

The next theorem shows that the stochastic processes F^ : [0,r] x fJ — > R"^, A^ e N, in (|148p converge in 
probability to the solution process X: [0, T] x f7 ^ £> of the SDE (|147p if the increment function </) : M'' x [0, T] x 
K"* —> R"* is (/i, CT)-consistent with respect to Brownian motion (see Definition 13. ip . 

Theorem 3.3 (Convergence in probability). Assume that the setting in Subsection \3.1\ is fulfilled and that 
(j): R^ X [0,r] X M™ — > R'^ is {fi, a) -consistent with respect to Brownian motion. Then 



lim P 

N-yoo 



sup llXt-rt^ll >e 
te[o,T] 



= (153) 



for all £ e (0, oo). 



The proof of Theorem l3.3l is presented in the following subsection. More results on convergence in probability 
and pathwise convergence of numerical approximation processes for SDEs with non-globally Lipschitz continuous 
coefficients can, e.g., be found in [45l [221 [211 ISll [59l [37l [11] and in the references therein. 

Thcorcm l3 . 3l proves convergence in probability for a class of one-step numerical approximation processes (jl48p 
in the case of finite dimensional SDEs with locally Lipschitz continuous coefficients /j, and a. The locally Lipschitz 
assumptions on /i and a ensure that solutions of the SDE (|147p are unique up to indistinguishability. We expect 
that it is possible to generalize Theorem 13.31 to a more general class of possibly infinite dimensional SDEs and 
also to replace the locally Lipschitz assumptions on /x and a by weaker conditions such as local monotonicity 
(see, e.g., Krylov |45], Gyongy & Krylov [22] in the finite dimensional case and Liu & Rockner |48l [49] in the 
infinite dimensional case) which ensure that solutions of the considered SDE are unique. 
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Proof of Theorem | 

For the proof of Theorem 13.31 we introduce more notation. First, we define mappings Y^^ : il -^ W^, n E 

'-nT/N 

T>dxm 



{0,1,..., TV}, N e N, through Y^ := Fj^/j^ for aU n e {0,1,..., iV} and aU N e N. In addition, we use 




the mappings p.: R'' ^ W^, a = iaij{x))^^^^ d}je{i m} ■ ^^ ^ ^ and a,:W-^ W, i G {l,2,...,m}, 
defined by p.{x) := fj.{x), (t(x) := cr{x) for all a; e Z?, by fl{x) := 0, ct(x) := for all x € D'^ and by ai{x) :— 
{al^^{x), . . . , ad,^{x)) for aU a; e R'' and aU i e {1, 2, . . . , to}. Next the assumption on </«: M'' x [0, T] x M™ -^ R'^ 
to be (/i, CT)-consistent with respect to Brownian motion and Lemma 13.21 ensure that there exists a sequence 
ty e (0, T], w e N, of real numbers and a sequence D^ C D, v £ N, oi nondecreasing relatively compact open sets 
with UveNDi, = D such that E[||(?!)(a;,i, Wt)|j] < oo for ah x e A), t e [0,ty] and aU w e N and such that (|149p 
and (I150p are fulfilled (see Definition 13. ip . Then let c^, e [0,oo), u G N, be a family of real numbers defined by 

(154) 

for all w e N. Due to Lebesgue's number lemma one can indeed show that c^ < oo for all u e N since fj,: D — > M'' 
and a: D ^f j^dxm ^-^^ assumed to be locally Lipschitz continuous and since D^ C -D, t" G N, is a sequence of 
compact sets. 

Roughly speaking, the consistency condition of Definition 13 . 1 1 reg uires the increment function (p: M'^ x [0, T] x 
M™ -^ W^ from Subsection l3.1l to be close to the increment function of the respective Euler-Maruyama approxima- 
tion method. For this reason, we estimate the distance of the exact solution X : [0, T] x fJ — ?► I? of the SDE (jl47p 
and of the Euler-Maruyama approximations (see Lemma 13.41 below) and we estimate the distance of the Eulcr- 
Maruyama approximations and of the numerical approximations Y^^ : il-^-M'^, ne {0,1,..., N}, A^ e N, (see 
Lemma [3.81 below). The triangle inequality will then yield an estimate for HXtit — Y^\\, n e {0, 1, . . . ,iV}, 
N G N (see CoroUarv 13.91 below) . For this strategy, we now introduce suitable Euler-Maruyama approximations 
for the SDE (fTT7)l . More formally, let Z^ : n ^ R'', n e {0,1, . . . , N}, iV e N, be defined recursively through 
Zl^ ■.= Xq and 

Z^^, := Z.^ + mZ^) ■ ^ + o{Z«) (m/(^^ - W^) (155) 

for all n e {0, 1, . . . , iV - 1} and all A^ e N. Furthermore, let Z^ : [0, T] x fi ^ M'*, A^ e N, be given by 

Zf = Z,f + /l(Z,f ) • (t - f^) + a{Z^) (W^ - W^) (156) 

for all t e [^, ^"^^^^ ], n e {0, 1, . . . , A^ - 1} and aU A^ G N. Finally, let t^ : 9. ^ [0,r], w,7V e N, and 
(5^ : r^ ^ {0, 1, . . . , TV}, u, TV e N, be defined by 

rf (w) ~ inf ({T} U {t e [0,r] : Xt{u) i D,} U {t e [0,T] : Zf (w) i i?,}) (157) 

and by 

5^{uj) := min({A^}u{n G {0,1, . . . ,N}: Z^ [uj) i D,] \J [n G {0, 1, . . . , A^} : y„^(c^) i D„}) 



for all cj G r2 and all w, A^ G N. Using the notation introduced above, the proof of Theorem 13.31 is divided into 
the following lemmas. 

Lemma 3.4. Assume that the setting in Subsection \3.1\ is fulfilled, that the setting in the beginning of Subsec- 
tion \3.3\ is fulfilled and that (p: R'' x [0,7"] x R™ — >■ R"* is (ij^, a) -consistent with respect to Brownian motion. 
Then 



lim I 

Af-yoo 



sup ||Xt-Zf II > e 

t€[0,T] 



(158) 



for all £ G (0, oo). 

The proof of Lemma [5^ is literally the same as the proof of Corollary 2.6 of Gyongy & Krylov [52] (replace 
assumption (ii) in [22 by the weaker assumption of the existence of an exact solution; see also Section 2 in [39|). 
The proof of Lemma 13.41 is therefore omitted. 

Lemma 3.5. Assume that the setting in Subsection \3.1\ is fulfilled, that the setting in the beginning of Subsec- 
tion \3.3\ is fulfilled and that (f>: R'^ x [0,T] x R™ — > R'' is {fi, a) -consistent with respect to Brownian motion. 
Then 

lim limsupP[Tf < T] = 0. (159) 



N- 
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Lemma 3.6. Assume that the setting in Suhsection \3.1\ is fulfilled, that the setting in the beginning of Subsec- 



tion \S.S^ is fulfilled and that 
Then 



for all w G N. 



[0,T] 





- 




- 


lim E 

\f— i-oo 


sup 

ne{0,l,...,<5^} 




^ n 



'^ is [jj.^ a) -consistent with respect to Brownian motion. 
(160) 



Lemma 3.7. Assume that the setting in Subsection \3.1\ is fulfilled, that the setting in the beginning of Subsec- 
tion \3.°A is fulfilled and that cj): M'^ x [0,T] x R™ — >■ M.'^ is (iJ^, a) -consistent with respect to Brownian motion. 
Then 

lim limsupP[J^ < A^] = 0. (161) 



TV- 



Lemma 3.8. Assume that the setting in Subsection \3.1\ is fulfilled, that the setting in the beginning of Subsec- 
tion \3.!A is fulfilled and that (f): M'' x [0,T] x R™ — > R'' is [ji, a) -consistent with respect to Brownian motion. 
Then 



lim 



sup 

■ie{o,i,...,N} 






= 



(162) 



for all e G (0, oo). 



The proofs of Lemmas l3.5H3.8l are given below. The next corollary is an immediate consequence of Lemma[331 
and Lemma [131 



Corollary 3.9. Assume that the setting in Subsection \3.1\ is fulfilled, that the setting in the beginning of Sub- 
section COI is fulfilled and that (p: W^ x [0,T] x R™ — )■ W^ is {(x, a) -consistent with respect to Brownian motion. 
Then 



hm P 



sup ||Xr 

rie{0,l,...,JV} 



y^||>e 



= 



(163) 



for all e G (0, oo). 

Using Corollary 13.91 the proof of Theorem 13.31 is completed at the end of this subsection. We now present 



the proofs of Lemmas I3.5H3.8I Let us begin with the proof of Lemma [ 

Proof of Lemma \3.5\ Subadditivity and monotonicity of the probability measure P show 

P[rf < T] < Pht G [Q,T]:Xt i dJ\ +p[rf <T,\/te [Q,T]: Xt G D 



< 



< 



3te[0,T]:Xti D^\ + p[t„^ < T, X,n G D,,, Z^^ G {D,) 

3te[0,T]:XtiD^ +P sup \\Xt - Z^\\ > dist{D^,{D,Y 

te[o.T] 



(164) 



for all v,w,N G N with w < v. Lemma [3.41 and the estimate dist(_Dtu, (Z?„)'^) > for almost all i; G N hence 
give 



lim sup lim sup P[rf < T] < P 3t G [0,T]: Xt ^ D^ 

for all w G N. The continuity of the sample paths oi X : [0, T] x 51 — > 13 therefore yields 
lim sup lim sup PK^ < T] < lim P[3t G [0,r]: Xt i DJ = 

D^OO JV-i-OO UJ-J-OO 



and this completes the proof of Lemma 13.5 



(165) 

(166) 
D 



Proof of Lemma \3.6[ Throughout this proof, we use the mappings AW^ : ft — >■ R™, n G {0,1, ..., iV — 1}, 
N G N, defined by AW^ := VF („+i)t - Wnx for aU n G {0, 1, . . .iV - 1} and aU iV G N. This notation, in 
particular, ensures 



{kA5^)-l 



l=Q 

fe-1 

= Xo + J2 Ms->i} Uzn ■ i + -oiz!") AWi 



(167) 



/=o 
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and 



fe-i 



^feA5« =Xo + J2 Ms->i} ^{Yi"^ W^^^i 



(168) 



1=0 



for all k G {0, 1, . . . , N}, N G N and all w € N. The triangle inequality therefore implies 















1=0 


fc-i 

)^ !{,«>,} (a(Z,^ 


)-a(y,^))Al^,^ 




J. k-l 
1=0 


m") - f • E 


<^(y,^,^,AW^,^) F.. 




+ 


k-l 
1=0 


'{Yi^) l\Wi^ + E 


0(y 


-,^,Ai^-)^^; 


-4>{Y,\^^^wn) 



for aU /c e {0, 1, . . . , TV}, TV e N with j^ < U and all w e N. The definition of c„ e [0, oo), u € N, (see pM)) ) 
hence yields 



7N 



rN 






fce{0,l,...,n} 
rp ri-1 



TV 



/=0 



!a5« ~ ^/A<5"|| + ^^P 

fee{oa,...,n} 



fc-i 



Y,M5«>i}{o{Z^)-o{Yi''))^Wl' 



- T f sup ||m(x) - f • E[0(x, ^, W^T/Jv)] ll) 



+ sup 

0<fc<Af 



fc-1 

Y.l{S^yl}{HYi^)^Wi'' 
1=0 



(r-,^,Aw^-)|^. 



-0(y,^^,Aw^,^; 



P-a.s. for aU n e {0, 1, . . . , A^}, A^ e N with ^ < U and aU u e N. The Burkholder-Davis-Gundy inequality (see, 
e.g., Theorem 48 in Protter 62 ) therefore implies the existence of a real number k E [0, oo) such that 

2 



sup 

ke{0,l,...,n} 



yN xrN 

^kAS^ ~ ^kAS" 



4:T^\c^l ^ II ^ ^ 



N 



1=0 



Ll(a;R) 
2 
Ll(a;R<i) 



'' llLl(a;R'') 
l|2^ 



1=0 

4T2 ( sup ll^(x) - f • E[(/.(x, ^, W^T/Jv)] 

JV-l 

4^5: ||l{,«>,}(a(ll^)AP^,^+E[0(i^^,^,AT4^,^)|^^] -0(r,^,f ,AW^,^)) 



(=0 



2 

Li(n;l 



for all n e {0, 1, . . . , TV}, TV e N with j^ < ty and all v e N. The estimate 



E 



ks^>n{aiY,^)AWr+E[HYi''.^,^Wn\:F^]-m'',^.^Wn 



E 



1{5^>;}IE 



a(r,^) AM/,^ + E[0(y,^, ^, AW^,^) 1^^] - ^"^ , f , AM/,^; 



J'l 



< sup E 



r{x)WT/N+K4>{^, JJ,Wt/n)] - Hx, ^,Wt/n) 



for alU e {0, l,...,iV - 1}, iV e N with ^ < t^ and all w £ N and the inequality IJAW/ 



N\\2 



!l1(0;R'") 



< 
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|AWj^||^2(n.R™) = ^ foi" all / e {0, 1, . . . , A^ - 1} and all A^ e N hence show 



sup 

fcG{0,l, ■■-,"} 



Z, 



N 
fcAi5' 



-Y, 



N 

fcA(5/; 



< 



4T^ 

ATmn 

N 



N 



E\\yN Y 

^n— 1 m 



Ll(n;R) 
2 
Ll(n;R 



EEl|lW>0(-^(^."')~<^»(^"^))llL.(O;» 



, /=0 i=l 



+ 4T2 sup ||^(a;)-f .E[(/)(x,^,W^T/Jv)]| 



4kA^ ( sup 



r{x)WT/N + E[0(x, ^, W^T/Jv)] - (/-(a^, ^,Wt/n) 



Ll(0;E'') 



for all n e {0, 1, . . . ,iV}, A^ e N with ^ < <„ and all w e N. The definition of c„ G [0,oo), w G N, (see (fT54l) ) 
hence gives 



sup 



7JV 

^fcA(S« ^ fcA^iy 



y, 



TV 



Ll(0;lfi 



< 



/ 4Tm2|c.p(T + At; 



N 






.1=0 



Ll(a;R'') 



- 4T^( sup \\fi{x) - f • E[cj){x, ^,Wt/n)] \ 



ATkN 

-TfT- I sup 



Ll(0;R'') 



for all n e {0, 1, . . . , A^}, A^ e N with jr <U and all v G N. Moreover, note that (|150p implies 



lini —^- sup ||EU(a;,t,H/t)]|| =0 



for all u e N. This and P^ then yield 



lim — • sup 

t\,0 \ t 2;g_D„ 



cT{x)Wt + E[(/)(a;, i, W^t)] " 0(^, ^, W^t) 



Li(n;] 



for all u G N. Combining (|150p . (|169l) . (|171l) and Gronwall's lemma then shows 

2 



lim 

N-¥oo 



sup 

n6{0,l,...,A'} 



^nA(5« ^nA<5« 







L1(0;] 



for all u G N. This completes the proof of Lemma 

Proof of Lemma \3.7\ Subadditivity and monotonicity of the probability measure P show 

PK < iV] < P[t^ < T] + P[rJ - T, 5,f < A^] 

< p[tJ < T] +f[z^^ G D.^, y^^ ^ i?„ 

< p[t^ < T] +p[\\z^. - r,f II > dist{a^,{D,r) 

for all v^w, N Si N with w < t;. Markov's inequality and Lemma 13.61 therefore yield 

'E[ll^f«-y, 



ATI 



lim sup lim sup P [5^ < Af] < limsupP[T^ < T] + lim sup lim sup . . , .c 

w— >-oo N-^00 TV— >-oo u— i-oo N-¥oo \ dlSt (£)„,, (Z?^,) 

= limsupP[T^ <T] 

JV-i-oo 

for all w E N. Combining (|174p and Lemma 13.51 completes the proof of Lemma 13.71 



(169) 



(170) 
(171) 

(172) 
D 



(173) 



(174) 



D 
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Proof of Lemma\KE The identity {S^ < iV} l±) {S^ = N} = n and Markov's inequality imply 



sup \\Z, 

ne{0,l,...,N} 



N 



Y^W >e 



<P[S^ <N] 



s:N 



N, sup \\Z. 
rie{oa,...,5«} 



N 



Y^W >e 



< 



[5^ < iV] + - • E 



sup 

ie{Q,i,...,s^} 



\n'N vN\ 



(175) 



for all u, A^ e N and all e G (0, oo). Lemma [5771 and Lemma [5^ therefore yield 



sup 

•ie{o,i,---,JV} 



I 7^ V"^ll > P 



lim 

Ar-»-oc 



for all e G (0, oo). This completes the proof of Lemma 

Proof of Theorem\KM In order to show Theorem K^ let X^: [0,r] x n 
stochastic processes defined by 



(176) 



D 



X 



N 



n + 1 ] Xr 

T - 



tN 



n X(^+iy 



^''j A^ G N, be a sequence of 

(177) 



for all te[^, ^^^^^], n G {0, 1, . . . , iV - 1} and aU N eN. Then we obtain 



sup Xf - 


-Yf 


= sup 


X^T - 


_YN 


e[o,T] 




n£{0,l,...,N} 







(178) 



for all A^ G N. Moreover, the continuity of the sample paths of X: [0,r] x fi — > D yields 



lim sup \\Xt-Xn\ =0 



(179) 



-a.s.. This implies limAr^oo ]P[supjg[o yj \\Xt — X^\\ > e] = for all e G (0,oo) and CoroUarv 13.91 hence shows 



lim sup I 

TV— >-oo 



^Af| 



sup llXf-r/'ll >e 
te[o,T] 



< lim sup I 

iV-i-oo 



< lim sup I 

AT-J-oo 



sup \\Xt~X^\ 
te[o,T] 



sup \\X^-Yr\\>e 
te[o,T] 



(180) 



■^N\ 



sup \\Xt~Xi'\\ >§ 
te[o,T] 



lim sup I 

Af— i-oo 



sup 
ie{o,i.....,N} 



X„T - Y, 



N 



> 



= 



n 



with moment bounds for 



for all e G (0,00 ). This completes the proof of Theorem l3.3l 

3.4 Strong convergence 

In this subsection, we combine the convergence in probability result of Theorem 
the numerical approximation processes in (I148P to obtain strong convergence. 

3.4.1 Strong convergence based on moment bounds 

This subsection presents strong convergence results under the assumption that the numerical approximation 
processes in (1148^ satisfy suitable moment bounds. Below in Subsections 13.4.31 and I3.4.2[ we will give more 
concrete conditions on the drift coefficient /i, the diffusion coefficient a and the numerical method which allow 
us to avoid these moment bound assumptions. The strong convergence results in this subsection use the following 
well-known modification of Fatou's lemma. For completeness its proof is given below. 

Lemma 3.10 (A modified version of Fatou's lemma). Let (il, J^, P) be a probability space, let {E^cIe) he a 
separable metric space and let Z^: 57 — > _E, iV G N, and Z: D, ^ E be T /B{E) -measurable mappings with 
limAr_j.oo P[c^_e(^7V, •^) ^ &] = for all e G (0, oo). Then ^\ip{Z)\ < liminf7v->oo E[(/3(Z7v)] for all continuous 
functions ip: E ^>- [0^ oo]. 
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Proof of Lemma \3.1(A First, let (^: i? — > [0, oo] be an arbitrary continuous function. Then let N{k) G N, A; G N, 
be an increasing sequence of natural numbers such that 



liminfE[(^(Zw)] = lim E[(^(Z^(fc))] . 



AT— >-oo 



(181) 



Next note that limfc_>.oo P[d£;(Z7v(fe), Z) > e] = for all e £ (0, oo) by assumption. Consequently, there exists 



an increasing sequence fc; G N, / G N, of natural numbers such that lim, 



;_s.oo ^N{k,) 



Z P-a.s.. The continuity of 



if hence implies lim^_^.oo (p(Zjv(fe,)) = ^{Z) P-a.s.. Combining this, Fatou's lemma and (J18ip then gives 



E[(p(Z)] ==E 



lim (p(Zjv(fe,)) 



< liminfE[(^(Z7v(fc,))] = lim E[ip{ZN{k))] = liminf E[(p(Zjv) 



The proof of Lemma FS.lOl is thus completed. 



D 



Let us now present the promised strong convergence results which use the assumption of moment bounds for 
the numerical approximation processes in (I148p . 

Corollary 3.11 (Strong final value convergence based on moment bounds). Assume that the setting in Subsec- 
tion \3.1\ is fulfilled, suppose that (f>: M'' x [0,T] x M™ — )> R"^ is {^, a) -consistent with respect to Brownian motion 
and let p G (0, oo) be a real number such that 



lim sup E 

Af->-oo 



Y, 



N\\P 



< OO. 



ThenE[\\XT\\P] < oo and 



for all q G (0,p). 



lim E 

N->-oo 



Xj 



Y, 



N\\Q 



= 



(182) 



(183) 



Proof of Corollary \3AJ\ First, let q G (0,p) be arbitrary. Next observe that inequality (jl82p implies that there 
exists a natural number No £ N such that 



sup E 



Yn 



N\\P 



< OO. 



Theorem 13. 31 and Lemma 13.101 hence imply that 



EniXTirl <liminfl 



Yn 



N\\P 



< sup E 

Ne{No,No + l,...} 



Y 



N\\P 



< OO. 



(184) 



(185) 



This together with inequality (|184p shows that the family of random variables \\Xt — ^j^l!'', N G {A^o,^o + 
1,...}, is bounded in LP/'?(f2;R) and, therefore, uniformly integrable (see, e.g.. Corollary 6.21 in Klenke 41 ). 
Theorem 6.25 in Klenke ^41^ hence implies (|183p . This completes the proof of Corollary 13. Ill D 

Corollary 3.12 (Strong convergence based on moment bounds). Assume that the setting in Subsection \S.l\ is 
fulfilled, suppose that (p: M.'^ x [0,r] x R™ — > R'' is {^, a) -consistent with respect to Brownian motion and let 
p G (0, oo) be a real number such that 



lim sup sup E 

JV-i-oo ne{0.1,...,N} 






N IIP 



< OO. 



(186) 



Then sup^grp yi E[||Xt||Pl < oo and 



lim sup E 



\Xt - Y, 



Nni 



(187) 



for all q G (0,p). 

Proof of Corollary \3.12\ Inequality (jl86p implies 



lim sup sup E 

N-yoo te[0,T] 



Y 



N\\P 



< oo 



(188) 
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and this yields that there exists a natural number Nq £ N such that 



sup sup E 

Ne{NQ,No + l,...}te[0,T] 



Y, 



N\\P 



< CX). 



(189) 



Theorem 13.31 and Lemma 13.101 therefore show that 



sup I 
te[o,T] 



Xt 



Moreover, the identity 1 



< sup lim inf E 

te[Q,T] w-»-oo 

ll 



Y, 



N\\P 



< lim sup sup ] 

iV-i-oo te[o,T] 



Y 



N\\P 



< c». 



(190) 



sup E 

te[o,T] 



Xt~Yt 



{sup,^[o,^, ||x,-?^"||<i} + l{sup,^[o,T] !l^.-n"ll>i} == 1 ^'^^ Holder's inequality imply 

Nlll] ^ „„„ „ -, _ _ \\v vNUl 



< sup E 
te[o,T] 

+ sup E 
te[o,T] 



{sup,e[o.^, iix,-y«![<i} IF* ^t 



-a _ II V V' 

^{™P=6[0,Tll|Xa-n"l|>l} 11^* ^t 



N\\1 



<E 



, (p — g) , 



JV||9 



>Af| 



sup \\Xt~Yt''\\ >1 
te[o,T] 



sup sup . 



Xt-Y, 



M\\P 



for aU N e {No, iVo + 1, . . . } and aU q e (0,p). The estimate \a + 6|^ < 2p |a|'' + 2^ |6|^ for all a,beR and the 
inequalities (jl89p and (jl90p therefore give 



sup E 
te[o,T] 



P^t-^t 



N\\9 



<E 



min ( 1 , sup Xt — Yj 

te[o,T] 



Af||? 



(191) 



+ 2P 



sup llXt-Fj^ll >1 

t€[0,T] 



sup sup E 

, j\/e{Afo,Afo+i,...}te[o,T] 



Yt 



M\\P 



■V 



sup |ix,~i>^|i>i 

te[o,T] 



sup E 
,te[o,T] 



X, 



< oo 



for all A^ e {A^'o, A'o + 1, . . . } and all q G (0,p). Moreover, Theorem 13.31 implies 



lim P 

N->-oo 



yNi 



sup \\Xt^Yt"\\ >1 

*G[0,T] 



= 



and 



lim E 



min I 1 , sup Xt — y^ 

te[o,T] 



N\\1 







(192) 



(193) 



for all q e (0, oo). Combining ([1M|) - P^ then shows p^T]) . This completes the proof of CoroUarv 151^ D 

Corollary 3.13 (Uniform strong convergence based on a priori moment bounds). Assume that the setting in 



Subsection \3.1\ is fulfilled, suppose that 



[0,T] 



motion and let p G (0, oo) be a real number such that 



Then E[sup(grQ2ni HXtp] < oo and 



lim sup E 

N-+00 



hm E 

N—^oo 



sup 

ne{Q,l,...,N} 



\V^ IIP 



"^ is {^, a) -consistent with respect to Brownian 

(194) 



< oo. 



sup \\Xt — Yf 
te[o,T] 



NUI 



(195) 



for all q G (0,p). 

The proof of Corollarv l3.13l is analogous to the proof of Corollarv 13 . 1 1 1 and therefore omitted. Corollarv l3.13l 
extends Theorem 2.2 of Higham, Mao & Stuart [50] in several ways. First, the moment bound condition on 
the exact solution in Assumption 2.1 in |30j is omitted in Corollarv 13.131 This assumption can be omitted in 
Corollarv 13.131 since (J194I) and Lemma [3.101 implv E[sup(grQ'ri ||Xt||P] < oo. Moreover, Theorem 2.2 in [30] 
proves uniform strong mean square convergence for the Euler-Maruyama scheme while Corollary 13.131 proves 
uniform strong L'-convergence with q e (0,p) for (/n, CT)-consistent one-step schemes of the form (|148p . 
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3.4.2 Strong convergence based on semi stability 

Roughly speaking, the next corollary asserts that semi stability with respect to Brownian (see Definition | 
together with consistency with respect to Brownian motion (see Definition I3.ip and with the a priori growth 
bound (J196I) implies strong convergence of the numerical approximation processes (|148l) to the solution of the 
SDE (jl47p . Its proof is a direct consequence of Corollary 12.91 and of Corollary 13. 121 and is therefore omitted. 

Corollary 3.14. Assume that the setting in Subsection \3.1\ is fulfilled, let a,r Cz (l,oo], p G (0, oo), 6 G (0,r], 
let V: M"^ —7- [0,oo) be a Borel measurable function with supj.gjjd ]_!y/-^-v < oo and M[V{Xq)] < oo, assume that 
0: R'^ X [0,T] X M™ -> R'' is {fi, a) -consistent with respect to Brownian motion, assume that R'^ x [0,0] x R'" 3 
{x, t,y) t—^ X + (f){t, X, y) e R is a-semi V -stable with respect to Brownian motion and assume that 

limsup sup (iV('-")('/''-'/(^'-»||?,^/A,|lLp.-(0:M<')) <oo. (196) 

Af->-oo 7ie{0.l,...,N} ^ ' 

r/ienlimsup^^^suptg[o.T]E[||Xt||P+||yt^p] < <x anrf limjv^oo suptg[o,T] IE[||^t-i^t^||''] ^ for all q E (0,p). 

3.4.3 Strong convergence of an increment-tamed Euler-Maruyama scheme 

In Subsection 13.4.11 strong conyergence of numerical methods of the form (jl48p has been proyed under the 
assumption of suitable moment bounds for the numerical approximation processes (|148l) . The next result proyes 
strong conyergence of an increment-tamed Euler method and imposes appropriate assumptions on the coefficients 
gander of the SDE piT]) . 

Theorem 3.15 (Strong conyergence of an increment-tamed Euler-Maruyama scheme). Assume that the setting 
in Subsection \3.1\ is fulfilled, let p G [3, oo), c, 70,71 G [0, 00), let fl : R"^ ^ R'^ , a : R'^ ^ jjdxm ^g Borel measurable 
functions with /IId = /i and a\o — a, let V G Cp{R'^, [l,oo)) with E[y(Xo)] < 00 and with 

{G^,^V)ix)<c-Vix), mx)\\<c\Vix)\i^\, \\aix)\\Li^^^j,.)<c\Vix)\i^^\ 
for all xeR'^ and let F^ : [0, T] x fi ^ R'', N eN, satisfy 

"^ ^^ ' rna.^{l,WiYnT,N)^'^^iY^T/N)iW(n+l)T,N-W^T,N)\\) 

for all t e {nT/N, {n + 1)T/N], n € {0, 1, . . . , iV - 1} and all N eN. Then we obtain 

lim sup E[\\Xt-Y/^\\'^]=0 (198) 

for all q € (0,oo) with q < 2-yi+4mai(7o,7i,i/2) ~ 5 «'^'^ limsup^NqSup^gRd \\x\\''/Vix) < 00. 

Theorem 13.151 is a direct consequence of Corollary 12.211 Corollary 13.121 and Lemma 13.281 The next result is 
a special case of Theorem 13.151 

Corollary 3.16 (Potencies of the Lyapunov-type function). Assume that the setting in Subsection \3 . l\ is fulfilled, 

let p e [3,00), q £ [1,00), c, 7o,7i G [0, 00), let fi: W^ ^ R"*, a: R"^ ^ jjdxm ^^ Borel measurable functions 

with P\d = M and a\D = cr, let V £ C^(R'', [1, 00)) with E[y(Xo)] < 00 and with \\p{x)\\ < c\V{x)\^^^'\ , 
Mx)\\LiR'^M^) < c\V{x)\i^^ and {g-^^aV){x) + ^||T^'(a;)a(x)||2^5(j,„_R) < c-V{x) for all x £ R'^ and let 
Y^ : [0, T] X rj ^ R'', N £ N, satisfy (|TW1) . Then limAr^oo sup^^fo j,] E[||Xt - Y^^f] = for all r £ (0, 00) with 

r < 271+4 max('7o.7l.l/2) " ^ ""'^ hm SUp.x^, SUp.gRd ^^ < OO . 

Corollary 13 . 161 follows immediately Corollary l2.221 Corollary 13. 121 and Lemma r3.28l The next corollary giyes 
sufficient conditions for strong L'-conyergence of the increment-tamed Euler-Maruyama method (jl97p for all 
q£ (0,00). 

Corollary 3.17. Assume that the setting in Subsection \S.1\ is fulfilled, let c £ (0, 00), let fi: R'^ —>■ R'', 

ct: R'^ — > Rf*^™ be Borel measurable functions with fllu = fj, and a\D = cr, let V: W^' — >■ [l,oo) be a twice 
differentiable function with a locally Lipschitz continuous second derivative, with linisup„x osup^jg^d yT^ < 00, 
'"«^'i ELi ll^''''(2;)||L(')(Rd.R) < c|y(a;)|[i"i/^l for X^d-almost all x £ R'' , with E[y(Xo)] < 00 and with 



sup 

xeTSL'' 



{gf.,^V){x) , r\\V'{x)a{x)\\' 



Vi^) \Vix)\ 



2 



< 00, sup 

xeM'' 



\lJ.{x)\\ + ||o-(2;)||l(R" 



(i + ii:^r) 



< 00 
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for all r G [0,oo) and let the sequence of stochastic processes Y^ : [0,T] x il — )• M"*, N eN, satisfy p97p . Then 
we obtain that limjv^oo supjgrg 7-1 ]E[||-'^t — ^t^ll''] = for ali q G (0, 00). 

Corollary 13.171 follows directly from Corollary 13.161 The next corollary specializes Theorem 13.151 to the 
function 1 + \\x\\p, x € M'', with p e [3, 00) appropriate as Lyapunov-type function and is the counterpart to 
Corollary 12.241 It is an immediate consequence of Corollary 12.241 Corollary 13.121 and Lemma 13.281 



Corollary 3.18 (A special polynomial like Lyapunov-type function). Assume that the setting in Suhsection VJ . 1\ 
is fulfilled, let 0,70,71 S [0,oo), p e [3, 00), assume E[||Xo||p] < 00, let /2: W^ -^ W^ , a: M'^ -^ R'^^™ be Borel 
measurable functions with Ji\D — /i, a\D — c and 

{x,mx)} + (£^||a(x)||^^(jj„ R.) < c (1 + \\xf) , (199) 

|l/2(x)||<c(l + |lx|j[^''+^l) and ||a(a;)|L(K™,R.) < c (l + ||x|| [^] ) (200) 

for all X e M'' and let Y^: [0, T] xfi ^ M'', N eN, satisfy p^ . Then we obtain that\[mN^<x, supt^[Qj.]E[\\Xt- 
r/V||'?l = for all q G (0, 00) with q < ,^—-, 2 _ 1 

tllj •' 'i'-viy "i 27i+4inax(7o,7i,l/2) 2 

3.5 Weak convergence 

Convergence in probability implies stochastic weak convergence. The next corollary is thus an immediate con- 
sequence of Theorem [ 



Corollary 3.19 (Weak convergence with bounded test functions). Assume that the setting in Subsection \3.1\ 
is fulfilled, suppose that (p: W^ x [0,r] x R™ — )■ M'' is {^, a) -consistent with respect to Brownian motion and let 
[E, ll'll^) be a separable W-Banach space. Then 

Jirn^ \\E[fiX)] - E[/(F^)] 11^ = (201) 

for all bounded and continuous functions f: C([0,T],R'*) — > E. 

Corollary 13.191 follows directly from Theorem 13.31 since the test functions in Corollary 13.191 are assumed 
to be bounded. The case of unbounded test functions is more subtle and is analyzed in the sequel. The 
next lemma proves weak convergence restricted to events whose probabilities converge to one sufficiently fast 
with possibly unbounded test function under the assumption that the discrete-time stochastic processes Y^/j^ , 
n E {0, 1, . . . , N}, TV e N, are 0-semi F-bounded with 1^: R'' — > [0, oo) appropriate. 



Lemma 3.20 (Semi weak convergence). Assume that the setting in Subsection \3.1\ is fulfilled, suppose that 
(p: R"* X [0,T] X R™ — >■ R*^ is {ij,, a) -consistent with respect to Brownian motion, let {E, \\-\\^) be a separable R- 
Banach space, let V : R"^ — ?■ [0, cx)) be a Borel measurable function and assume that the sequence Y^,j^ : ft — > R'^, 
n G {0, 1, . . . , N}, N £ N, of discrete-time stochastic processes is 0-semi V -bounded with respect to ¥, i.e., let 
^N G <^n{Y^), N (z N, be a sequence of events such that 

limsup sup E[ln„T^(F„T/A,)] < oo and limsupP[(nAr)'=] = 0. (202) 

Af-i-oo 7ie{0,l,...,N} Af-foo 

Then E[V{Xt) + ||/(^t)|U] < oo and 

^lim^E[lo„||/(XT)-/(FT^)|U] =^linr^||E[lo„/(XT)] -E[lo„/(r/)] 11^ = (203) 

for all continuous functions f : M.'^ — > E which satisfy limsup^ j,]^ sup^g^d /|f y/ f\r < oo. 

Proof of Lemma [3.2(A Let /: R'' — > i? be a continuous function which satisfies limsup^/.^ sup^gjjd (i+y/i-f)^ *^ 

xeS.'^ (l + V{x)Y " 

\\f{x)\\^<7^{l + V{x)Y (204) 

for all a; S R'' . Next observe that Theorem 13.31 and assumption (|202l) show 



oo. The assumption limsup^y.]^ supj.gRd Jf y^^lfxr < oo ensures that there exists a real number r € (0, 1) such 
that rj :— sup^g^d .Hy,^|^N^^ < oo. We thus obtain that 



limsupP[||XT - lo„F^|| > e] < limsupP[||XT - 1o„^t|| > f ] + limsupP[ln„ ||Xt - Y^\\ > 

AT— >-oo N-^oo Af— >-oo 

< limsupP[(r2Ar)^] =0 

Af-i-oo 
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(205) 



for all e G (0,00). Lemma [3. 101 hence implies that 



E[V{Xt)] <liminfE[F(ln„r/)] < V^(0) + limsupE[lo„y(r2^)] < 00 

Af-J-oo N-^oo 



and Jensen's inequality and ()204[) therefore yield that 

E[\\f{XT)\\E] < r; • e[(1 + V{Xt)Y] < V ' {1+E[V(Xt)]Y < 00. 



(206) 



(207) 



In addition, observe that Theorem 13.31 implies limAr_j.oo IP[||/(-'^t) ^ /(i^T')!!^; > e] = for all e G (0, 00). 
Therefore, we obtain that 

>Wmi_ ^ .1 _r, (208) 



hm P[1oJ|/(Xt) - /(F/)|U > £] = 



for all e G (0,cxd). Next note that estimate (j204p ensures that 



limsup ln„/(y5i^) 



N- 



\L^/^(n;E) 



< rj ■ 1 + lim sup E 



N- 



,v{m 



< c». 



(209) 



Inequality (j209p . the fact that ^ > 1 and Corollary 6.21 of Klenke 41 show that the sequence l^jv ll/(i^)ll-Ej 
A'^ € N, of J-"/S(K)-measurable mappings is uniformly integrable. This together with (|207l) yields that the 
sequence Injvll/l-^r) -/(?5^)l!s, iV G N, is uniformly integrable. Equation (|208p hence shows (I203p . The proof 
of Lemma 13.201 is thus completed. D 



3.5.1 Convergence of Monte Carlo methods 

In Lemma 13.201 the convergence holds restricted to a sequence of events whose probabilities converge to one. 
In Proposition 13.221 below, we will get rid of the restriction to these events and prove convergence of the 
corresponding Monte Carlo method on an event of probability one. For proving this result, we first need a 
minor generalization of Lemma 2.1 in Kloeden & Neuenkirch 22] ■ More precisely. Lemma 2.1 in Kloeden & 
Neuenkirch [32] proves that convergence in the p-th mean with order /3 G (0,c») for all p G (0,oo) implies 
almost sure convergence with order /3 — e where e G (0,oo) is arbitrarily small. The next result is a minor 
generalization of this result and, in particular, proves for every arbitrary fixed p G (0,oo) that convergence in 
the p-th mean with order /3 G (- , 00) implies almost sure convergence with order /3 — 1/p — e for every arbitrarily 
small £ G (0,00). 

Lemma 3.21 (L^-convergence with order (3 G (1/p, 00) implies almost sure convergence). Let M G N, let 

{E,£, 11) be a measurable space and let Yn : E ^ R, N d {M, M + 1, . . .}, be a family of £/B{S.) -measurable 
mappings. Then 



sup AT" • \Yn\ 

Ni£{ALM+l,...} 



< 



LP(E-] 



E^'"" 



./3)p 



.N=M 



sup N^-\\Yn\ 

N<£{M,A4+1,...} 



LP{E;B 



(210) 



for all a,(3 gM. and all p G (0,oo). In particular, if s,uppffztM,M+i,...}{^'^ \\Yn\\lp{E;R)) < 00 for one p G (0, 00) 
and one /3 G M, then 



sup (iV" • IYatI) i d^i < 00 and sup (TV" • 1^^!) < 

E [Ne{ALM+l,...} I Ne{M,M+l,...} 



CX3 fi-a.s. 



(211) 



for all a G (— oo,/3— i). Moreover, J/supjvg{M,j\/+i,... jl^*^ II^a^IIlp(£;;R)) < 00 for all p G (0,00) andonefi^M., 
then 

Y 

sup (iV" • |yAr|) W^ < 00 and sup (iV" • 1^^!) < 00 fi-a.s. (212) 

E [Ne{M,M+l....} J Ne{M,M+l....} 

for all a G (~oo, /?) and all p G (0, 00) . 
Proof of Lemma [K2l[ Note that 

sup {N" ■ \Yn\)\ df, = [ I sup {N'^P ■ lYNf)] dfi 

E I Afe{Af,M+l,...} I JE I A'e{A/,Af+l,...} I 



< 



•'E y^^M ) N=M ^ •'^ 



(213) 



< 



E^^" 



-I3)P 



\N=M 



sup Nf^P ■ / \YN\^dfi 

NG{M,M+l,...} Je , 
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for all a,^ e M and all p e (0,oo). This proves inequality (|210p . The assertions in (|21ip and (|212p then follow 
immediately from inequality (|210p . The proof of Lemma 13.211 is thus completed. D 



Note also that the last assertion in Lemma [3.211 i.e., the assertion in (|212p . is essentially Lemma 2.1 in |32] 
generalized to arbitrary measure spaces. 

Proposition 3.22 (Convergence of Monte Carlo methods). Assume that the setting in Subsection \3.1\ is fulfilled, 



suppose that (p: 



[o,r] 



t^ is [n, a) -consistent with respect to Brownian motion, let n G N, a G 



[n + l,oo), let V: M"^ — )• [0,oo) be a Borel measurable function, assume that the sequence Y^/j^'- ^ -^ K'^, 
k e {0, 1, . . . , N}, N d N, of discrete-time stochastic processes is a-semi V -bounded and for every N d N let 
y '™ : [0, T] X rj — >■ R, m G N, be a sequence of independent stochastic processes with Py-w = Py-w,™ for all 
N, m e N. Then 



lim 

AT-j-oc 



E[/(Xt)] 



^N 



EZ=^fiY^n 



m- 



= 



l/MI 



which satisfy limsup^ ;itia2 sup^g^d /-^ ' y/^U 



< oo. 



P-a.s. for all continuous functions f : R** - 

Proof of Provosition VJ .2^ Let / : R'' — > R be a continuous function with limsup^yiiiA2 sup^-gRd (i+yi^))!- 
Next let M G N be a natural number and let A^q g *^''^'^ 



l/WI 



(214) 



< oo. 



sup I E 

Ne{M.,M+l,...} 



^{(n"/«)'=M0,i....«}6^«}^^"^^ 



T^) 



+ Ar".p 



i,...,N} ^ TV € N, be a sequence of sets such that 

(215) 



(^fcT/Ar)fce{0,l,...,A'} ^ {^n) 



< OO. 



Such a sequence of sets indeed exists since the sequence (?j^/jY)fce{o,i,...,JV}5 A^ G N, of discrete-time stochastic 
processes is assumed to be a-semi T^-bounded. We define now two families iljv G J^, N G N, and ilN,m € ^, 
A^, m e N, of events by 



riN ■■= |(^fcT/Ar)fce{0,l,...,Ar} G Ajvj 



and ^w,™ := |(yfeT/")fce{04,...:Af} ^ ^A'} (216) 

for all N, m e N. Note that ^[^1^] = ^[^N.rn] for all A^, m e N. Moreover, observe that Lemma ^?M implies 

hm E[1(,,„).|/(Xt)|] =0 and lim |E[1o„/(Xt)] - E[lo„/(r/)] | = 0. (217) 

In the next step note that the condition limsup^^iiA2 supj,g^d (iIyVXy < °o shows the existence of a real 
number r £ (0, ^^^) with r] := sup^gRd ,iJy,Sy < oo. The Burkholder-Davis-Gundy inequality (see, e.g.. 
Theorem 48 in Protter [62]) then shows the existence of a real number k G [0, oo) such that 

Elli {in.^^fiY^n - E[lo„/(r/)] } 



TV" 



TV" 



<i^ Ellin«.>/(^"'^ 



< 



2K\\ln^fiY^ 



N\ 



lLl/'-(n;B 



/ATr; 



< 



Ll/'-(n;R) 

-E[io„/(y-/)]||U(^, 

2K77||ln„|l + y(r/)r 



1/2 



Ll/'-(0;l 



(218) 



/TV" 



2^7?{E[ln^(l + V(y jy))]}'' ^ 2Krj{l+E[ln^ViY^)]y ^ 2Kr^{l + c) 



for all N e{M,M + 1,...}. Therefore, we obtain 



sup I N-^ 

Ne{M,M+l,...} 



TV"/2 



E[la,/(rjy)]} 



TV"/2 



N" 



< oo. 



(219) 



Ll/>-(0;] 



The condition ^ > r and the second inequality in (|21ip in Lemma [3.211 hence yield 



li^ / E:^i{in..,„/(y^^-)-E[io,/(y/)]} \ ^ 



JV-i-c 



TV" 



(220) 
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P-a.s.. Furthermore, observe that 

hm Ir „o ij-f<" ro^ -ici(a;) = lr„oo ,,oo u'f" fn^ )-=]i^) ^ ^U 
for all w e O. In addition, the condition a > n + 1 gives 



iinsup^_ 



.(u~:,(a„,„)=)](^) (221) 



N=M 

and the Borel-Cantelli Lemma thus shows 



oo oo oo 



< oo 



N=M 



N=M 



limsup (U,^li(iljv,m)' 

AT-i-OO 



Combining (I22ip and (|223l) with the estimate 



Y.H^Mn.^^YM{y^n\ A^--^\^^^^ 



< 



EllAfiyTni 



< 



us==«u^"i(njf,„)- 



T 



= 0. 



(max,„g{i_2^...,Ar"} l(o„,„)<=) 



N" 



for all A^ £ N results in 



TV" 



lim 



^N 



EZ^,Mn..^r\fiYTn\ 



-a.s.. In the next step observe that the triangle inequality implies 



= 



E[/(^t)] 



^N 



E:=i/(^t''") 



AT" 



< 



< 



|E[/(Xt)] -E[1o„/(Xt)]| + |E[1o„/(Xt)] -E[1s,„/(F/)] 



^Af 



E:=iln„,„/(yT'") 



TV" 



E;U(in„,„-i)/(r, 



Af,mN 



TV^ 



+ 



E[io„/(r/) 



|E[ln„/(XT)]-E[lo„/(r/)] 



TV" 



TV" 



(222) 



(223) 



(224) 



(225) 



(226) 



for all TV e N. Combining (pT7| . (j^ . ((^^ and (^^ finally implies (EH]). The proof of Proposition [S22] is 
thus completed. D 



3.5.2 Convergence of the Monte Carlo Euler method 

Next we prove almost sure convergence of the Monte Carlo Euler method. This result generalizes Theorem 2.1 
in |33| which assumes fj, to be globally one-sided Lipschitz continuous, a to be globally Lipschitz continuous 
and /i and a to grow at most polynomially fast. Corollary 13.231 is a direct consequence of Theorem 12.131 of 
CoroUarv 12.61 and of Proposition 13. 221 and its proof is therefore omitted. 

Corollary 3.23 (Convergence of the Monte Carlo Euler method). Assume that the setting in Subsection \3.1\ 

is fulfilled, let p G [3, oo), c, 70,71 G [0,oo) be real numbers with 71 + 2(70 V 71) < p/4, let fi: M.'^ — )■ R'^, 
be Borel measurable functions with Pl\d — /i, (t|d — <j and (j){x,t,y) — fL(x)t + &{x)y for all 



a: 



pdx? 



G K'', i e [0,r], y G M™, let V G Cl{W^, [l,oo)) with ¥\V{Xq)\ < 00 and with 



20+11 



{Gf.,aV)ix)<C-V{x), \\fl{x)\\ <c\Vix)\l — 



\a{x)\\ 



,M^)<c\V{x)\i^F- 



for all a; G M'^ and for every N ^ N let Y^'"^ : [0, T] x 51 — S- R, tti G N, be a sequence of independent stochastic 
processes with Pj-w = Py-w.m for all TV, m G N. Then 



lim 

Af— i-oo 



E[/(^t)] 



^N 



j:z=^fiY^n 



TV2 







(227) 



-a.s. for all continuous functions f: 



1/(^)1 



which satisfy limsup^^i]^ sup^-gRd (i+vMY ^ ^^' 



More results for approximating statistical quantities of solutions of SDEs with non-globally Lipschitz contin- 
uous coefficients can, e.g., be found in Yan [77] and Milstein & Tretjakov [55] . 
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n e {0, 1, . . . , iV - 1}, iV e N, by AW^ := W^n+i)T/N - WnT/N for all n G {0, 1, . . . iV - 1} and all iV e N. 



3.6 Numerical schemes for SDEs 

The purpose of this subsection is to present a few examples of numerical schemes which are (/i, cr)-consistent with 
respect to Brownian motion. Theorem 13.31 then shows that the approximation processes (jl48p corresponding to 
such schemes converge in probability to the exact solution of the SDE (|147l) . We now go into detail and describe 
the setting that we use in this subsection. Throughout this subsection, assume that the setting in Subsection 13. II 
is fulfilled and let il:W^ — >■ K'' and tr: R*^ — > W"''^''- be two arbitrary Borel measurable functions satisfying 
ii{x) — ^(x) and d-{x) — (t{x) for all x E D. Moreover, define mappings Y^^ : J7 -^ R'', n S {0, 1, . . . , A^}, 
A^ e N, by Y^ := 1;^/^^ for all n e {0,1,..., A^} and aU A^ e N and define mappings AW^ : n ^ R", 

A^ e 

Using this notation, we get from (|148l) that 

for all n E {0, 1, . . . , A^ — 1} and all N £ N. The following subsections provide examples of Borel measurable 
functions 0: M'' x [0,T] x R™ — >■ R'^ and numerical approximation schemes of the form (|228p . respectively, which 
are (/i, CT)-consistent with respect to Brownian motion. 

3.6.1 Some simple numerical schemes for SDEs 

Let 770: R'' X [0,T] X R" ^ R'*, 771 : R'^ x [0,T] x R™ ^ R'^^'^ and 772 : R"* x [0,r] x R"* ^ R'^'^'^ be Borel 
measurable functions. The next lemma then gives sufficient conditions to ensure that schemes of the form 

Y^l, = Y,^ + ^o{Y,:',^,AW,^)+m{Y,:',^,AW,^)KY:')^+V,{Y:',^,AW^)aiY:')AW:' (229) 

for all n e {0, 1, . . . , A^ — 1} and all A^ e N are (/i, (T)-consistent with respect to Brownian motion. 

Lemma 3.24. Assume that the setting in SuhsectiorMlMs fulfilled, let p.: W^ -^W^,a: W^ ^ R™^'*, 771,772: R'^x 
[0,r] X R" -> R''^'* and 770 : R'^ x [0,T] x R™ -^ R'' be Borel measurable functions with J1\d = M and a\D = a 
and assume that 

limsup ( supE[||772(a;,i,VFt)cr(a;)W^t||] ) < 00, (230) 

t\o \xeK J 

limsup ( sup e[ ||77i(a;, i, Wt) ~ IW^^^.^ + Wmi^, t, Wt) - I\\l^^,^ 1 ) = 0, (231) 

limsupf ^- supE[||77o(x,t,T^t)|l] ) = limsup ( - • sup ||E[77o(x, t, W^t)] 11 ) =0, (232) 

limsupf - • sup ||E[772(x,t,W^t)CT(a;)VFf]|| ) =0 (233) 

t\o \t xeK J 

for all compact sets K d D. Then 0: R'^ x [0,T] x R™ — > R'* given by (f>{x,t,y) — rio{x,t,y) + rji{x,t,y)fl{x)t + 
ri2{x,t,y)(j{x)y for all x G M.'^ , t G [0,r], y G R™ is {^, a) -consistent with respect to Brownian motion. 

Proof of Lemma \3.24\ The triangle inequality and the Holder inequality give 

-^ ■ sup ¥.[\\a{x)Wt - (f){x,t,Wt)\\\ 
\Jt xeK 

< 4= • supE[||77o(x,i,P^,)ll] + 4= • snpE[\\r„ix,t,Wt)K=c)t\\] 

Vt x^K \t xeK 

+ \- su^¥.[\\{I ~ 7^2{xXWt))a{x)Wt\\] (234) 

vr xeK ^ ' 

<i^.supE[||77o(a-,i,W^OII]+Vi(supE[||7,i(a;,i,Wt)|U(R.)] ) (sn^Mx)\\ 
\Jt xeK \xeK ' ^ -" 



xeK 



\\Wt\\ mn;R^) 



—n-^ ( sup 11772(2;, i,W"f) -/||L2(j^.i(Rd)) I ( sup ||cr(a;)||L(R'",R<i) I 

vr \xeK / \xeK / 
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for all t e (0,T] and all compact sets K C D. Combining (gM]), (P?T|) and ([^^ then shows pi^)) . It thus 
remains to establish (|150l) in order to complete the proof of Lemma 13.241 To this end note that 

limsup 1 sup ||/i(a;) — j • E[(/)(x, i, Wt)] || I < limsup 1 j ■ sup ||E[7]o(a;,i, Wt)] 
t\o \xeK J t\o \ xeK 

+ (limsup sup E[||?7i(a;,i,Wt)-/|U(Kd)] ) ( sup ||^(x)|| ) (235) 
\ t\o xeK I \xeK 



+ limsup J- sup ||E[?72(a;,i, VFt)cr(a;)pyt]| 
t\,o \ xeK 

for all compact sets K C D. Inequality ((2551) and equations dMH), (125^ and (1255)) then show p^ . This 
completes the proof of Lemma 13.241 D 

In the remainder of this subsection, a few numerical schemes of the form (|229p are presented. 

The Euler-Maruyama scheme In the case i]o{x,t,y) — 0, rii{x,t,y) — I and ri2{x,t,y) — I for all x E K**, 

t e [0,r] and all y e M™, the numerical scheme (I229P is the well-known Euler-Maruyama scheme 

Y.^+, = Y^ + fJ.{Y^)^ + a(r„^)AI^„^ (236) 

for all n S {0, 1, . . . , A^ — 1} and all A^ g N (see Maruyama [551). Of course, this choice satisfies the assumptions 
of Lemma 13.241 Combining Lemma 13.241 and Theorem 13.31 thus shows that the Euler-Maruyama scheme (|236p 
converges in probability to the exact solution of the SDE (I147p . In the literature convergence in probability 
and also pathwise convergence of the Euler-Maruyama scheme has already been proved even in a more general 
setting than the setting considered here; see, e.g., Krylov (45| and Gyongy [21]. Strong convergence of the Euler- 
Maruyama scheme, however, often fails to hold if the coefficients fj, and a of the SDE (|147p grow more than 
linearly (see !36 ) and therefore, we are interested in appropriately modified Euler-Maruyama schemes which are 
truncated or tamed in a suitable way and which therefore do converge strongly even for SDEs with superlinearly 
growing coefficients. Although strong convergence fails to hold, the corresponding Monte Carlo Euler method 
does converge with probability one for a large class of SDEs with possibly superlinearly growing coefficients; see 
Corollary 13.231 above. 

A drift-truncated Euler scheme In the case r]o{x,t,y) = 0, rii{x,t,y) — ^^xfi t]\i(x)\\) -^ ^^'^ ^2(2;, i,j/) = / 
for all X eR'^,t £ [0,T] and all y e M™, the numerical scheme (|229l) reads as 



max(l,^|i/2(y„A')||; 



r„%l = r„^ + ,rV"„-L;v^n^ + ^iYn)^W^ (237) 



for all n e {0, 1, . . . , A'^ — 1} and all N e N. In the case where the drift p. in (|237p is of gradient type and the 
noise is additive, i.e., a'{x) = <t(0) and p.{x) = —{'VU){x) for all x G Mf^ and some appropriately smooth function 
U : W^ — )■ M, the scheme (|237p has been used as proposal for the Metropolis- Adjusted-Langevin- Algorithm 
(MATLA; see Roberts & Tweedie [63]). The choice in (j237p also satisfies the assumptions of Lemma 13.241 
Combining Lemma [3.241 and Theorem 13.31 hence proves that the drift-truncated Euler scheme (|237p converges 
in probability to the exact solution of the SDE (|147p . If the diffusion coefficient a in (|237p grows at most 
linearly, then moment bounds and strong convergence of the drift-truncated Euler scheme (j237p can be studied 
by combining Theorem 12.131 Lemma 12.181 Corollary 12.91 Lemma 13.241 and Corollary 13.141 More precisely, 
Thcorcin l2.13l and Lemma 12.181 can be used to prove, under suitable assumptions on p. and a (see Theorem l2.13l 
for more details), that the drift-truncated Euler scheme (|237p is a-semi F-stable with respect to Brownian motion 
with a G (0,00) and T^: K'^ — >■ [0,oo) appropriate. Combining this and Corollarv 13 . 141 with the fact that (I237P is 
(/x, (T)-consistent with respect to Brownian motion according to Lemma [3 . 241 finallv proves, under the additional 
assumption that a grows at most linearly, strong convergence of the drift-truncated Euler scheme (1237p . 

A drift-tamed Euler scheme A slightly different variant of the drift-truncated Euler scheme (j237p is a drift- 
tamed Euler-type method considered in {34]. More precisely, in the case r]o{x,t,y) — 0, rii{x,t,y) — Y+tlP7vnT^ 
and r]2{x, t,y) = I for all x e M'', t e [0, T] and aU y £ R™, the numerical scheme (j229p reads as 



^"""+1 = ^n + , t:-Zn.u + ^(y:')aw:^ (233) 
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for all n e {0,1,..., N— 1} and all N E N. li D = K"*, if a is globally Lipschitz continuous and if fi is continuously 
differentiable and globally onc-sidcd Lipschitz continuous with an at most polynomially growing derivative, then 
strong convergence of the drift-tamed Euler scheme (|238p with the standard rate ^ has been proved in ^34j . 
Moment bounds and strong convergence of the drift-tamed Euler scheme (j238p in a more general setting can be 
obtained in precisely the same way as illustated for the drift-truncated Euler scheme (j237p . 

The Milstein scheme hi addition to the setting described in the beginning of Subsection l3.6l assume in this 
paragraph that a — {(Jij)if:^i d},je{i,....m} — {'^j)je{i m} ■ K'^ — > R''^™ is continuously differentiable. In the 
case rji{x,t,y) — r]2{x,t,y) — I and 

^ d m , d m 

Voix, t,y) = ^J2J2 {■£:^^)i^) ■ ^k,j{x) ■ y^ ■ Vj - -^^ {-£^a^)ix) ■ ak^x) (239) 

k—l i,j — l k—1 i—1 

for all X — {xi, . . . , Xd) £ M.'^, t £ [0, T] and all y = {yi, . . . , j/™) e M™, the numerical scheme (|229p reads as 



y,r+i = y:^' + hy:" ) i + ^(Yn ) aw^^ 



N.j 



^ ■ - ■ ■ - (240) 



k—l ij—l 
d rn 

EE(afr-0(i^.r)-M(>^„") 



rp d 



2N 

k=l 1=1 



for all n e {0, 1, . . . , A^ - 1} and aU iV e N where {AW^-'\ ..., AW^-"") = AW^ for all n e {0, 1, . . . , A^ - 1} 
and all N £ N. This choice satisfies the assumptions of Lemma 13.241 if and only if the commutativity condition 
(see, e.g., (3.13) in Section 10.3 in Kloeden & Platen [15] ) 



d d 



air*0(^) • '^k,j{x) = > {^^j){x) ■ akAx) (241) 



k=l k=l 



for all X G M"* and all i,j £ {1, 2, ... , m} is fulfilled and in that case, (I240p is nothing else but the well-known 
Milstein scheme (see Milstein [55] or, e.g., (3.16) in Section 10.3 in [43]). Combining Lemma l3.24l and Theorem l3.3l 
thus shows that the scheme (j240p converges in probability to the exact solution of the SDE (jl47p provided that 
(|24ip is fulfilled. In the literature, almost sure convergence with rate 1 — £ for £ G (0, 1) arbitrarily small and 
thus also convergence in probability of the Milstein scheme has already been proved in [31]. However, as in the 
case of the Euler-Maruyama scheme (see (|236p '). strong convergence of the Milstein scheme often fails to hold if 
at least one of the coefficients /i and a of the SDE (|147p grows more than linearly (see [35]). Nonetheless, if the 



drift term is tamed appropriately as in (|238p , then strong convergence of the corresponding drift-tamed Milstein 
scheme has been established in Gang & Wang [TS] for a class of SDEs with possibly superlinearly growing drift 
coefficients. 

Balanced impHcit methods Milstein, Platen and Schurz [ST] introduced the following class of balanced 
implicit methods. Let cq, ci, . . . , c™ : R'' ^- M.'^^'^ be Borel measurable functions such that the matrix 

m 

/ + co(a;)i + ^c,(x)|yj|eR'^^'^ (242) 

i=i 

is invertible for all x £ R'', t £ [0,T] and all y = (j/i,...,?/™) £ R™. This condition is, e.g., satisfied if the 
matrices cq{x, t,y), . . . , Cmic, t, y), {x, t, y) e R'^ x [0, T] x R™, are positive semi-definite. The associated balanced 
implicit method is then given through 

VN 

= y:' +{i+ co(r„^)i + e;1i '^.(^n"') I Aw^^^-|)"'(A(i;f')^ + hy:')aw:') (243) 

where {AW^'\ ..., AW^'"") = AW^ for ah n € {0, 1, . . . , A^- 1} and all N £N. For the case of at most linearly 
growing coefficients /x and a and uniformly bounded matrices cq, ci, . . . , Cm, Theorem 4.1 of Schurz [55] implies 
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uniformly bounded moments of the balanced implicit method (j243p . Moreover, under additional assumptions 
such as global Lipschitz continuity of /2 and (t, Theorem 5.1 of Schurz [BB implies strong mean square convergence 
of the balanced implicit method with convergence order ^ . In order to apply the above theory, write (j243p in 
the form (j229p with rio{x,t,y) = and rii{x,t,y) — 772(2;, i,y) = (/ + co{x)t + J2T=i ^ji^)\yj\)~^ ^'^^' *^1 ^ ^ ^"^j 
t G [0, T], y ~ (2/1, . . . , ym) G K™. Lemma l3.24l can be applied to derive conditions on the functions cq, ci, . . . , Cm 
which imply (/i, (T)-consistency with respect to Brownian motion. Theorem 13.31 then shows convergence in 
probability of the balanced implicit method (|243l) . Moreover, similar as described above for the drift-truncated 
Euler scheme, moment bounds and strong convergence of the balanced implicit method (j243p can be studied by 
combining Theorem 12. 131 Lemma [2. 181 CoroUarv 12.91 Lemma [3.241 and Corollarv l3.14l 

3.6.2 Comparison results for numerical schemes for SDEs 

The next lemma considers two numerical approximation schemes of the form (|148p and shows that if one of the 
two schemes is (/x, cr)-consistent with respect to Brownian motion and if the two schemes are close to each other 
in an appropriate sense (see (I244p '). then the other scheme is also (/Lt, (7)-consistent with respect to Brownian 
motion. Its proof is straightforward and hence omitted. 

Lemma 3.25 (A comparison result for consistency). Assume that the setting in Subsection \3.1\ is fulfilled and 
let (j>: M.'^ X [0, T] x M™ — >■ W^ be a function which is (/i, a) -consistent with respect to Brownian motion and which 
satisfies 

^^" '{x,t,Wt)-Hx,t,Wt)\\\)^0 (244) 



lim sup — • sup E 
t\o \t xeK 

for all compact sets K C D. Then (p: M'* x [0, T] x R™ — > M'' is (/i, a) -consistent with respect to Brownian motion 
too. 

The next lemma, in particular, illustrates that convex combinations of schemes that are (/^, cr)-consistent 
with respect to Brownian motion are (/x, <T)-consistent with respect to Brownian motion too. More precisely, the 
next lemma shows that schemes of the form 

Y^+, = r„^ + m • '^i (i^„^, ^, AI^„^) + r,, ■ MY.^, ^,AW^) (245) 

for all n £ {0, 1, . . . , N} and all A^ e N are (/i, a)-consistent with respect to Brownian motion provided that 
?7i,?72 € K are real numbers with r/i + ri2 = I and provided that 01,02: R'^ x [0,T] x R™ are (/i, cr)-consistent 
with respect to Brownian motion. Its proof is clear and therefore omitted. 

Lemma 3.26 (Generalized convex combinations of numerical schemes for SDEs). Assume that the setting 
in Subsection \3.1\ is fulfilled, let (f>i,(f>2'- R'' x [0,T] x R™ — >■ R'' be {fi, a) -consistent functions with respect to 
Brownian motion and let ?7i,?72 G K be two real numbers with ?7i +772 = 1- Then 4>: R"^ x [0,r] x R™ — > R'' given 
by (l){x,t,y) = 771 • (piix, t,y) + 772 ■ 4'2{x,t,y) for all x G R'', t £ [0,r], y G R™ is {^, a) -consistent with respect 
to Brownian motion. 

Finally, the next lemma gives a simple characterization of (/i, cr)-consistency with respect to Brownian motion. 
Its proof is straightforward and hence omitted. 

Lemma 3.27 (A characterization of consistency). Let T G (0,oo), d,m £ N, let D C M."^ be an open set and 
let fi: D — > R'' and a: D ^ jjdxm ^g locally Lipschitz continuous functions. A Borel measurable function 
(p: R'^ X [0, T] X R™ — >■ R'^ is then {n, a)-consistent with respect to Brownian motion if and only if there exists a 
Borel measurable function (p: R"^ x [0,T] x R™ — > R'' which is [fj,, a) -consistent with respect to Brownian motion 
and which satisfies 



and 



\im[^^-supE \\(l){x,t,Wt)-(f>ix,t,Wt)\\ )=0 (246) 



\im [ \ ■ sup \\E[cPix,t,Wt)] -E[^{x,t,Wt)]\\] =0 (247) 

'NO \ xeK / 

for all compact sets K d D where W : [0, T] x J7 — > R™ is an arbitrary standard Brownian motion on a probability 
space {rt,T,F). 
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3.6.3 Taming principles for numerical schemes for SDEs 

Let 0: R'' X [0,r] x M.™ ^ M'^ be a function which is (/i, cr)-consistent with respect to Brownian motion and 
consider a sequence Z^ : {0, 1, . . . N} x J7 — !• R'^, N G N, oi stochastic processes given by Zq = Xq and 



^n+1 — ^n 



7N T 



AW^) 



(248) 



for all n S {0, 1, . . . , A^ — 1} and all N E N. The next lemma then gives sufficient conditions to ensure that 
numerical approximation schemes of the form 



Y^'i, = r„^ + r?(y„^, ^, Aiy„^) 0(y„^, ^, aw^„^) 



(249) 



for all n e {0,1,..., A^ — 1} and all A^ e N are (/x, CT)-consistent with respect to Brownian motion where 



V- 



[0,T] 



l>dxd 



is a suitable Borel measurable function. 



Lemma 3.28 (Increment taming principle). Assume that the setting in Subsection \3.1\ is fulfilled, let rj: 



[0, T] X R" 



pdxd 



be a Borel measurable function, let (f>: R x [0,r] x 



(fj,, a) -consistent with respect to Brownian motion and assume that 

1 



lim sup 



sup E 



{x,t,Wt) 



< 00 



be a function which is 



(250) 



lim sup I — • sup I 
t xeK 



\r]{x,t,Wt) " l\ 



L(E<') 



t\0 

for all compact sets K C D. Then 0: R"* x [0,T] 

X e R , t S [0,T], y E R™ is {fj,, a) -consistent with respect to Brownian motion. 

Proof of Lemma \3.28[ Holder's inequality, equation (j250p and equation (|251l) imply 



(251) 



l'^ given by (j){x,t,y) = r]{x,t,y) 4){x,t,y) for all 



lim sup — • sup E 



< lim sup I — • sup E 

t x£K 



t\,0 



< lim sup 



(x,t,Wt)-(j>{x,t,Wt)\ 
\\ri{x,t,Wt) - IWuK-^) ■ 



sup^eKlE[||r/(a;,i,VFf)-/|| 



t 



\^{x,t,Wt)\\ 

vl/2 



(252) 



lim sup - • sup E 

, t\fl t xGK 



{x,t,Wt)\\' 



= 



for all compact sets K C D. Inequality (j252p and Lemma [3.251 then complete the proof of Lemma [3.281 

Let us illustrate Lemma [3.28l bv an example. More precisely, in the special case 77(2;, t, y) 
for ana;€R'', tE [0,T], y El 



D 



the scheme (I249P reads as 



^{l,t\\4,{x,t,y)\\) 



V^ 
^n+1 






max(l,^||0(r,f^,f,Aiy,f)||) 



(253) 



for all n G {0, 1, . . . , iV - 1} and all N eN. 

If we now additionally assume that the Z^ , A^ e N, in (12481) are Euler-Maruyama approximations, i.e., that 
(^(x, t, y) = fi{x)t + a{x)y for all x eW^, t E [0,T], y E R™, then ((^551) reads as 






fi{Y^)^ + -a{Yj:^)AW,{ 



max(l,^||/2(>;^)^ +a(r/)AI^,f II) 



(254) 



for all n E {0, 1, . . . , Af — 1} and all A^ e N. This increment-tamed Euler-Maruyama scheme clearly satisfies the 
assumptions of Lemma 13.281 Theorem 13 . 31 hence shows that the scheme (|254p convergences in probability to the 
exact solution of the SDE (jl47p . Note that this scheme is frequently studied in this article. Strong convergence 
of the scheme (|254l) is studied in Subsection 13.4.31 above (see also Section |4] for a list examples of SDEs in case 
of which the scheme (|254p has been shown to converge strongly) . 
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In the case of spatially discretized semilinear stochastic partial differential equations, another choice for 
the increment function (/>: R'' x [0,T] x M™ -^ M.'^ in (|248p naturally arises. More precisely, suppose that 
fl{x) = Ax+F{x) for all x S K'^ where A e R'^^'^ is a dx d-matrix with det{I-tA) ^ for all t g [0, oo) and where 
i^ : R'* — > M'' is a Borel measurable function and suppose that ^(a;, i, y) ^ {I — t.A)^ {x + F{x)t + a{x)y) — x for 
all X G M.'^, t £ [0,T] and all y S M"*. The approximations processes Z^, A^ e N, in equation (|248p thus reduce 
to the linear implicit Euler approximations 



Z„^+, = (/ - ^Ay' {Zf, + F(Z,f )^ + a(Z„^)AM^,f ) 



zN 



{I - I, A)-' {Z^ + F{Z^)I, + a(Z„^)AM/„^) 
for all n e {0, 1, . . . , iV — 1} and all A^ e N and the scheme in p53p then reads as 

(/ - ^Ay' [Y^ + F{Y^)^ + B{Y„^)AW,f) 



Z. 



N 



(255) 



-'ri+l 



Y 



N 



Y 



N 



max(l,^||(/-^^)-'(F„^ + F(r/)^+i?(y/)AW^„^) 



yNl 



(256) 



for ah n e {0, 1, . . . , iV - 1} and ah N eN. 

In dUHl), dinSl), (PMl) and dUS]), respectively, the increment function (^: M'' x [0,T] x M™ ^ M"* is tamed in 
a suitable way so that the scheme does not diverge strongly (see [3S]) and moment bounds can be obtained (see 
Subsections 12.1.31 12. 1.41 and 12. 2.3p . Instead of the increment function, also the whole scheme can be tamed in an 
appropriate way. This is subject of the next lemma. More precisely, the next lemma gives sufficient conditions 
to ensure that schemes of the form 



Yj:i,^v{Y:',^,Aw:'){Y,^ 



{Y,^,^,AWr^ 



(257) 



for all n £ {0,1,..., A^ — 1} and all A^ e N are (/i, CT)-consistent with respect to Brownian motion where 



[0,T] 



77: M'' X [0, T] X R"^ — >• R'^^'' is a Borel measurable function and where 0: 1 
which is (/i, cr)-consistent with respect to Brownian motion. 

Lemma 3.29 (Full taming principle). Assume that the setting in Subsection \3.1\ is fulfilled, let rj: 



fdxd 



be a Borel measurable function, let 1 



with respect to Brownian motion and assume that 



and 



[0,T]: 



limsup I sup E ||0(a;,t, IVt)|| 
t\,o \xeK 



l'^ is a function 

X [0, T] X 
W'' be a function which is {fj,, a) -consistent 

< 00 (258) 



lim sup — • sup E 



Hx,t,Wt) 



ni> 







(259) 



for all compact sets K <Z D. Then (, 
for all X £R'^,te [0, T] and all y £ 



l"^ X [0,T] X M'" -> R'^ given by (j){x,t,y) = ri{x,t,y){x + (j){x,t,y)) - x 
is (fi, a) -consistent with respect to Brownian motion. 

Proof of Lemma V3.2S\ Holder's inequality, the estimate (a + b)"^ < 2a^ + 26-^ for all a,b £R and equations (|258p 
and (P^ imply 



lim sup I - • sup E 

t\,0 \ i x£K 



< lim sup — • sup E 

t\a \t xeK 



< 2 lim sup — • sup ] 

t xeK 



t\o 



sup ||a;|| 

yX£K 



[x,t,Wt) 
\r]{x,t,Wt)-^ 

\v{x,t,Wt) 



(x,t,Wt) 



{x,t,Wt) 



lim sup sup E 



{xA,Wt) 




(260) 







t\0 x£K 

for all compact sets K C D. Inequality (|260|) and Lemma [3.251 then complete the proof of Lemma [3.291 D 

As an example of Lemma [3.291 consider the special case rj{x,t,y) — l[o,i] (i^||a; + (j){x,t,y)\[) ■ I^d for all 



x£R'^,t£ [0,r], y£ 



and the scheme (I257p reads as 



Y^ — 11 

-^n+l ^!yn+4,{Y«,^,AW'^)<^] 



(y„^ + 0(y,f,i,Aw/„^)) 



(261) 
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for all n E {0, 1, . . . , A^ — 1} and all N G N. A similiar class of approximatfons has been proposed in Milstein 
& Tretjakov [55] in which the truncation barrier -^ in the indicator set in (I26ip is replaced by a possibly large 
real number R G (0, c») which does not depend on N £ N. 

3.6.4 Linear implicit numerical schemes for SDEs 

The next lemma gives sufBcient conditions that ensure that linear implicit numerical schemes of the form 

_i (262) 



for all n G {0,1,..., A^ — 1} and all TV G N are (/x, cr)-consistent with respect to Brownian motion where 
A: R'' X [0,T] X M™ ^ R''^'' and &: R'^ x [0,r] x R™ ^ M"* are suitable Borel measurable functions (see 
Lemma [3.301 for the detailed assumptions). 

Lemma 3.30 (Linear implicit numerical schemes for SDEs driven by Brownian motions). Assume that the setting 
in Subsection \3.1\ is fulfilled, let Dy CC D, v G N, be a nondecreasing sequence of relatively compact open sets 
with Uye^Dy = D, let A-.R'^x [0, T] x R™ -^ R'^^^'^ and 6: R'' x [0, T] x R" ^ R'' be Borel measurable functions 
and let t^ G (0, T], u e N, be a sequence of real numbers such that /— yl(a;, t, y) £ R is invertible for all x G Dy, 
t e [0, ty], y e R™ and all v G N. Moreover, assume that R'* x [0, T] x R™ 9 {x, t, y) ^ A{x, t, y) x + b{x, t, y) G R"^ 
is {n, a)-consistent with respect to Brownian motion, that (j>{x, t, y) ^ [I — A{x, t, y)) [x + b{x, t, y)) — x for all 
X G Dy, t G [0, ty], y G R™ and all v € N and that 

limsup - • sup E 



limsup — • sup ] 

t\,0 \t x<£D^ 



iI-A{x,t,Wt))-'A{x,t,Wt) =0, (263) 



\A{x,t,Wt)x + b{x,t,Wt)f] ] <oo (264) 



for all w G N. Then 0: R'' x [0,r] x R'" — > R'' is {fi, a) -consistent with respect to Brownian motion. 
Proof of Lemma [K3d[ Note that 

(l>{x, t, y) - {A{x, t,y)x + b{x, t, y)) 

= ((/ - Aix, t, y)y^ - /) {x + b{x, t, y)) ~ A{x, t, y) x 

= {I- A{x, t, y))-^ A{x, t, y) {x + b{x, t, y)) - A{x, t, y) x (265) 

= ((/ - A{x, t, y)y^ - /) A{x, t,y)x + {I- A{x, t, y)y^ A{x, t, y) b{x, t, y) 

= {I- A{x, t, y)y^ A{x, t, y) {A{x, t,y)x + b{x, t, y)) 

for all X G Dy, t G [0,ty], y G R™ and all w G N. Combining equation (j265p . Holder's inequality, equation (|263L 
inequality (j264p and Lemma 13.251 then complete the proof of Lemma 13.301 D 

4 Examples of SDEs 

In this section, we apply the strong convergence results of Section |3] to a selection of examples of SDEs with 
non-globally Lipschitz continuous coefficients. In Subsection 14.11 we first describe the general setting in which 
these examples appear and then treat each example separately in the subsequent sections. 

4.1 Setting and assumptions 

The following setting is used throughout Section H) Let T G (0, oo), d,m G N, let {il,T,¥) be a probability 
space with a normal filtration {Tt)te[Q,T] and let W = {W^^\ . . . ,W'^"''>): [0, T] x 17 -J> R" be a standard 
(-7^t)te[o,T] -Brownian motion. Moreover, let Z? C R"* be an open set, let fj.: R"* -^ R'' and a: R'^ -^ R'^'"" be 
Borel measurable functions such that /i|_D : D — > R'' and a\D- D — > W^^™ are locally Lipschitz continuous and 
let X = {X'-^\ . . . ,X^'^'>): [0,r] X il — > I? be an (J't)tg[o j-j-adapted stochastic process with continuous sample 
paths satisfying E[||Xop] < oo for all p G [1, oo) and 

Xt = Xo+ f KXs) ds + [ c7{Xs) dWs (266) 

Jo Jo 
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P-a.s. for all t E [0,T]. Below it is often not necessary that all absolute moments of the initial random variable 
Xq of the SDE (|266p are finite but for simplicity of presentation this assumption is used in the sequel. 

Our goal is then to approximate the solution process X : [0,T] x fl ^ D of the SDE (|266p in the strong 
sense. For this we concentrate for simplicity on the increment-tamed Euler-Maruyama approximations in Sub- 
section [3A3l More precisely, let Y^ : [0,T] x fi — > R"^, N G N, be a sequence of stochastic processes defined 
through Y(^ := Xo and 

^ ^^ ' max(l,f||A.(y^/^)^+a(r^/^)(W^(„+i)T/iv-l^„T/iv)||) 

for all t e (^, N ] ' " ^ ■f'^' ^, ■ ■ ■ ,N — 1} and all N E N. In the next sections, we present several examples 
from the literature for the SDE (I266P and we show that the numerical approximation processes (|267l) converge 
strongly to the solution process of each of the following examples. To the best of our knowledge, this is the 



first result in the literature that proves strong convergence for the stochastic van der Pol oscillator (|269l) . for 
the stochastic Duffing- van der Pol oscillator (|272p . for the stochastic Lorenz equation (|275l) . for the stochastic 
Brusselator (|277p . for the stochastic SIR model (|283|) . for the SDE (|290|) from experimental psychology, for the 
Lotka-Volterra predator-prey model (|294l) and in the case a > /3^ (see Subsection 14.10.31 below for details) also 
for the instantaneous variance process (I309p in the Lewis stochastic volatility model. 

4.2 Stochastic van der Pol oscillator 

The van der Pol oscillator was proposed to describe stable oscillation; see van der Pol [751 and the references 
therein. As, for instance, in Timmer et. al |72| . we consider a stochastic version with additive noise (see also 
Leung [35] for more general stochastic versions of the van der Pol oscillator and also equation (4.1) in Schurz [57] 
for a generalized stochastic van der Pol oscillator with multiplicative noise). More formally, assume that the 
setting in Section |4T1 is fulfilled, let a, /3, 7, 5 G (0, 00) be real numbers and suppose that d — 2, m = 1, D ^M."^ 
and 

K S ) = ( a(7-(x0')x,-<5., ) ' K S ) = ( ? ) ^^''^ 

for all X — {xi, X2) G K^. Then the SDE (|266p is the stochastic van der Pol equation 

(i)^2W(2)_.^(i) ,. ^ ,... (269) 



dXp) = fa (7 - {X't'')'' ) Xi""' - 5Xl'> dt + l3 dWt 



for t e [0,00). In an abbreviated form, the SDE (I269|) can also be written as 

Xt-a{l-iXtf)Xt + SXt = l3Wt (270) 

for t G [0, 00). The diffusion coefficient a: M'^ ^- M^ in (|269p is globally Lipschitz continuous. The drift coefficient 
^ : R^ — )■ M'^ in (|269p is not globally Lipschitz continuous and also not globally one-sided Lipschitz continuous. 
However, the drift coefficient in (j269p satisfies the global one-sided linear growth condition 

(x,/i(x))<c(l + ||a;||2) (271) 

for all a; e M'' and some c G [0,oo). Therefore, Corollary [3T7| applies here with the Lyapunov-type function 
V: M^ — > [1,00) given by V{x) = 1 + \\x\\^ for all x ~ (a::i,X2) € R^ and we obtain limAr_>.oo sup^grQ^^^i E[||Xt — 
Yf\\i] =Ofor allge (0,oo). 

4.3 Stochastic DufRng-van der Pol oscillator 

The Duffing equation is a further model for an oscillator. The Duffing-van der Pol equation unifies both the 
Duffing equation and the van der Pol equation and has been used, e.g., in certain aeroelasticity problems; see 
Holmes & Rand [3T| . We also refer to Section 2 in Arnold, Namachchivaya & Schenk-Hoppe T for more details on 
the physical background of this equation. As, for instance, in Schenk-Hoppe [64j . we consider a stochastic version 
with an affine-linear noise term (see also Section 9.4 in Arnold [3] and Section 13.1 in Kloeden & Platen 43 ). 
More precisely, assume that the setting in Section ITT] is fulfilled, let {ai,a2,ct3), (/3i,/32, /^s) E M^ be two triples 
of real numbers with a^ > and suppose that d = 2, m ~ 3, D — M.^ and 

Xl \ f X2 \fxi\fOOO 



X2 ) V "1^1 ~ '^2a;2 - a'iX2 {xi) - {xi) ) ' " \ X2 J \ /SiXi l32X2 lis, 
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for all X = {xi, X2) G K^- Then the SDE (I266P is the stochastic DufBng-van der Pol equation 



(1) _ Tr(2) 






xr> dt, 
^(1) 



for t e [0, 00). In an abbreviated form, the SDE (I272p can also be written as 



.(1)n3 



di 



Xt - aiXt + a2-'i^t + a^XtiXt 



{Xt 



PiXtM^^ + p2XtM'^'' + /^sM^f ^ 



(272) 



(273) 



for t £ [0,00). The diffusion coefficient <t in (J272I) is globally Lipschitz continuous. The drift coefficient /i in 



(j272p is not globally one-sided Lipschitz continuous and also fails to satisfy the global one-sided linear growth 
condition (|27ip . However, Corollarv 13 . 1 71 applies here with the Lyapunov-type function F : R^ — >• [1, 00) given by 
y(a;i,X2) = l-l-(a;i) -|-2 (0:2) for all a; = (a;i,a;2) £ M^ and we hence obtain limAr_).oo sup^g [0^7-] E[||X4—yf^||'^] = 
for all q € (0, 00). 



4.4 Stochastic Lorenz equation 

Lorenz [50j derived a three-dimensional system as a simplified model of convection rolls in the atmosphere 
and this equation became famous for its chaotic behaviour. As, for instance, in Schmalfufi [53], we consider 
a stochastic version hereof with multiplicative noise. Assume that the setting in Section 14.11 is fulfilled, let 
(ofi, a2, as), (/3i, /32, /^a) G R'^ be two triples of real numbers and suppose that d = m = 3, D = M.^ and 




ai {x2 ~ xi) 

= I a2Xi - X2-' X1X3 
X\X2 - OizXz 



.1 \ 


/ /?i:ei 













hX2 





.3/ 


V 





^■iX-i 



(274) 



for all X ~ [x\,X2,X'i) E M.^. Under these assumptions, the SDE (|266p is thus the stochastic Lorenz equation 






aiXf ^ - aiX 



(1) 



dt + |3lX^^^Uw^^^\ 



a2X, 



(1) 



X 



(2) 



X,«Xf' 



dX 



(3) 



x,«xf) 



aaX, 



(3) 



di + l32Xf^ dWp'^ 

(3) ^w(3) 



(275) 



dt + P^Xl"' dWt 



for t £ [0,00). The diffusion coefficient in (|275p is globally Lipschitz continuous. The drift coefficient in (|275p 
is not globally one-sided Lipschitz continuous but fulfills the global one-sided linear growth condition (|27ip . 
Therefore, Corollarv 13 . 1 71 applies here with y: R^ _-, [l,oo) given by V{x) = 1 + j|xjp for all a; e E^ to obtain 
limjv^oo supt^[o,T] E[ll^t - Yt'^W"] = for aU q e (0, 00). 



4.5 Stochastic Brusselator 

The Brusselator is a model for a trimolecular chemical reaction and has been studied in Prigogine & Lefever [5T] 
and by other scientists from Brussels (Tyson [73] proposed the name 'Brusselator' to honour the pioneering 
work of this group). The following stochastic version hereof has been proposed by Dawson [M] (see also 
Scheutzow [SS]). Assume that the setting in Section [4.11 is fulfilled, let a,S E (0, 00) be two real num- 
bers, let gi,g2'- [0, 00) — >■ R. be two globally Lipschitz continuous functions with 51(0) = .92(0) = and 
sup3.g[o_3c-) |52(a;)| < 00 and suppose that d = m = 2, D = (0, 00)^, fi{x) = and a{x) = for all x E D'^ 

and 

xi \ ^ /" 5~ {a + l)xi+X2-ixif \ ^1' xi \ ^ I' gi{xi) 

X2 J \ axi-X2-{xif J' V ^2 / V ^ 52(2:2) 

for all X — {xi, X2) E D. Then the SDE (|266p is the stochastic Brusselator equation 



M 



(276) 



dX,« = 



'5-{a + 1)X['^ + Xi'\xi'^f] dt + 9i{xi'^)dW,^'^ 



dX. 



(2) 



^(1 



axr-xr(x, 



(2)^v(l)^2 



dt- 



-g2{X^^^)dW. 



(2) 



(277) 



for t E [0, 00). The diffusion coefficient in (j277p is globally Lipschitz continuous. The drift coefficient in (|277p is 
not globally one-sided Lipschitz continuous and also fails to satisfy the one-sided linear growth condition (|27ip . 
Here the squared sum of the coordinates is a Lyapunov-type function on the state space (0,c»)^. However, in 
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order to apply Corollary 13.171 in this example, we need to construct a Lyapunov-type function V^: R^ — > [l,oo) 
on the space M^. For this let </>: M — ?► [0, 1] and ip: M.'^ ^ [0, 1] be two infinitely often differentiable functions with 
4>{x) — for all X G (— oo, 0], with (l>{x) — 1 for all x G [1, oo) and with ip{xi, X2) = (l){xi) ■ (j){—X2) + (j>{—xi) ■ (j>{x2) 
for all x — (a;i,X2) G M^. Then consider V: M'^ — > [l,oo) given by 

5 2 

y(xi,X2) ^ - + {xi + X2) - 2 -xi • a;2 • -0(a;i,a;2) 

5 2 

2 (278) 

(Xi +X2j 

5.1,, ,,2 

x\ 



2 , 2 

||.x|| +2a;ia;2 (1 -'(/'(2:i,a;2)) 



a;iX2 > 

(0:1X2 < 0) A (min(|a;i|, |a::2|) > 1) 
otherwise 



for all X = {xi, X2) G K^ and observe that V^: M^ — > [1, 00) is infinitely often differentiable and satisfies 



xei 



ELillv-'-'c 



(1+lb 



sup -" a.iiJr ^'" "' < 00. (279) 



Next note that the fact 2 + a — 2a > ^ for all a e [0, 00) implies 



V{xuX2) = f + ||a;||^ + 2x1X2 (1 - tp{xi,X2)) 

Jl + ||x||' : (xiX2 > 0) V (min(|xi|, |x2|) > 1) 

- 1 5 , ||^||2 



x|| -2|xi||x2| : (xiX2 <0) A(min(|xi|,|x2|) < 1) 
>| + ||x||'-2max(|xi|,|x2|) (280) 



, 2 
2 



i + |min(|xi|,|x2|)P 



2 + |max(|xi|, |x2|)| - 2max(|xi|, |x2|) 



^ 1 + |min(|xi|,| x2|)|^ + |max(|xi|, |x2|)|^ _ 1 



2 2 

for all X = (xi, X2) G M^. In addition, observe that ||^'(x)it(x)||^,][j2 r) = for all x € I?^ and 

|r(x)a(x)||^(jj,^jjA /(xi+X2)(|5i(xi)| + |52(a;i)|) 



sup ^\;/,^"'"^ < 4 • sup ^^ ' ^\, ^ , ,' '"^^^^'^ < 00 (281) 

for all X = (xi,X2) G Z?. Combining (J279p - (l28ip shows that Corollary 13.171 applies here with the function 
V: M^ -^ [1,00) given in ((278)) and we therefore get limjv_^oo suptgfg yj E[||Xt - Y/'W'^] = for aU q e (0, 00). 
Please also note that other choices than (I278P are possible for the Lyapunov-type function y : M'^ — > [1, 00) and 
that the choice (|278p simply ensures that the smoothness and growth assumptions of Corollarv l3.17l are met. 

4.6 Stochastic SIR model 

The SIR model from epidemiology for the total number of susceptibles, infected and recovered individuals has 
been introduced by Anderson & May [2]. The following stochastic version has been studied first by Tornatore, 
Buccellato & Vetro [73). Assume that the setting in Section HTTj is fulfilled, let a, (3,^,6 G (0, 00) be real numbers 
and suppose that d = 3, m = I, D = (0, oo)'^, /i(x) = <7{x) — for all x G D'^ and 

-axiX2 - 5xi + (5 \ / ^1 \ / -/3xiX2 \ 

axiX2 - (7 + S)x2 , ''^ ^2 = /3xiX2 (282) 

7X2 - (5X3 / \ ^3 / \ / 

for all X = (xi,X2;X3) e D. Under these assumptions, the solution process {St,It,Rt) '■= i^t ' ^t ^-^t )' 
t G [0, 00), of the SDE dM]) fulfills 

dSt^[- aSJt - SSt + S]dt~ pSJt dWt 

dit = [aSJt - (7 + 5)It] dt + pSJt dWt (283) 

dRt^ [lIt-5Rt]dt 

for t G [0, c»). In this example, both the drift coefficient and the diffusion coefficient grow superlinearly. 
Moreover, the drift coefficient also fails to satisfy the one-sided linear growth condition (|27ip . We now construct 
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an appropriate Lyapunov-type function V^: M"^ — > [1, cx)) for this example so that Coroharv 13 . 1 71 can be apphed. 
Let (/): E — > [0,1] and tp-.M.'^-^ [0,1] be two infinitely often differentiable functions with (p{x) = for all 
X € (— cxD,0], with (f>{x) — 1 for all x € [1,cxd) and with 'ip{xi,X2) — (t>ixi) ■ (f>{—X2) + (f>{~xi) ■ 4>{x2) for all 
X = (xi,a;2) G M^. Then consider T^: M'^ — > [l,oo) given by 



V{xi,X2,X3) = - + {X1+X2) + (Xs) ~2- Xi- X2-1p{xi,X2) 



(284) 



for all X — (a;i, X2, x-^) e M^ and observe that V : 



[1,00) is infinitely often differentiable and satisfies 



sup 



i:LiIK'"(^)IL(.)(,,3,, 

(i+Nii) 



< 00. 



(285) 



Next note that the fact 2 + a^ — 2a > ^ for all a G [0, 00) implies 

V{xi,X2,xz) = f + ||a;||^ + 2xiX2 (1 - '0(a;i, 2:2)) 



Ml 



2|X1||X2| 



(a;iX2 >0)V(min(|a;i|,|a;2|) > 1) 
(xiX2 < 0) A (niin(|a;i|, |a;2|) < 1) 



> 



> 



■ + ||a;|| - 2max(|2:i|, |x2|) 
5 + |a;3|^ + |min(|a;i|, 1x21)1^ 
l + |a;3|' + |min(|xi|,|x2|)|'- 



2 + |max(|a;i|, |x2 
2 



|max(,|2;i|, \X2\ 



1 



2max(|a;i|, |a;2| 
.112 



(286) 



2 2 

for all X = {xi,X2,xz) € K'^. In addition, observe that |y(x)(T(a;)| = for all x £ D'^ and 

\V'{x)a{x)\ = -{-^V){X1,X2,X3) ■ /3 ■X1-X2+ {-^V){xi,X2,X3) ■ 13 ■ X1-X2 

= /3|xi| |X2| \{^V)iXi,X2,X3) - {^V)iXi,X2,X3)\ (287) 

= /3 |xi| 1x21 |2(a;i + X2) - 2{xi + X2)\ = 

for all X — (xi, X2, X3) e D and we thus get 

\V'{x)a{x)\^0 (288) 

for all X e M^. Combining (|285p - (|288p shows that Corollary [XTZl applies here with the function F : R^ -J> [1, 00) 
given in (|284p and we therefore get limjv_i.oo supjgrg j-i E[||Xf — y^^H^] — for all g S (0, 00). Please also note that 



other choices than (|284l) are possible for the Lyapunov-type function V: R — > [l,oo) and that the choice (j284p 
simply ensures that the smoothness and growth assumptions of Corollarv 13 . 1 71 are met. 



4.7 Experimental psychology model 

The motivation for the following example from experimental psychology is explained in Section 7.2 in Kloeden 
& Platen [43 . Assume that the setting in Section HTTI is fulfilled, let a,6 G (0,oo), /3 e R be real numbers and 
suppose that d = 2, m = 1, Z? = R^ and 



fJ- 



Xl 
X2 



13^X2 



(.X2) {5 + iaxi 
-xiX2{S + 4:axi) + ^-f^ 



Xl 
X2 



-13X2 
(3X2 



(289) 



for all X — (a;i,X2) G R^. The SDE (|266p is then the Stratonovich stochastic differential equation 



dX, 



(1) 



(xP)^(J + 4aX«) 



dt - px[^^ o dWt 



dX, 



(2) 



^(1) 



Xl'^Xl^^ {5 + ^aXY^) dt + PX\^^ o dWt 



(2) 



(290) 



for t G [0,00). The diffusion coefficient in p90p is globally Lipschitz continuous. The drift coefficient in (|290p is 
not globally one-sided Lipschitz continuous but fulfills the global one-sided linear growth bound (I27ip . Therefore, 
Corollarv 13.171 applies here with y: R^ — >■ [l,c») given by V{x) = 1 + ||x|j for all x E R^ and we hence get 
limjv^oo suptg[o,T] E[||Xt - YfWi] ^ for all q e (0,oo). 
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4.8 Stochastic Ginzburg-Landau equation 

The Ginzburg-Landau equation is from the theory of superconductivity and has been introduced by Ginzburg 
& Landau [20] to describe a phase transition. As, for instance, in Kloeden & Platen [33], we consider its 
stochastic version with muhiphcative noise. More precisely, assume that the setting in Section 14.11 is fulfilled, 
let a, /3, 5 € (0, oo) be real numbers and suppose that d = m ^ 1, D ~ R, ^(x) = ax — 6x^ and (j{x) — fix for 
all a; G R. Under these assumptions, the SDE (|266p is the stochastic Ginzburg-Landau equation 

dXt = [aXt - 5Xf] dt + pXt dWt (291) 

for t e [0, (X)). Corollary 13.171 applies here with F: R — > [l,oo) given by V{x) — 1 + x"^ for all x e K to 
obtain limjv_i.oo supj^rp yi E[|Xt — ^^^1''] — for all q S (0, oo). This example has a globally one-sided Lipschitz 
continuous drift coefficient and a globally Lipschitz continuous diffusion coefRcient. So, the convergence results, 
e.g., in [211 1301 [SI1I331 HH] apply here (see Section [1] for more details). 

4.9 Stochastic Lotka-Volterra equations 

The Lotka-Volterra predator-prey model and the Lotka-Volterra model for competing species have a quadratic 
drift term. Here we study the following more general Lotka-Volterra model as considered, for instance, in 
Section 7.1 in Kloeden & Platen [33] (see also Dobrinevski & Frey [E]). Assume that the setting in Section WA\ 
is fulfilled, let A, c = (ci,...,Cd) G K'^, v = {vi,...,Vd) G (0, oo)'', B G W^^"^ and suppose that d = m, 
D = (0, oo)"^, //(cc) = (j{x) = for ah x e D" and that 

M \ 

\ Bx) <0 (292) 

\ V ' ^d J I 

for all a; G M'' and 

xi \ I xi \ ( ^^ \ ( ^^^^ ^ 

Ml : = Ua + Bx), "{ ■ ]^\ "■ ' ^^^^^ 

Xd J \ Xd ) \ Xd ) \ CdXd 

for all X — {xi, . . . , Xd) G D. The SDE (|266p is then the d-dimcnsional stochastic Lotka-Volterra system 

/X« \ /c,X« \ 

dXt^ [ • . \iA + BXt) dt+\ ■■ \ dWt (294) 

V ■ x^^^V V ■ cx^'^V 

for t G [0, oo). Note that the drift coefficient of the SDE (|294[) contains a quadratic term and is therefore 
not globally Lipschitz continuous. Moreover, it depends on the matrix B whether the drift coefficient of the 
SDE (|294D is globally one-sided Lipschitz continuous. In order to apply Corollarv l3.17l a Lyapunov-type function 

V": M'' ^- [l,oo) needs to be constructed. For this let 0: M ^- [0,1] and V: M^ ^ [0,1] be two infinitely 
often differentiable functions with (f>{x) — for all x G (— oo,0], with (f>{x) = 1 for all x G [l,oo) and with 
'ip{x, y) — (j){x) ■ (f>{—y) + 4>{—x) ■ 4>{y) for all a;, y G M. Note that iIj{x, y) = -0(2/, x) for all x, y G M and '>j^{x, y) = 
for all X, 2/ G M with x • y > 0. Then consider V : M."^ -^ [1, oo) given by 




F(x) = 1 + |(u, x)|^ - ^ Vi -Vj ■ Xi ■ Xj ■ il;{xi,Xj) 

(295) 



:,je{i,2,...,d} 



= I +^\viXi\'^ + ^ Vi-Vj-Xi-Xj-{l-il;{xi,Xj)) 

i=l i,je{l,2,...,d} 

for all X = (xi, . . . , Xd) G M'' and observe that V : M'' ^ [1, oo) is infinitely often differentiable and satisfies 



-P { ^^"^l^''^^ ) < oo. (296) 
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In the next step observe that 

V{x) > 1 + mm(wi, 

> 1 + min(wi, 

> 1 + min(wi, 
for all a; G M** and this shows 

In addition, observe that 
Next inequality (I292p gives 



>«d; \m\ 



Y^ l[04)(min(|a;i|, Ixj])) \v^\ \vj\ \x^\ \xj\ 



t.je{1.2,-,d} 
i+2 



,Vd)\x\ - |max(t;i,...,Wd)| 



^ max(|a;,|,|xj|) 
,ije{i,2.....d} 



(297) 



2 J2 



,Wd)||a;|| - |max(wi,...,Wd)| d \\A 



^GR'' I ^{^) 



< oo. 




V'ix)a{x)\\ 



L(R<',1 



< OO. 



sup 

xeR'' 



1 + 11^11 
V{x) 



< oo 



(298) 



(299) 



(300) 



Combining ipM]) . (EMI), (PM)) and ([50111) shows that CoroUarv lXTTl applies here with the function F: R'^ -J> [1, oo) 
given in p95l) and we therefore get liniAr_j.oo sup^^jQ ^.j E[||Xt — Ft^lj*] = for ah q G (0, oo). Finally, let us 
describe two more specific examples of the stochastic Lotka-Volterra system (|294p . 



4.9.1 Stochastic Verhulst equation 



In addition to the assumptions above, let in this subsection 77, A G (0, 00) be real numbers and suppose that 
d — m — 1, A — 7] + ^ and B = —X. The stochastic Lotka-Volterra system (j294p thus simplifies to the 
one-dimensional Stratonovich SDE 



dXt^ 



irXt-x- {XtY 



dt + c- Xto dWt 



(301) 



for t G [0,00). Equation pOip is known as stochastic Verhulst equation in the literature (see, e.g.. Section 4.4 
in Kloeden & Platen [43 ). Note also that (|292l) is fulfilled here with v ^ 1, for instance. This example has a 
globally one-sided Lipschitz continuous drift coefficient and a globally Lipschitz continuous diffusion coefficient. 
Hence, the convergence results, e.g., in [351 [301 EZl IMl H] can be applied here (see Section [H for more details). 

4.9.2 Predator-prey model 

In addition to the assumptions above, let in this subsection a, /?, 7, 5 G (0, 00) be real numbers and suppose that 
d = 771 = 2, A = (a, —6) and 

• -/3 
7 



B = 



(302) 



The stochastic Lotka-Volterra system (|294p is then the two-dimensional SDE 

dX^/^ ^ X^'^ {a~(3- XP) dt + c- X['^ dwi^^ 
dXp) = Xp) (7 • X['^ -S)dt + C2- Xf ) dwP 



(303) 



for t G [0, 00). Note that (|292p is fulfilled here with v = (7,/3), for instance. The deterministic case (ci = C2 — 0) 
of this model has been introduced by Lotka [51] and Volterra [76] . 

4.10 Volatility processes 

There are a number of models in the literature on computational finance which generalize the Black-Scholes 
model with a stochastic volatility process. To unify some squared volatility processes of these models, we 
consider the following SDE. Assume that the setting in Section HT] is fulfilled, let a,a,P,6 G (0, 00), b G [5,00) 



54 



be real numbers with a+1 > 26 and S > l{i/2}(&) ■ 2 ^^'^ suppose that d = m = I, D = (0,cx)), /i(.x) ~ 5 — a-a;° 
and a{x) = (3 ■ x^ for ah a; € (0, 00) and ^(x) — a{x) — for ah x G (—00, 0]. The SDE p66p then reads as 

dXt = [(5 - a • [XtT] dt + l3- (Xt)' dWt (304) 

for i e [0,00). The assumptfon 6 > l^i/2}{b) ■ ^ ensures the existence of an up to indistinguishabihty unique 
strictly positive solution of p04p . In the sequel, the application of Corollarv l3.18l for the SDE (|304p is illustrated. 
For this define po G (l, 00] and qo € (—^,00] by 



00 



: b£ (0,1]U(0,^) 



PO ■= i 2a+0^ ., • (305) 



^°'tf : otherwise 



and 

90 



Po 1 

4{b-2+max(a,2h-l,3/2)} 2 

00 



Po€(l,3) 

Poe[3,oo) (306) 

Po = oo 



and note that this definition ensures po > qo- Next observe for every t e (0, cx)) that E[||Xt||P] < oo if and only 
if p < Po- In addition, note that Corollary 13.181 implies limAr_>.oo sup(g[Q_3-] E[||Xt — Y't^H^] = for all q G (0,oo) 
with q < qQ. Let us illustrate this by three more specific examples. 

4.10.1 Cox-Ingersoll-Ross process 

In addition to the assumptions above, suppose that a = 1 and b — ^. The SDE p04p is then the Cox-Ingersoll- 
Ross process 

dXt = [S- aXt] dt + l3y^tdWt (307) 

for t g [0, oo) which has been introduced in Cox, IngersoU & Ross [TU] as model for instantaneous interest 
rates. Later Heston [26] proposed this process as a model for the squared volatility in a Black-Scholes type 
market model. Here we have po ^ qo = oo and hence get limjv->.oo supj^rQ yi IE[||^t — Y^^H''] = for all 
q € (0, oo). Both the drift and the diffusion coefficient grow at most linearly. Therefore, strong convergence of 
the Euler-Maruyama approximations is well-known (see, e.g., Krylov [45] and Gyongy [21) for convergence in 
probability and pathwise convergence respectively). Strong convergence rates of a drift-implicit Euler method 
for the SDE p07p are established in Theorem 1.1 in Dereich, Neuenkirch & Szpruch 1161. 

4.10.2 Simplified Ait-Sahalia interest rate model 

In addition to the assumptions above, suppose that a ~ 2 and 6 < |. Under these additional assumptions, the 
SDE (pOi)) reads as 

dXt=^[S-a- {Xtf ]dt + l3- {Xtf dWt (308) 

for t G [0, oo). A more general version hereof has been used in Ait-Sahalia [P for testing continuous-time models 
of the spot interest rate. Here we also have po = 9o = oo and therefore liniAr^oo supjg[o,T] lE[ll-''^t ~ ^t^ll'^] = 
for all q G (0, oo). A strong convergence result for this SDE is Theorem 6.2 in Szpruch et al. [7T| . 

4.10.3 Volatility process in Lewis stochastic volatility model 

In addition to the assumptions above, suppose that a — 2 and 6 = |. The SDE (I304p is then the instantaneous 
variance process in the Lewis stochastic volatility model (see Lewis [TT] ) 

dXt=[5-a- {Xtf ]dt + P- {Xt)i dWt (309) 

for t G [0, oo). Here we get po = ^"t/ G (1, oo) and go = max(^ ~ 5' O) = inax(^^, O) G [0, oo). In the case 

a > /3^ we thus get limAr_j.oo supjg^o.T] ^[ll^t ~ ^t^y] = for all q G (0,(7o)- The stochastic volatility model 
associated to p09p is a.k.a. 3/2-stochastic volatility model (see also [551 [55]). Note that Lemma [2.291 can be 
applied here to prove moment bounds of a linear drift- implicit Euler method for the SDE (|309p . 
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4.11 Langevin equation 

A commonly used m.odel for the motfon of molecules in a potential is the Langevin equation (see, e.g., Section 2.1 
in Beskos & Stuart 6 for more details). Corollarv l3 . 1 7l does not apply to completely arbitrary Langevin equations 
but requires the following assumptions on the potential. Assume that the setting in Section |4J] is fulfilled, let 
e, c e (0,oo) be real numbers, let U : W^ — >• IR be a continuously differentiable function and suppose that 
£) = K'^, d = TO, ijl{x) = -{vU){x) and a{x) = \/2e/ for all x G W^. Moreover, let F: R'^ ^ [l,cx3) be a twice 
differentiable function with a locally Lipschitz continuous second derivative, with limsup„x q sup^-gRd yr^ < oo 

and with {s7U{x),VV{x)) > -c • V{x) and X^Li \\V^'Hx)\\L(OiR''M) < c\V{x)\^^-'^/''^ for Agd-almost ah x e R'^. 
Under these assumptions, the SDE (I266p reads as 

dXt =: -(vC/)(Xf) dt + V2edWt (310) 

for t e [0, oo). Corollarv 13.171 then implies hmAr^oo supjg^o.T] ll^[ll-''^t ~ ^t^y] = for all q £ (0, cx)). If the force 
(— Vt/) in the Langevin equation is globally one-sided Lipschitz continuous and satisfies suitable growth and 
regularity conditions, then the strong convergence results, e.g., in [311 [301 IHZl IMl HH] apply here (see Section [T] 
for more details). 
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